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Some Results of Sweptback Wing Structural 
Studies 


A. L. LANG* anp R. L. BISPLINGHOFF7 
Massachusetts Institute of Technology 


SUMMARY 


Static stress and deflection tests have been conducted in the 
laboratory on a 45° sweptback single-cell beam. The results of 
these tests are compared with simplified theories for predicting 
stresses and deformations in swept shell beams. The agreement 
between theoretically and experimentally determined deflections 
is closer than the stresses, but both agree within acceptable engi- 
neering accuracy. Both cases of skewed ribs and ribs perpendicu- 
lar to the spars are included. 
deformation is unavoidable in sweptback wing static tests, a 
method is described of accounting for this effect in the theory. 


Since some degree of root support 


NOTATION 


corner flange area (except as used in matrices), 
in.” 

wing chord perpendicular to spars, in. 

disturbance influence function due to the effect of 
sweep 

influence function associated with a straight wing 
clamped at the root 

Young’s modulus, Ibs. per sq.in. 

shear modulus (except as used in matrices), Ibs. 
per sq.in 

depth of cross section, in. 

cross-sectional moment of inertia, in.‘ (also de- 
notes unit matrix) 

length of straight wing, in. 

bending moment (vector perpendicular to spars), 
in.Ib. 

concentrated shear force, lb. 

torque (vector parallel to spars), in.Ib. 

cover thickness, in. 

web thickness, in. 

potential energy of external loads, in.Ib. 

internal strain energy, in.Ib. 

vertical displacement, in. 

beam slope, rad. 

torsion-bending parameter, in. 

sweep back angle, deg. 

warping angle of cross section, rad 


Presented at the Structures Session, Nineteenth Annual Meet- 
tg, 1.A.S., New York, January 29-February 1, 1951. 
» *Engineer, Aero-Elastic and Structures Research, M.I.T. 
| tAssociate Professor of Aeronautical Engineering, M.1.T. 


flange normal stress, lbs. per sq.in 
shear stress, lbs. per sq.in. 
twist of cross section, rad. 


(I) INTRODUCTION 


. es CALCULATION OF THE STRESS distribution in 
the vicinity of the root of a sweptback shell beam 
is an extremely difficult problem. Solutions to the 
exact equations are obtainable by numerical methods, 
but they are much too cumbersome for the designer. 
An approach is needed which is sufficiently comprehen- 
sive to yield the significant effects of sweep, yet simple 
enough to be usable by the structural designer. Ex- 
perimental data have shown that the effects of the 
sweep which are of greater significance to the structural 
engineer are the concentration of normal and shear 
stresses at the root of the rear spar, and the coupling in- 
duced by the sweep between bending and torsional de- 
flections. 

Little theoretical and experimental data have been 
published by the sweptback wing structural problem. 
Levy' described a method of computing stresses and de- 
flections in swept wings with ribs parallel to the air 
stream. This method consisted of breaking the struc- 
ture into its elements and determining the internal 
forces by solving the equilibrium equations plus ad- 
ditional equations determined from the principle of 
minimum strain energy, Thompson and Wittrick” * 
present a scheme based on the work of Hadji-Argyris for 
calculating the stress distribution in a swept tube with 
rigid ribs parallel to the air stream. Goodman‘ presents 
a computation of the stress distribution in a sweptback 
sheet based upon the methods of reference 1. Zender 
and Libove® and Zender and Heldenfels® give excellent 
and useful experimental data on the behavior of a 45° 
sweptback box beam subjected to bending and torsion. 
In their latest paper, Heldenfels, Zender, and Libove’ 
present methods for the analysis of a 45° swept box 
beam with ribs perpendicular to the spars and with an 
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elastic carry-through structure. This paper, which has § wi 
a = recently come to the attention of the authors, employs § flat 
a \ z methods and reaches conclusions similar to those pre. | Th 
sented in the present paper. These eight papers, out. § fla 

\ ; side of the M.I.T. Report® upon which the present paper § vet 
> 2 —— fe am is based, seem to comprise the main published theoretica] J ma 
and experimental literature on the swept shell beam § roc 
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SECTION A-A 


problem. 

The purpose of the present paper is to compare known § pre 
methods of analysis with experimental data and to pre. J effi 
sent simplified methods that can be applied to a single. J 2t 
cell shell beam. In the first part of the paper, laboratory } ing 
experiments on a 45° sweptback beam are described in yie 
which both stresses and deflections were measured, } fle 
The remainder of the paper is devoted to analysis of the C 
test specimen and comparison of theoretical and ex. 
perimental results. At first, a modification of the 
method proposed by Levy in reference 1 is applied, 











\ 
» Secondly, a method based upon the principle of mini- 
| ae mum potential energy is applied. Finally, a simplified 
| a 1 approach is suggested which can be used in preliminary 
| design. 
| | Th 
| ! i2v2 (II) Test SPECIMEN 
eve | The plan form and cross section of the sweptback ;, 
wing specimen used in the present investigation are 
sini illustrated in Fig. 1, and a photograph is shown in Fig. , 
a E.R ees pro. 2. The specimen is a stressed skin, single-cell box swept 
back at 45°. The two spars are parallel and of equal ( 
cross section along the span. Top and bottom surfaces 
are identical. Each spar flange is composed of two */, ( 
by */, by '/s 24-ST extruded angles back to back and 
separated by the spar webs. Three equally spaced '/; 
by '/2 by '/i§ 24-ST extruded angles are used as longi- snp 
tudinal stiffeners on the upper and lower skins. Six dia 
solid ribs parallel to the wind stream are spaced 6+/2in. opt 


along a line normal to the airplane centerline. The - 





webs and ribs are 0.051-in. 75 ST sheet, and the upper low 
and lower skins are 0.032-in. 75 ST sheet. The overall ligh 
length of the wing specimen measured parallel to the refl 
spars is 60 in. The section inboard of the inner rib is 
built around a heavy steel attachment fitting as shown 
by Fig. 3. 7 
SECC 
(III) Static Tests the 
The objectives of the static tests of the sweptback _ 
wing were to obtain experimental data for purposes of = 
comparison with theory. In the present paper, the ex- obt 
perimental determination of the following are included: cus 
(a) Nine by nine matrix of flexibility influence co- effe 
efficients. “th 
(b) Normal and shear stresses in the flanges and webs has 
due to a single load applied midway between Stations 2 —, 
and 3 (cf. Fig. 1). die 
Fig. 3 illustrates the arrangement employed in the I 
static tests. Vertical deflections were measured by dial plic 
asst 


gages attached to a fixed supporting structure above the 


Fic. 2. Test specimen. 
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RESULTS OF SWEPTBACK 
wing. Stresses were measured at 43 points on the 
fanges and webs by means of electric strain gages. 
The gages were mounted longitudinally along the 
fanges to measure the normal stresses and at 45° to the 
vertical on the webs to measure the shear stresses. The 
majority of the gages were located in the vicinity of the 
root. 

The wing was loaded by a single hydraulic jack. The 
procedure for determining the matrix of influence co- 
efficients consisted of applying a 1,000-Ib. load at points 
2 to 10 (ef. Fig. i) and recording the nine dial gage read- 
ings for each point of load application. This procedure 
yields 81 numbers, which form the square matrix of 
flexibility influence coefficients for the wing denoted by 


Clox 9: 
: & { 6P a] 6P \ 
€ 10K } 9 ’ Cor ee C 10R {> 


{fr 46P 
) Cor eee 


The elements of these matrices are defined as follows: 


Cr’ = front spar deflection at root due to unit load 
at Station 1 

Cy’ = rear spar deflection at root due to unit load 
at Station | 

Cy*” = front spar slope at root due to unit load at 
Station | 

ON or enar s ; : ‘ 

Ce“ = rear spar slope at root due to unit load at 


Station 1 


Linear root deformations were measured with 0.010-in. 
dial gages, and slopes were measured by means of an 
optical system. In the latter, a plane first surface 
mirror was attached to a bar pinned to the upper and 
lower flanges of the spar at the root. By means of a 
light beam, a hairline was projected onto the mirror and 
reflected to a screen. 
(IV) IDEALIZED SECTION 

The swept-wing problem entails an investigation of 
secondary normal and shear stresses in the vicinity of 
the root. Because of the complexity of the theory 
necessary to investigate such effects, it is common to 
simplify the cross section in order to reduce the labor in 
Such simplifications are dis- 
cussed, for example, in reference 9. In addition, the 
effective cross section will depend upon the type loading 
to which the structure is subjected. This effect, which 
has been neglected here, requires further investigation 
since it is not possible to separate bending and torsion 
effects near the root of a swept wing. 


obtaining a solution. 


In constructing simplified sections in the present ap- 
plication, the cover skin is assumed fully effective. This 
assumption seems reasonable due to the relatively low 


WING 


STRUCTURAL STUDIES 707 

In conducting static tests on sweptback wings at 
M.I.T., it has been found that it is difficult to design 
and construct a root fitting that provides anything ap- 
proximating complete root fixity. The fitting finally 
selected shown by Fig. 3 did not provide sufficient root 
fixity so that its elastic effects could be neglected com- 
pletely. Root fitting deformations have a significant 
effect on the wing deflections and a somewhat lesser 
effect on the stresses. Corrections for both effects 
should be introduced in the comparison between theory 
and experiment. In the present method of correction, 
experimentally determined data on root deflections are 
introduced into the theoretical computation of the 
This experimental data is 


stresses and deflections. 


represented by four column matrixes as follows: 


7 r ” P\ if P >. ori 
ee «ce Gee ac eee” 3c Oe I 


magnitude of the applied loads in the static tests. Fig. 4 
shows the actual wing section and the simplified section 
employed in the analyses. 


(V) THe Levy METHOD 


(a) General Remarks 


In general, for complex structures such as are some- 
times found in the root regions of sweptback wings, the 
general principles described by Levy’ are often the only 
tool that can be used to obtain a satisfactory solution. 
This is especially true when the yielding of elastic ribs 
forms an important part of the problem. Although the 
scheme used by Levy in setting up the problem of a 
swept wing with ribs parallel to the air stream leads to 
incorrect results, the case of ribs perpendicular to the 
spars can be readily handled. 





Arrangement for deflection and stress measurement of 
test specimen 


Fic. 3. 
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SIMPLIFIED FOUR 
FLANGE SECTION 


Simplification of cross section perpendicular to spars 


ACTUAL SECTION 


Fic. 4. 











Structural segments and loading for analysis of swept- 
back wing by the Levy method 


Fic. 5 


b) Application to a Swept Wing with Elastic Ribs 

Perpendicular to the Spars . 

As an illustration, the 45° sweptback wing described 
in Section (II) is analyzed as though its ribs are perpen- 
dicular to the spars. Fig. 5 illustrates schematically the 
assumed loading and the manner in which the elastic 
ribs divide the wings into segments. Using the methods 
of reference 1, the wing may be broken up into struc- 
tural units, and each unit can be put into equilibrium 
with internal forces. This results in the diagram shown 
in Fig. 6. In this diagram, each sheet and stringer is in 
equilibrium, and 15 unknown internal forces a, . . . is 
have been introduced. The stress distribution is com- 
pletely known when these 15 values are determined. 
By inspection of Fig. 6, it is clear that nine equations of 
equilibrium can be written as follows: 


iPs.. § Po} = [7] }a,.. . a} (1) 


The system has six redundants, and six additional equa- 
tions must be derived by applying the principle of 
minimum strain energy. 

In writing an expression for the total internal strain 
energy in the structure, three fundamental formulas are 
Two of these are given in reference | for the 
In computing the in- 


required. 
flanges and rectangular panels. 
ternal strain energy in the triangular panel, it is as- 
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sumed that the shear flow is uniform over its area. This 
approximation results in an expression for strain energy 


as follows: 
U; = 1.3(¢°A/Et) D) 


where A is the triangular area and ¢ is the skin thick. 
ness. 

By means of the three fundamental formulas, the 
strain energy can be expressed in terms of the internal 
forces in the following bilinear form, 


U, = (1 2E)l[a,... ais] F] $0, .. . ys} 3 
Since it is necessary to express the strain energy in terms 
of the external forces and the six redundant members, 
the following equation is constructed from Eq. (1 


{ ’ rev 
,2,..- Ay, = G ) Pe are Py, Q2, As, As, Ai, Ay4, Ays; 


Substituting Eq. (4) into Eq. (3), 


U = ] 2E P, a Py, Qe, As, Ag, Ay, A114, A15 x 
H ,P a Py. (eo, A5, Ag, Ay, Ai4, a5} 5 


where [/7| = [G]’|F)|G] and [G]’ is the transpose of 
G]. 
Applying the principle of minimum strain energy, 


fy Bey AS, 


é, V0, > 


0U,/Oa; = 0 5 = 14, 15 6 


the internal redundant forces are obtained, 


Hep] {P:... 


j ’ 
) 2, As, Ag, Air, A14, Gis, = o 


land |H,p 


| Hpp|9 x 9 
Hap|s x 9 
Substituting the result given by Eq. (7 
stress distribution can be seen to be 


are found by partition- 


[palo x | 
‘Pte 6 x 6 


into Eq. (4), the 


The matrices | /7,, 
ing | /7] as follows, 


fI| = 


‘ ‘ ’ ' 
12, ..-.455 = S| [P.... Pr} Y 


G | an x9 7 | 
"4 | —(H..)— [Hep] 


Substituting the result given by Eq. (7) into Eq. (5) and 


using the relation, 


where 


S| = 


U, = '/-[P2... Pwl(C] {P2... Pw} 10 
where [C] is the matrix of flexibility influence co- 
efficients, there results for |C] 

[C] = (1/E) [[Hee] — (Hea) [Hua)~ [Hor)] (1 


The 15 by 9 matrix for the internal load distribution 
and the 9 by 9 matrix of flexibility influence coefficients 
can be obtained from the partitioned members of the 
strain energy matrix [/7}. 
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WING 


STRUCTURAL STUDIES 


c) Corrections for Root Deformations Ry = [Tr] fai... ass! (16) 
Because of the impracticability of constructing a root Re = [Te] ja... ay} (17) 
fitting that completely suppresses root deformations, it Vr = [Sy] far... aw} (18) 
is necessary to include these effects when comparing Ve = [Se] far... ar} (19) 
theory with swept-wing experimental data. In the 
present case, this is done by introducing experimentally where 
determined Toot stiffnesses into the analysis carried out Rr = flange force in front spar at root 
above. This results in corrected expressions for Eqs. Ry = flange force in rear spar at root 
9) and : 1). . I» = vertical shear on front web at root 
[ sing the experimentally determined column matrices Ve = vertical shear on rear web at root 
from Section (III), the linear and angular deflections at 
the root of the front and rear spars are Try] = [2-8 0 2-6 0 2-4 0 2-2 O 2 
. ~ 8P . 6P 0-2] 
Sy = [P.... Pu] {Cap ... Cor™} (12) re = e-8 A ee ee 
i =o 1... Pal 1... Co”? (13) 2 a . iia 
_ [ P. P ] ee aP c aP | 14 
Oe if Dee Ww} ) Ger oes 10F ‘ ) + 7 a 
on = (Po... Pw] 1 Con”? .-- Cun?! (am [Se] = [0 000000000 0 0 
1 O 0} 
The root reactions can be expressed in terms of the in- [So] =(0 0000000000 1 
ternal forces, 0 1 O 
The total internal strain energy in the wing plus root fitting is given by 
U, = Ui, + (1/2)[har hag br be) {Re Re Vr Vr! (20) 
where l’; is given by Eq. (3). Introducing Eqs. (12) to (19) into Eq. (20) yields 
U = (] 2E)(P» $a Pr, dle, A5, As, Ayty A14y ay5|{H7] | P, 08 Po, Qe, Ap, As, Aity A14, a5} (21) 
where |/7] = [77] + [J%)] and [J4\] is given by 
i= Ele... Cor™.&66066 0] [Tr] +h 1G... Cos”. 66006 0) [Te] + 
Cor”... Cre”, 0, 0, 0, 0, 0, 0} [Sr] + {Cap . . . Cros”, 0, 0, 0,0, 0,0} [Sp] [G] 


The matrix |//,| has the following form, 


[Hpplo x 9: [Hpalo x 6 
[11] = | eee | 22 
[Hap |¢ x 9 [O}s x 6 


In general, {/7,] will not be symmetrical, but can be 


made symmetrical by putting 


Qi = ay = (Ay + Ay)/2 (23) 
where 
a, = an element of [/7,;] symmetrical 
A, = an element of {/7,] unsymmetrical 


Referring to the results given by Eqs. (9) and (11), it 
can be seen that the corrected stress and influence co- 
elicient matrices are given by 


(J 9x9 
io] = «il a = | (24) 
— [Hoa] Hop] 
C) = 1/E [(Ape] — (Apa) [aa] '[Aar]] (25) 
It should be observed that [H,,] = [Ha]! so that 


lurther matrix inversion is unnecessary when including 
toot deformations in the analysis. 


d) Comparison of Experiment and Theory 


The results of applying Eqs. (24) and (25) to the 45‘ 
sweptback wing are compared with experimental data 
from the test specimen. This comparison is not com- 
pletely justified since the rib direction assumed in the 
analysis is perpendicular to the spars, whereas the rib 
direction in the test specimen is parallel to the air 
stream. 


Figs. 7 and 8 compare flange normal stresses and web 
shear stresses for a unit shear force applied midway be- 
tween Stations 2 and 3. Fig. 9 compares front and 
rear spar deflections for a unit load applied at Station 2. 


In comparing these results, it can be seen that the 
agreement in deflections is somewhat better than the 
agreement in stresses. It is important to observe that 
the assumption of constant shear stress over each bay 
does not permit the theory to predict the critical root 
shear stress in the rear spar. On the other hand, the 
assumption of a linearly varying normal stress yields 
results that compare favorably with theory at the root 
of the rear spar. 
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(VI) THe Metnop oF MINIMUM POTENTIAL ENERGY 
(a) Application to a Swept Wing with Ribs Perpendicular 

to the Spars 

When it is assumed that all ribs are rigid except the 
root rib, the method of minimum potential energy can 
be conveniently applied. It can, of course, be applied 
also to the elastic rib case, but it is somewhat more 
tedious than the method of minimum strain energy 
(Levy method). 

Fig. 10 illustrates the loading that is assumed in the 


present application. The wing shown in Fig. 10 is 


U; = Ve. a, eee 


j 


U1.» 


The potential energy of the external loads in matrix 
form is given by 


U, = Py T; ene 74] fv, ee 0} (27) 


—([Pi... 


The total potential energy is given by, 


oe 


[K] }0 oe 
Inverting the matrix [K], Eq. (30) becomes, 


- U4, A; Ps G4, Pi--- 


‘ 
pUr .. - Ua; 0; sas 04, Gy... M4, Yi... G4, Vo, Ao, vo} = 


The matrix [C] can be partitioned as follows, 


» — [lCu]s xs | [Culs x u on 
ili far x0! (Caln x | ans 


where 
[Cie]s a [Ca] x 8 


The 8 by 8 matrix in the upper left-hand corner» 
[Cus x s, is of most interest in problems involving cal- 
culations of deformations of the wing and in problems 
in dynamics. The elements of this matrix represent the 
linear and angular displacements at the i‘ section due 


to unit force or torque at the j™ station. 


AERONAUTIC 


= ) { 7 
$4, Vo, 2, Yo; = Py ma ia P, T; ie 
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assumed to have nineteen degrees of freedom which 
completely define its strained configuration. These con. 
sist of a linear and two angular displacements at the 
juncture of the elastic rib with the front spar (Yo, a, and 
vo), and four displacements at each of the four loading 
stations. These four displacements are slope, angle of 
twist, warping, and vertical displacement denoted by 
a, ;, gi, and v;, respectively. Fig. 11 illustrates positive 
displacements. 

The internal strain energy in terms of the linear and 
angular displacements can be written in the following 


matrix form, 


« G4y Voy Ao, Yo] [K] 4 


Pay 9 ( Syn 
- U4, A; oa 04, Gy... - Ay, Pi. - + Oty VO, Ao, Yos (26) 


VU=U,+ U, (28 
Applying the principle of minimum potential energy, 
(29) 


there is obtained, 


.1,,00000000000; (30 


(C] }Pi...Ps T1...71,0000 0000000} (31) 


In determining the matrix [Cj |s y s, it is not necessary 
to invert the entire [A] matrix. It is observed that the 
[K] matrix can be partitioned as follows: 


[K] = eee s | [Kule x "| 


where 


[Ka |u ig [Ki2]s x un’ 

The desired matrix of flexibility influence coefficients, 
[Cu]s x s, can be written in terms of the partitioned ele- 
ments of the [A] matrix as follows, 
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[Culsxs = [[Kulsxs = [Kielsxu [Kee}u yu! [Kalu xsl “" (34) 
Some simplification in the analysis is obtained if it can be assumed that the wing is composed of an infinite num- 
ber of closely spaced rigid ribs outboard of Station 1 (cf. Fig. 12). Fig. 12 illustrates the wing loaded at a station 
located a distance x = ¢ from the most inboard rigid rib. The deformations of the wing at a distance x are denoted 
by v(x) and @(x), and at the first station by 1, a, 6, and g;. It is convenient to write the deformations at the 
point x as follows: 
1x00 10 : 

fo(v\A(v)t =< f \ {a(x\A(xv) t ar 

V(x)O(x) 5 = = v 6 d(x )O(: (3é 

01 = [Jot Ka | tromied + [95] foe 
where 3(x) and 6(x) denote the deformations at the point x when it is assumed that the wing is rigidly clamped at 
the first station. The term A(x) denotes the angle of twist at point x due to unit warping at Station 1, and is given 
by 

_ BAE (b h ; cosh Bx — 1 — 
K(x) = — _ (sinh fx — (36) 
8Gb \E. ty coth 6/ 
where 
B* = 8G/[AE(6/t. + h/tw) 
The internal strain energy can be written in two parts, 

U,; = |] /2[v, ay A, gi |[Ky] 10; al A, ¢i} + U a(x), 6(x), o(x) | (37) 
where [, is the internal strain energy associated with a(x), 0(x), d(x). The matrix [K,] is derived in Appendix 
(A). 

The potential energy of the external loads is, 
v1 + gay v() 
U. = —(PMT}| a — [PMT] &() (38) 
6+ K(é¢ 6(£) 
\pplying the principle of minimum potential energy, there results two separate sets of equations, 
[Ky] {rarPi¢r f = [A ] }PMT\ (39) 
where 
1 0 O 
oe E I g 
4) =|°* 
4I=!9 9 4 
0 0 K¢é) 
and 
x St 2>é 
di(x) l | x P da(é) 
—— = —| Mx + Pxié — + — a(x) = o(£) + ~(x — &) | 
dx El a a, wR wreg *-* | 
, 2T L (tub — teh , 1 /tob — teh (40) 
?) ———————- + ( ) a0) 6’(x) = ( were? ik. ( 
©) = Gonecd + th) * 6 Nth + ted) ™ ©) = lee + tb) OO 


4 tb + *) ; 
ao" (x = 2a(x — << ro 7 o"(x — 23(x = 0 | 
" (x) — BG(x) (= Lb e"(x) — B(x) 


AEbh* 
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The solution of Eq. (39) yields the deformation of the Substituting the results of Eqs. (41) and (42) into the 
system assuming that outboard of Station 1 the wing is deformation Eq. (35), there is obtained, 
completely rigid in bending and twisting but elastic in 


warping. The solutions of Eqs. (40) yield the elastic bo(x)O(x)} = [[Csr(x, £)] + [Cp(x, 2) ]] | PMT! 
deformation of the wing outboard of Station 1 when it is (43 
assumed that the wing is completely clamped at Station 
3 where 
The solution of Eq. (39) can be written as, 
ae i 
frig: = [R] {PMT} (41) art I= 1g yp ee 
where ' ane , . 
is the matrix of influence functions assuming that the 
[A,J—'[A] = [R] wing is rigidly clamped at Station 1. Each term repre 
sents the linear or angular deflection at x due to a unit 
The solutions of Eqs. (40) can be written as, force or moment at &. 
Sa(x)O(x)} = [C(x, £)] | PMT} (42) [C.(x, £)] = lL « 8 0 _ (as 


0 0 1 A(x) 
where [C(x, £)] represents a matrix of influence func- 
tions resulting from solutions of Eqs. (40), and is a is the matrix of influence functions assuming that the 
well-known result"’ for straight wings with warping re- wing is, rigid in bending and twisting but elastic in 
straint at the root. warping outboard of Station 1. 
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Fic. 6. Sweptback wing broken down into structural units each in equilibrium. 
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(VII) SIMPLIFIED SOLUTION BY THE METHOD OF aaa 


8 orimanpaticipimmanoamin : : a 
MINIMUM POTENTIAL ENERGY [ REAR SPAR ——— THEORY ASSUMING 
7 ELASTIC RIBS 






a) Application to a Swept Wing with Closely Spaced 6| 





EXPERIMENT | 
Se ae . . Mc ® FRONT SPAR x 
Rigid Ribs Perpendicular to the Spars 5| » ime L REAR SPAR © 
o | 
The method described in Section (VI) suggests a P g| al 
smewhat more simplified approach to the problem. " : 
The outer section of the wing is assumed to be a box FRONT SPAR 
with closely spaced rigid ribs and joined to the inner 2! 
part of the wing at Station 0 as shown by Fig. 14 in- I 
stead of at Station 1 as shown by Fig. 12. The dis- = satiiieilics _i | = oe ee 
lacement of the juncture of the inner an1 outer sectior oS ¢ SO ee ee 
a. J wale ne STATION ALONG SPAN (INCHES) 
is assumed to be described by % and a and the dis- Fic. 7. Normal stresses in flanges of sweptback wing result 


placements at a point x in the outer section by v(x) and — ing from unit shear force at center of Section 2-3 (obtained by 
= 2 - i - oo Levy method). 

4x). With this arrangement, the continuity conditions , 

are obviously not satisfied at Station 0; however, the 





essential features of the problem are preserved. 


























The f = : ¢ — a ae 4—-—— — a r - —_—_ 
The following transformation is introduced, THEORY ASSUMING ELASTIC RIBS 
1/2 x/2 10 EXPERIMENT: FRONT. SPAR x, REAR SPAR ® 
\ cs A/4 </\f 
10(x)O(x) 5 = k b Pll } Vor | + E Heenan) 3 } | { | 
@ REAR SPAR | = 
, ; ; } , H ; 
146) me. 
: ; ; 2 j | 
where #(x) and @(x) denote the deformations at a point =, a. 
rif it is assumed that the wing is rigidly clamped at P | 
° * ¢ . ° . “A a x | 
Station 0. The factor A(x) is given by Eq. (36). | — 
; ; ‘ ‘ iP . FRONT SPAR 
The internal strain energy is written in two parts, ; 
U; = (1 /2) [voao| [Ky] } Yoaty } + U [6(x) O(x) o(x) | O} F 
yd -0.5U—_1._1 1 1 1 ee ee ee ah — 
(47) 0 4 8B 12 16 20 24 28 32 36 40 44 48 52 56 60 


STATION ALONG SPAN (INCHES) 


and the potential energy of the external loads is, ac- : —" : é ‘ 
i Fie. 8. Shear stresses in webs of sweptback wing resulting from 


cordingly, unit shear force at center of Section 2-3 (obtained by Levy method ) 
(¥/2) + E(ao 2) vté) 

Ul. = — [PMT) ao/s — [PMT] | a(é) (48) 
(Up b) + K(é) ao 6(é) 


Applying the principle of minimum potential energy, two separate sets of equations are obtained, 


1/2 0 1/6 - 
/2 1/8 bal {PMT} (49) 


[Ay] | Toa} = | 
and the results given by Eqs. (40). 
The deformations at any point x, due to loads at point £, are given by 
| 2(x)0(x) } = [[Csr(x, &)] + [Cp(x, &)]] {PMT! (50) 


where [Csr(x, &)] is the matrix of influence functions assuming that the wing is unswept, and 


1/2 x/2 : 1/2 0 I1/b 
[C.(v. == zs 1 5 
LCr(, §)] k b Tl [Ay] } 2 1/2 i anes 


are correction influence functions to be superimposed upon the unswept wing results to account for the effects 
of sweepback. 

The stresses in the wing can also be regarded as a superposition of the stresses in an unswept wing plus correction 
terms accounting for the effects of sweep. Considering the important rear spar stresses, the normal stress in the 
tear flange at point x due to loads P, M, and 7 at point £ is given by 


COrs(X) = osr(x) + aof(x) (52) 


where 
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(x) k —x |] BE (° h )/ ‘ah 8 cosh =) ‘Pur} 7 
OsT = = —~ . sin \ ies i De 
, 2Ah 2Ah 8Ghb\t. tp)\" * ~ cosh BE . 


is the normal stress distribution in the rear spar assuming the wing clamped at x = 0, and 


op(x) = (hBE/4) [(cosh Bx/coth B€) — sinh Bx] [S,] {PMT} (54 


IK aR ‘2 0 1/b | 7 | 
“" Le/2 1/2 K(é) [Seliyes 


is the correction normal stress distribution for the effects of sweep. 


in which 


The rear spar shear stress can be expressed in a similar manner, 


Trs(X) = tTs7r(x) + Tr(x) (55) 
where 
(x) | l BAF (° *) - sinh 4 ‘ pyr} 
Tst (xX) = |] —- - — cosh Bx — l ; 56 
ws hte  GGhbty \te by coth p/ | ' ; ™ 
is the shear stress distribution in the rear spar assuming the wing clamped at x = 0, and 


r(x) = hptAE (= S — cosh sx) [Sthys {PMT} (57 
St, coth BE 
is the correction shear stress distribution. The normal and shear distributions above have been computed on the 
assumption that the applied torque does not introduce normal stresses at the point of application. 
If the torque and moment vectors are perpendicular and parallel, respectively, to the air stream, the following 
transformation may be introduced into all the previous equations developed for torque and moment vectors 
parallel and perpendicular, respectively, to the spars. 


l 0 0 
{PMT} =|0 cosf sinI'| {P’M’T’} 58 
0 -—sinI cos. 
b) Comparison of Experiment and Theory (VIII) ConcLupInG REMARKS 


As an illustration of the validity of the simplified 
theory, Figs. 15 and 16 compare theoretical and experi- 
mental deflections and spar normal stresses for a single 
shear load midway between the spars at the tip. The 
theoretical curves have been corrected for the effects of 


In the discussions above, theories based upon the 
assumption of ribs normal to the spars have been com- 
pared with experimental results from a wing with ribs 
parallel to the air stream. The influence of the rib direc- 
. tion does not seem to be important except for the pos- 
root deformation. sible presence of a root rib normal to the spars which 
closes the triangular section. In such cases, the rib 1s 





heavily loaded, and it should be included in the analysis 
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ribs parallel to the air stream are used and if the most i1- 
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board rib is somewhat removed from the root, it cal 
probably be assumed that the wing outboard of this rib 
behaves like a box with closely spaced rigid ribs normal 
REAR SPAR to the spars. 

As a final illustration of the effects of sweep 0 
hig | stresses and deflections, Figs. 17 and 18 have been pre 
4 8 !2 16 20 24 28 32 36 40 44 48 52 56 60 pared for a four-flange beam with double symmetry. 


STATION ALONG SPAN (INCHES) 3 acu pee ; _ 
5: : ks j , Fig. 17 illustrates the vari: 7 f the ratio ol 
Fic. 9. Spar deflections of sweptback wing resulting from unit es ates the variation with oa + 
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force at Station 2 (obtained by Levy method). rear spar normal and shear stress at the root of the 
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RESULTS OF SWEPTBACK 
swept wing to the corresponding stresses in the straight 
wing for a single shear force applied remote from the 
root and midway between the spars. Fig. 18 illustrates 
similar ratios for tip deformations—namely, tip deflec- 
tion due to a concentrated shear force at the tip mid- 
chord and tip twist due to a concentrated torque at the 
wing tip. The curves of Figs. 17 and 18 are based on 
the simplified theory described in Section (VII), and 
the ratios are developed in Appendix (B). 

It is important to emphasize some simplifying as- 
sumptions employed in the present paper. The problem 
of shear lag in the cover skin and the effects of buckling 
have been neglected. In order to include shear lag 
effects, it is necessary to analyze a cross section that 
more closely approximates the true structure than the 
simplified four-flange box beam. To include shear lag 
greatly complicates the analysis, and a discussion of this 
procedure is presented in reference 7. The effect of 
shear lag is a reduction of structural efficiency resulting 
in increased deflections and increased stresses in the 
corner flanges. Skin and stringer buckling introduce 
nonlinearity into the problem and are usually accounted 
for by approximations. The inclusion of shear lag and 
buckling effects in the swept-wing problem must be a 
subject for further research. 
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Fic. 10. Sweptback wing subjected to concentrated torques 
and vertical loads 














O° Sie ~ Sip 


Fic. 11. Displacement notations. 





ELASTIC RIB 





@, v (x) 
a (x) 
a Bix) we) 
$ (x) a(€) a! 
8(€) 
¢ (€) 


_ 


2. Sweptback wing with closely spaced rigid ribs 
outboard of Station 1. 


Fic. 








716 JOURNAL OF THE 


AERONAUTICAL 





7— ELASTIC: RIB 





T 


dmg 


INNER BOX 





OUTER BOX 


Fic. 13. Structural division of sweptback wing for 
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determination of energy matrix [A,] in Appendix (A 


APPENDIX (A). 


Development of Energy Matrix |K, | 


The energy matrix [A,] of Eq. (37) may be deter- 


mined by the following procedure. The wing structure 


of Fig. 12 is divided into two parts as shown in Fig. 13. 
Here, the inner box is completely elastic with finite rib 
spacing, whereas the outer box is assumed to be elastic 
in warping only. The inner box is subjected to the 
indicated loads, and by applying the Levy mtethod of 


Section (V), the internal strain energy is written as 
U, = (1/2). PMTM,)[C] {| PMTM,} (Al 


The matrix {C] may then be inverted to express the 
inner box strain energy in terms of displacements as 


r ‘ Te ad | j ‘ 
U, = (] 2)[v1a19,¢1|[C] : 0419161 (A2) 


Using the results of reference 10, the strain energy in 
the outer box due to warping ¢g; at the inboard end is 
written as 


U; = (1/2)((h?BAE) /4 coth B/]¢1" (A3) 


The total strain energy in the inner and outer boxes is 
the sum of Eqs. (A2) and (A3), which yields 
U; = 


Ne Quan l 
] 2[v10161¢)] [Ky] V9 ¢1 5 (A4) 
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where 


[Ki] = [C]—-! + [K’] 


0 0 0 O 
0 0 0 O 
(K’]=!|0 00 0 
h*BAE 
0 O 


{ coth £/ 


APPENDIX (B) 


Stress and Deflection Variation with Sweep 


In general, the important design stresses of the swept- 
wing problem are the normal and shear stresses at the 
root of the rear spar. For the test specimen of Section 
(II), these stresses have been plotted versus sweepback 
angle in Fig. 17 as the ratio of swept-wing stress to cor- 
responding straight-wing stress. The curves are based 
on the simplified theory of Section (VII), and are ob- 
tained as follows. 


> 


From Eqs. (53) and (54) we obtain 


osr(O) + of (0) + By 2+ btanl | 
os7(0) x Del 2/btanT)] + Bi 
(Bl 
and from Eqs. (56) and (57), 
tsr(O) + rp(O) + B*y) 2+ btan’T ) 
rsr(0) r $ ([2/(b tan )] + gf 
(B2 


By introducing into Eqs. (B1) and (B2) appropriate 
values for the parameters 8, £, and b, the curves of Fig 
17 may be duplicated for other four-flange, doubly 
symmetrical cross sections. 
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Fic. 15. Mid-chord deflection of sweptback wing resulting Fic. 16. Normal stresses in flanges of sweptback wing re 
sulting from unit shear force at center of Section 2-3 [obtained 


from unit shear force at center of Section 2-3 [obtained by the 
simplified method of Section (VIT)]. 


The deflection curves of Fig. 1S are the result of a 
more complicated derivation. The deflection ratios may 
be obtained as follows. 


From Eqs. (40) and (51), respectively, 


Cor P(E, &) = (2/3AEh?) + (1/2ht,G) — (B3) 
Cy? (g, &) = 
(x Ps h?BAE ow 
"a 4 coth 3/ a Se 
(B4) 


h*BAE 


Kp, K ) — Kp,” 
, ( Mat 4 coth B/ this 


where, from Eq. (47), 
Kp Kp 
Ku, (Kuve v 


K,] = h?BAE ) 


4 coth £/ 

For a concentrated shear load midway between the 
spars at the wing tip (x = &), the ratio of swept-wing tip 
deflection to the corresponding deflection of a straight 
wing 1s 

~ wPrp ¢ "oP(s »\1/f° P(e »¢ - 
Csr (€ & + Cr (& &))/Csr  (&, &) (B5) 
Similarly, from Eqs. (40) and (51), respectively, 


6T /» 


Csr (&, &) = (/ GJ 


l/.. 
plhve + 


Kr( Kure + 


{1 — (n?/ Bl) | (B6) 


h*BAE 


» 
—- K(t)Kp, + K°(&)Kp, 
$ coth =) call : 


b 
h?BAE ” 
4 foe 
4 coth 3/ 
(B7) 
where 
n = [(b/t.) — (ht/te)|/((b/t.) + (h/tw) 
J = 2b7h?/[(b/t.) + (h/ty) | 
K(€) is defined by Eq. (36). 
For a concentrated torque at the tip, the ratio of 
swept-wing tip twist to the corresponding twist of a 


Straight wing is 


by simplified method of Section (VII 
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;. 17. Variation of rear spar root stress with sweepback angle 
for a shear force at the center of Section 2-3. 
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Fic. 18. Variation of tip deflection and twist with sweepback 
angle for a shear force and torque, respectively, at Section 2-3 
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The ratios (B7) and (BS8) are plotted versus sweep 
angle in Fig. 18 for the experimental wing of Section 
(II). 
flange, doubly symmetrical cross sections by introducing 
appropriate parameters into the equations. 


Similar curves may be developed for other four- 











Gyro Pickoff Indications at Arbitrary 
Plane Attitudes 


LEONARD BECKER* 


Lear, Incorporated 


SUMMARY 


The error sensing device most often used in present-day aircraft 
having autopilot control—free or displacement gyros—generally 
gives pickoff indications at arbitrary plane attitudes for small de- 
viations from these attitudes, which are functions not only of the 
departures about the axes that these pickoffs are supposed to de- 
tect, but also the departures about the other two axes. These 
functional relationships are calculated for a common arrangement 
of vertical and directional gyros for an arbitrary attitude of the 
plane and various special attitudes which might occur in maneu- 
vers. In the general case, the resolution of the pickoff signals to 
feed correctly the autopilot servos for stabilization purposes is 
indicated. Gyro interactions which would impair stabilization 
are shown to exist. 


INTRODUCTION 


— USE OF DISPLACEMENT or free gyros as attitude- 
error sensing devices for autopilot control of air- 
These gyros give accurate indica- 
However, 


craft is well known. 
tions for small deviations from level flight. 
in the case where the autopilot is used for stabilization 
while extended maneuvers are taking place, the gyro 
error signals will generally not be accurate indications of 
departures of the plane about its own axes, since these 
inaccuracies arise as a geometrical consequence of the 
inertial properties of the gyro spin axes, and the ar- 
rangements of pickoffs and gimbals. The importance of 
obtaining correct indications of the plane’s perturba- 
tions about its axes, however they may happen to be 
oriented in space, is due to the fact that control mo- 
ments, used to stabilize the plane, are always applied 
about the airplane axes. Furthermore, as stability 
criteria demand that the applied control moments be 
certain functions of the perturbations about the mean 
orientations of the airplane axes, it becomes necessary 
to interpret the apparent perturbations about the space 
axes, represented by the gyros, in terms of the true per- 
turbations, about the plane’s axes, in order that the 
correct control moments may be applied. It is the 
purpose of this paper to calculate the consequences of 
large plane departures from level flight on the pickoff 
readings of a set of commonly arranged vertical and 


directional position gyros. 


PICKOFF AND DIRECTION COSINE RELATIONSHIPS 


Vertical and directional displacement gyros are 
usually free gyros with two gimbal rings which permit 
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the gyro spin axis to remain fixed in direction while the 
gyro case, integral with the airplane, rotates in an arbi 
trary fashion within certain limits. 

Schematically the vertical gyro may be represented 
as shown in Fig. 1. 

The arrangement of the directional gyro is similar to 
that of the vertical unit, except that its spin axis is 
horizontal. 

The gyros are that, 
initially, the orientation of the inertial, airplane, gimbal, 
and spin axes is as shown in Fig. 2. 

The axes may be defined as follows: 


assumed to be mounted so 


X, Y, Z are inertial axes 

Xu» Vay Zq are the airplane symmetry axes through the 
c.g. 

A,,, Asq are spin axes of vertical and directional gyros 

A oor, Aign are Outer and inner gimbal axes of vertical 
gyro 

Aoga, Aiga are Outer and inner gimbal axes of direc- 
tional gyro 


The axis, Aa, can be made to coincide initially with Y 
and Y,, since there will always be a heading at which 
the pickoff signals will be zero. Any deviation from 
this heading will evidence itself in pickoff signals. This 
convention, while greatly simplifying the analysis,” in 
no way decreases its generality. 


ANALYSIS 


If the airplane departs from its initial position, the 
orientation of the gimbal axes relative to the other axes 
will change. From properties and constraints of the 
gyro assemblies, the ensuing gimbal axes’ orientation 
and pickoff angles may readily be ascertained by noting 
the following relations: 


A.—Dnue to structural constraints 
_ (1) Aog: remains coincident with x, 
(2) Aga remains coincident with z, 
(3) Aig always remains | to Ay, and xq 
(4) Aiga always remains | to A,g and 2, 
B.—Dnue to the geometry of the gimbals 


(1) p: = outer pickoff v.g. (roll) = X between yz, 
A ign 

(2) p2 = inner pickoff v.g. (pitch) = { between A,, 
and line which is mutually L to Aoge, Aw 

(3) ps; = outer pickoff d.g. (yaw) = X between 


A igd» Va 
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!) From initial orientation of the axes, it is then 
seen that the initial reading of p; is 0; 
pe is 0; Ps is 90° 


Suppose that the airplane has assumed some large, 
arbitrary attitude. The relation between the airplane 
axes and the inertial axes may then be expressed by the 
transformation matrix: 
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The new orientation of the axes, as defined by the 
transformation matrix [], can then be represented by 
the system shown in Fig. 3. 

From the properties of the gyro assemblies, it can be 
deduced that: (a) X, Y, Aig, are coplanar (due to in- 
ertial rigidity); (b) Xa, Ya, Aigag are coplanar (due to 
structural constraint); (c) Aigg L to plane formed by 
Aq and 2,; (d) Aig, L to plane formed by A;, and x,. 

The pickoff angles /;, p2, p3, which are the only avail- 


Va Va “a 
y ] ; ] able data defining the position of the plane in space, 
: : , can be shown to have the following values in terms of 
[L) = Y ls ls ls the direction cosines: 
A ls ls ly tan Pr => le ly 
—sin po = /; (1) 
where the /’s represent the direction cosines. cot p3 = /2/l5 
Conversely: 
l, = [cos pz (Q®? — cos*p;)*]/Q 
where 
Q = (1/cos ps3) | (cos pe sin p; sin pi — cos pz sin pz)? + cos® p; cos? po} 
l, = (cos pi cos p2)/Q 
lI; = —sin peo 
l, = {(sin po — sin p; cos ps: tan ps3)? + cos? p.(cos* p; — sin® piQ*)}° */Q (2) 
ls = (cos p; cos pe tan p3)/Q 
le = sin pi cos pe 


l; = [cos pi cos pe (tan? p; — Q?)'*/O 


ls = (sin ps — sin pi cos pe tan p3)/Q 
ly _ 


cos p; Cos pe 


In determining these expressions, the relations 
Lh? + 1? + i? = 1 


etc., were used, as well as the orthogonality relations, 


such as 


hls + lls + lle = O 


CALCULATION OF PROPER SIGNALS TO BE INSERTED 


INTO STABILIZATION SYSTEM 


Suppose now that the airplane suffers small yaw, roll, 
and pitch deviations, y, ¢, 0, respectively, from its new 
attitude. In order to stabilize the plane about its new 
attitude, and since the control moments are applied 
about the airplane axes, it is necessary to compute Apy,, 


Aps, Ap; for these small deviations and determine how 
to insert the appropriate error signals into the control 
systems. The steady signals arising from the angles 
Pi, Pe, 3 are assumed to be cancelled out by appropriate 
means. To accomplish this computation, the following 
procedure might be used: 

For a small rotation, y, about z,, the rotation matrix 


is 


Case 
2." Vo’ Zu 
ts [ cos ¥ sin y 0 
[¥] = Me — sin .. cos W 0 
Za 0 o | 1 
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where x,’ and y,’ are the perturbed positions of x,, y, 


The relation between x,’, y.’, 3, and X, Y, Z may be 


found from the matrix multiplication, [y}] & [ZL], the 
result of which is 
ad a 3, 
Y 1 cos y + l sin y L 
; —/, sin Y +/; cos v 
; ls cos +/) sin y 
Wxj=¥| 7er) T l 
—/; sin p +/; cos p 
: /; cos p +/; sin p , 
vA pe ls 
—/, sin y +/, cos 
Similarly, for a rotation, ¢, about x, 
Case II 
\ Va w Sy ” 
\ l 0 0 
>] = y, 0 cos @ sin @ 
Z, O —sin @ cos } 
and by the matrix multiplication, [@] x |Z], 
” ” 
ie Va e 
vy] /, cos @ /; sin @ 
‘ —l; sin @ +/; cos @ 
/; cos @ /5 sin 
@]xX [L]=Y 1 : 4 
—/s sin ® +/: cos @ 
Zz 1 —/, sin @ lg sin @ 
+/; cos @ +1, cos @ 


and for a rotation, 0, about y,, 


Case III 
ere Ya gt 
a cos 6 0 sin 6 
6] = Ya 0 l 0 
Z. —sin 6 0 cos 6 


and by the matrix multiplication, [6] & [ZL], 


AAA yee 


Xa Ya ~a 


+/, cos 9 +/; cos 6 


' —/; sin 6 ‘ +/,; sin 6 
8) ir) =} +/, cos 6 +/, sin 6 
Pa na —/s sin 0 +/, cos 6 
, —l, sin 0 +/, sin 6 
Z sin 


+/3 cos 6 +/, cos 6 


Since only small rotations were made in each case 
from the arbitrary attitude denoted by [L], the new 





rol 
Af 


wh 





oe 


Z may be 


[L , the 


io 
x) 


¢ 


case 


new 


GYRO 
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direction cosines detailing the perturbed attitude may 
be denoted as /; + Ah, ete. 
turbed direction cosines with corresponding elements of 
vy) X [LZ], and also noting that for small enough y, cos 
y ~ l,sin y = y, the following expressions are obtained: 


Comparing these per- 


Case I 
Al = —lIqw Al, = hw Al; = 0 
Als = —leyp Al, = lap Als = 0 
Similarly 
Case II 
Al, = 0 Ak = —-lo Al = lo (3) 
Al, = 0 Als = —Ish Als ls 
Als = 0 Als _ —Ig Aly = led 
Case III 
Al, = —16 Al; = 0 Al; == 1,0 
Als = —/10 Al; = 0 Als 156 
Al, = —1/96 Als = 0 Als = 1,6 





The quantities Ap;, Ap,, Ap; may be found by dif- 


ferentiating Eq. 1. 


Then 
Api = (cos*p, 1g?) (LoAlg = 1gAly) 
Ape = — Al; cos pe r (4) 
Aps = (sin®ps/Is2) (lols — IsAle) | 


In order to express the Api, Aps, Ap; in terms of the 
airplane’s angular deviations 6, ¢, and y from the arbi- 
trary attitude, Eqs. (2) and (3) are used to evaluate (4). 


The results are: 


Case I Api = (—sin p2 cos pi)¥/cos po 
Ap2 = sin pw ; (5) 
Ap; = W J 
CaselI Ap, = —¢ ) 
Ap, = 0 
Aps - t (6) 


—cos p;(sin p2 cos ps — sin Pp; cos p2 sin p3)¢ 


cos p; cos pe 


Case III Ap, = (sin ps sin p1)0/cos peo 
Apo (cos p,) 0 
Aps = : (7) 
sin p;(sin p» cos ps — sin Pp: cos p2 sin ps3) 6 | 


cos P1 cos po 


Since the rotations are small, the results of the three 
rotations may be added, with the resultant signals, 
Api, Apo, Ap; expressible as 


Api = ay + hid aa C10 


Ape = ao + bop + 0:6 (8) 
Aps asw + bsp + c38 J 


where the coefficients are given in Eqs. (5), (6), and (7). 


be | 
bo 
_ 
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Solving for y, @, 6: 


Ap by Ci 


Ap» be Co | 
Aps bs C3 
y= D 
where 
ay b, Cy 
D = Qe be Ce 
az bs ¢ 
ay, Ap C1 
de Ap» Ce 
a3 Abs C3 
2) = 
D 
ay b, Ap; 
de be Ap» 
3 bs A 3 
ae a ) Ap: 7 
D 


Generally, the solutions will be of the form 


vy = d; Ap, + dz Aps + ds Aps | 
o= é Ap + eo Ap» os e3 Aps F (9) 
0 = fi Api + fe Ap2 + fs Aps | 


and an appropriate computer will solve for the y, ¢, and 
6 error signals from the intelligence Ap;, Apo, Aps. 

In the event that the pickoff readings are in the form 
of voltages, which are sinusoidal functions of the angu- 
lar departures, pi, p2, p3, Eq. (9) may be corrected by the 
following procedure: 

Let v; be the reading of the /; pickoff, for instance, so 
that 


2% = Vi sin fp; (V1 being a constant) 


and for small departures, Ap;, there is obtained the ex- 


pression 
Api = Av (Vi cos pi) 


and similarly for Ape, Aps. 

Substituting these expressions in Eqs. (5), (6), and 
(7) and in a similar manner as above solving for y, ¢, 
and 6, the solutions would be in the form 


y = d,;’ Av + dy’ Avo + d;' Avs ) 
d = &'Ar, + e2’Ave + €;’Avs 
Gg = fi’ Avy 4 fo’ Ave + fs’ Avs J 


(9a) 


and again an appropriate computer will solve for the 
y, ¢, and @ error signals from the Av, Ave, Av; readings. 


TURNING ABOUT A VERTICAL AXIS 


Let it be assumed that with the plane in the arbitrary 
attitude denoted by [Z] it is desired to make a change in 
heading by turning about a vertical axis. If the turn is 

, 


through an angle \ so that x,, Ve, 2. goes to Xq’, Va’, Za’, 
the resulting set of direction cosines is 
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a Yo Zq’ tan pi’ = |,/l, 
sin po’ = —/); 

» l A lye l; cos A ee . . . 
x 41, — \ i pee ] . . cot p3;’ = (/2cos\ — |; sin X)/(/; cos X + J, sin A) 

2 sin ) —/; sin —/; sin 

wie = = cot (p3; + Aps) 
Vy ls CoS X l,cosX | —l/gcosX = (cot p; cot Ap; — 1)/(cot p; + cot Ap; 

—/; si ls Si —l; sin X , 

eee | Yee a It will be noted by comparing the above with Eq. (1 

vA ls ls ly that turning about a vertical axis, after the plane as. 


sumes some arbitrary attitude, does not change the /, or 
and from Eq. (1), the new pickoff readings, p,’, p2’, ps’ 2 pickoff readings. The heading angle, H, may be com. 


are given by the relations below: puted from the pickoff readings, thus: 
—/, —=—Coe Pi cos Ps 
tan Hf = = § $ . ° 9 9a ae 9 2 ' l/s 
lh } (cos pz sin p3 sin pi — cos p3 sin pz)? + cos? pi (cos* p2 — cos? pz); “* 
SoME SPECIAL ATTITUDES Se Va 2a 
The result of perturbations from some special atti- x 11 0 0 
tudes into which a plane might be maneuvered will now : : 
be considered. [o>] = } 0 COS >> sin >> 
Zz 0 —sin dp COS dp 


(A) Coordinated Turn (Pitch Attitude Change Negligible) 


If initially the plane is flying straight and level and 
the new attitude is obtained by first rotating the plane and then a turn X\ about a vertical axis, the rotation 


into a bank of angle ¢, the rotation matrix is matrix becomes 
Xa : Va : Za : 
! 
A cos X —cos ¢ sin Xd —sin gd» sin 
[A] X [@] = Y —sin \ COS gp COS X —sin @, cos X 
A 0 —sin dp cos ¢, 


so that the direction cosines become 


l, = cosX ly; = —cos @» sin A l; = —sin @, sin 
ly = — sin X l; = COS dp COS A lj = —sin ¢, cos dX 
ls == () ls = —sin dp ly = COs dd» 
From Eq. (1) Case I 
tan pr = l, ly = —tan od»; pi aa —, (small yaw) P _ ‘ 
—sin po = l; = 0; pe = ( pi a 
tan p; = /;/l, = —cos ¢/tan X Ap. = —sin gd: (10) 


Then, applying Eqs. (5), (6), and (7), one obtains Ap = 





(n 
tic 


bu 


m. 
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Case II Ap; = —tan ¢»-0 (13) 
(small roll) 
in d) Api = —$¢ Thus, as may be seen from Eq. (12), a small pitch 
Ap: = 0 (11) disturbance, 6, will actuate both elevator and rudder. 
Ap; Aps = [(—sin ¢» tan d)/(cos’ g + tan? d) }¢ Rudder motion in turn will give rise to a y disturbance 
Ea. (1 Case III which will again cause the elevator to act (Eq. (10)), 
, (small pitch) giving rise to a possible instability. This kind of inter- 
lane as- “ ; 
he: Ap, = 0 action cannot involve the roll channel because of the 
the /; or : . - on 9 
‘Seal Aps = cos $6 (12) Api = 0 expressions in Eqs. (10) and (12). A roll de- 
Aps = [(—sin ¢» cos ¢»)/(cos? ¢ + tan? dA) ]6 parture will give a rudder signal from Eq. (11), causing a 
y disturbance, but this y disturbance cannot in turn 
WaT ‘ 9 7 ‘ i ‘ ° . 
Eqs. (10), (11), and (12) show that if Ap,, Ap2, and cause an aileron motion. 
Aps are the error signals actuating the control surfaces 
(no computer as given by Eq. (9) is used), no interac- Eq. (13) shows that the effect of @ on Ap; can be ex- 
tions will exist between pitch and roll or yaw and roll, tremely large at high bank angles. Finally, Eq. (12) 
but that pitch to yaw to pitch interactions will exist at shows that the pitch gyro sensitivity, insofar as the 
headings other than \ = 90° or 270°. AtA = O° or 180° pitch control signal is concerned, is decreased by the 
maximum interactions occur where factor, cos dy» 
(B) Climbing Turn 
If the new attitude is obtained by first making a pitch rotation, 6,, the rotation matrix is 
Xa , Va Za , 
X cos 6, 0 sin 4, 
(o,] = Y 0 l 0 
Z —sin 0, 0 cos 6, 
rotation 
and then a bank @,, the matrix of which is 
Xe’ Te” ed 
as” 1 0 0 
[oo] = Va 0 Cos ¢» sin dy 
a 0 -sin dp COs dp 
the resulting matrix is 
%_' %" Ag 
X cos 6, —sin 6, sin gd» sin 6, CoS d» 
[o>] X [0,] = Y 0 COS o> sin ¢» 
Z —sin 6, —cos 6, sin ¢» cos 6, cos d» 
If in addition a turn \ is made about a vertical axis, the matrix becomes 
Pd ed z,/"" 
X cos 6, cos —sin 0, sin @, cos X — cos @¢» sin A sin 6, cos ¢, cos A — sin @» sin A 
[M|] = Y | —cos @, sin dX cos ¢ cos A — sin 4, sin ¢» sin —sin ¢» cos \ — sin 8, cos ¢» sin 
ibs —sin 6, —cos 4 sin @» cos 4, COS d» 
(10) = ~ —— —- 
so that from Eq. (1) and comparison of [Z] and [M] 
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ls 
tan p; = = — tan ®, = — qd» 
ly 
—sin p2 = /; = —sin 6,, pe = 0 
tan ps = /;/l, = (cos d, cos A — sin 8, sin @», sin \)/(—cos 8, sin d) 


and, solving for Ap;, Ape, Ap; from Eqs. (5), (6), and (7), 


Case I 
Ap, = —tan 6, cos dy 
Ap. = — sin dry aa 
Aps _ y 
Case II 
Api => = (01) 
Ap. = 0 “ 
Aps = [(—cos? ps sin 6, — sin ¢» cos 6, sin pz Cos p3)/(COS dp Cos 8y)] > 
Case III 
Ap; = tan 4, sin dy 0 
(16) 


Apz = cos dy 8 


Ap; = [(sin @, sin p3 cos ps + sin @» cos 6, sin? p;)/(cos d» cos 8,) |A 


It may be seen that in this case, if no computer is used, roll and yaw, yaw and pitch interact through the gyros, 
but that roll and pitch still do not interact. The pitch gyro sensitivity is, however, decreased by the factor cos ¢, 


ALTERNATIVE ARRANGEMENT OF VERTICAL GYRO 


In the preceding discussion, the axis of the outer 
gimbal of the vertical gyro was made to coincide 
initially with the X-axis. If, instead, it is made to coin- 
cide with the Y-axis, then a similar discussion will hold. 

The details will not be repeated, but in the case of the 
coordinated turn the result will be 


CaseI Ap, = — tan ® ¥ 
Ap: = 0 (17) 
App = ¥ 
Case lI Ap, = 0 
Ap» = ro) (18) 
Ap; = —[(sin ¢ tan A)/(cos*¢, + tan?Ar)]o 
Case III Ap, = —86 
Ap, = 0 (19) 
Ap; = —[(sin ¢ cos \)/(cos*@ + tan?A)] @ 


Note that in this case Ap; represents the pitching error 
and Ap, the rolling error. Interaction again will exist 
between pitch and yaw if no computer is used, but more 


violently, because of the tan ¢, factor in Eq. (17) as com- 
pared to the sin ¢, factor in Eq. (10). On the other 
hand, none of the gyro sensitivities, insofar as the nor- 
mal control signals are concerned, are lowered by the 
maneuver, ¢, as was the pitch sensitivity in the pre- 
vious arrangement. 


CONCLUSIONS 


Gyro interactions will in general occur between yaw 
pitch and yaw-roll if no computer as given by Eq. (9) is 
The computer necessary, however, would be a 
practical. One 


used. 
complex affair and 
method of reducing the gyro interaction during a co 
ordinated turn would be to cut off the direction gyro 
during such a turn and utilize some other sensing de- 
vice, such as a rate gyro (displacement signals not being 
necessary for the period of the turn), for stabilization 
purposes. It should also be mentioned that in actual 
flight tests such interactions have been noticed but are 
not disturbing for departures from level flight below 


probably not 


about 45°. 
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Effect of Turbulence Level of Incident Air 
Stream on Local Heat Transfer and 
Skin Friction on a Cylinder 


W. H. GIEDT* 
University of California at Berkeley 


ABSTRACT 


At the conclusion of a previous investigation! of local heat- 
transfer and local pressure coefficients around a cylinder, the 
need for simultaneous local skin-friction coefficient measurements 
was indicated. In addition, the results of the skin-friction meas- 
urements of Fage and Falkner? suggested that the influence of the 
turbulence level in the free air stream on the local heat transfer 
should be investigated. Data were therefore taken and are pre 
sented herein for the local heat-transfer coefficient, pressure co- 
efficient, and skin-friction coefficient on a cylinder below and in 
the critical range of Reynolds Number. The effect of artificially 
increased turbulence in the free air stream on these local quanti- 
ties was also determined. Data for these quantities, taken under 
identical flow conditions and on models of the same geometry, 
have not been available previously. The results are partly ex 
plicable in terms of available theory and experience and in part 


present new effects requiring further consideration. 


NOMENCLATURE 


A = area of heat-transfer surface, ft.” 

D = outside diameter of cylinder with nichrome ribbon 
around it, ft. 

= unit heat-transfer coefficient at cylinder surface at 


he = 
angle @, B.t.u./hr.-ft.2 °F. 
= current flowing through nichrome ribbon, amps. 
k = thermal conductivity of air, B.t.u./hr.-ft.? (°F. /ft. ) 
static pressure at any point on the cylinder surface, 
Ib. /ft.? 
t = free-stream static pressure, lb./ft.? 
Ap = pressure reading of surface tube above local static 
pressure, lb. /ft.? 
= rate of heat transfer by convection from cylinder 
surface at any point defined by the angle 6, 
B.t.u./hr. 
R = electrical resistance of nichrome ribbon at 70°F., 
ohms /ft. of length 
= temperature of nichrome ribbon at any point 
around the cylinder, °F. 
= free air-stream temperature, °F. 
U = free air-stream velocity, ft./sec. 
U = velocity at edge of boundary layer, ft./sec. 


average value of turbulent velocity fluctuation be- 
hind grid, f{t./sec. 

= distance around cylinder circumference measured 

from stagnation point, ft. 

= distance from nichrome ribbon measured perpen- 
dicular to cylinder surface, ft. 

slope of velocity distribution at the edge of the 

at the stagnation point, 


boundary layer 


ft./sec. /ft. 
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9 


= absolute viscosity of air, lb.-sec. /ft.” 


M 

0 = angle at center of cylinder measured from the for- 
ward stagnation point, deg. 

p = mass density of air in free stream, lb.-sec.?/ft.‘ 

To = shear stress at wall of cylinder, Ib. /ft.? 

to/pol'o? = (C,/2)y, = intensity of skin friction at the cylinder 
wall based on free-stream velocity 

to/py»U? = (C,/2)y = intensity of skin friction at the cylinder 
wall based on local velocity at the edge of the 
boundary layer 

v = yu/po = kinematic viscosity, ft.?/sec. 

c = 4 - Pe = local pressure coefficient 

(1/2) po Uo? 

Cp = cylinder drag coefficient 

(Vwule = heD/ka = Nusselt Number based on hg 

(Nre) = Dl opp = Reynolds Number based on cylinder 
diameter 

( Np,) = uCp/k, = Prandtl Number for air 

Ae = hae/(poUC,) = Local Stanton Number 


INTRODUCTION 


Ws ON THE INVESTIGATION to be described began 
with measurements of local skin-friction co- 
efficient distributions on a circular cylinder normal to 
an air stream. These data were desired to supplement 
previous heat-transfer and pressure coefficient measure- 
ments on a cylinder! in order to determine if some 
readily recognizable relationship existed between the 
heat-transfer and skin-friction coefficients in the tur- 
bulent boundary-layer region. The results of these 
measurements, which were made with a surface tube 
and in the same air stream as used for taking the pre- 
vious heat-transfer data, showed that, throughout the 
entire range of Reynolds Number (70, 800 to 219,000), 
laminar separation occurred before the 90° point, and 
in no case was there a turbulent boundary layer present 
on the back half of the cylinder as the heat-transfer and 
pressure distributions had indicated. 

This testing was conducted in a 3-ft. open-throat-type 
tunnel (see Fig. 1) with the test section exactly in the 
center of a cylinder, which extended 2 in. on each end 
beyond the edge of the wind stream. Further measure- 
ments were made after installing some plates on the 
cylinder parallel to the air flow about 12 in. on either 
side of the center and placing a '/.-in. mesh wire screen 
across the wind-tunnel nozzle exit. A_ noticeable 
change in the pressure coefficient distribution occurred, 
and, at the higher Reynolds Numbers, a small turbulent 


boundary layer appeared. This indicated that, al- 
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Heat-transfer cylinder as mounted in the wind tunnel 
Photograph taken from inside the tunnel diffusor 
Note cord net over the tunnel nozzle 


Fic. 1 
tor testing. 
section looking upstream 
exit 





Enlarged view (about 2.5 times) of surface tube and 
static pressure taps in the lucite cylinder. 


Fic. 2. 


though the cylinder extended completely through the 
air stream, inflow of some sort had occurred with no 
side plates. A further check on the characteristics of 
the air flow around the cylinder with no side plates was 
obtained by determination of a form drag coefficient 
from the pressure distribution curves. Since in this 
region of Reynolds Number form drag constitutes in 
general 98 per cent or more of the total drag (p. 428, 
reference 3), these values should have compared favor- 
ably with accepted data. They failed to do this, and it 
was therefore considered desirable to repeat the heat- 


AERONAUTICAL 








SCIENCES NOVEMBER, 1931 
transfer and pressure distribution along with skin. 
friction measurements under improved airflow con- 
ditions and variable air-stream turbulence. 

This does not necessarily invalidate the previous data 
but indicates that caution is necessary in its interpreta- 
tion. An examination of the work of Schmidt and 
Wenner,‘ which is the other most recent investigation 

their 
This will 


of this nature, indicated that air stream also 


possessed serious limitations. be discussed 


in more detail later. 


EXPERIMENTAL EQUIPMENT AND PROCEDURE 


The extended test program involved the use of two 
models: a heat-transfer cylinder fitted with electrically 
heated surface ribbons, as shown in Fig. 1, and a cylinder 
fitted with a surface tube, as shown in Fig. 2, for the skin- 


friction and pressure measurements. They were 
mounted, as shown in Fig. 1, in the University of Cali- 
fornia 3-ft. wind tunnel. The environment for the 


tests was, therefore, exactly the same for the heat- 
transfer and skin-friction test cylinders. A more de- 
tailed description of their construction as well as the 
test technique employed is given in the Appendix. 

The primary parameter involved in the test program 
was the variation of the turbulence level of the incident 
air stream. Data were taken first with a 16-mesh damp- 
ing screen placed upstream of the wind-tunnel nozzle in 
order to minimize the turbulent fluctuations occurring 
in the normal tunnel stream. The estimated turbulence 
level, u’/Uo, was negligible, being of the order of less 
than 1 per cent.’ In contrast to this, data were then 
taken with greatly increased tunnel turbulence, accom- 
plished by placing a 1.25- by 1.25-in. rope net (*/j¢-in. 
diameter shot cord) over the nozzle exit. The cylinder 
leading edge was placed about 24 in. behind this net, the 
test arrangement being shown in Fig. 1. The esti- 
mated turbulence level in this case was of the order of 
+ per cent.° Data were also taken at the maximum 
velocity obtainable in the normal tunnel stream insofar 
as this nominally extended the Reynolds Number range 
of the investigation. 

Heat-transfer data and the pressure and skin-friction 
data were taken separately but at as nearly the same 
Reynolds Number as could be realized by the air-stream 
speed control. All Reynolds Numbers are based on free- 
stream velocities calculated from the rise in pressure at 
the cylinder stagnation point referred to the static 
pressure indicated by a pitot tube in the air stream 10 
in. upstream from the cylinder leading edge with the 
stagnation point pressure coefficient assumed to be unity. 


PRESENTATION AND DISCUSSION OF RESULTS 


Local Heat-Transfer, Pressure, and Skin-Friction 
Distributions 
Representative curves of the experimental results 
throughout the test range of Reynolds Number for local 
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TURBULENCE LEVEL OF 
values of heat-transfer, pressure coefficient, and skin- 
friction coefficient distributions are plotted in Figs. 3, 4, 
and 5. Referring to Fig. 3, it is noted that all of the 
heat-transfer curves, with the exception of the one for 
Reynolds Number equal to 99,300, are what might be 
described as characteristic of the critical range—i.e., 
they have two maximums and two minimums (the 
critical range refers to the Reynolds Number region in 
which the rapid drop in the cylinder drag coefficient 
occurs). The data for Reynolds Number equal to 99,- 
300 are definitely associated with the pre-critical range 
with laminar separation as is indicated by the magni- 
tude of the overall drag coefficient (see Fig. 6) and also 
by the local skin friction as shown on Fig. 5. The coin- 
cidence of the location on the cylinder surface of the 
minimum heat-transfer coefficient and the zero value of 
the skin-friction coefficient, which presumably indi- 
cated laminar separation, should be noted. Flow con- 
ditions for Reynolds Number equal to 167,500 appear to 
be of a transitional nature. The skin-friction data do 
not show the presence of a turbulent boundary layer, 
but neither does the skin-friction coefficient decrease to 
zero as rapidly as for the case with Reynolds Number 
equal to 99,300. The curves of Fig. 5 for the other three 
values of Reynolds Number show definite skin-friction 
magnitudes in the region 6 > 90°, and such boundary 
layers as exist in this region can be presumed to be tur- 
bulent in nature. 


The Fig. 4 reveal 
definite inflection points in the 90° region for Reynolds 
Numbers of 90,300, 172,000, and 212,000, as have been 


previously noted by Fage and Falkner and assumed by 


pressure coefficient curves of 


them to indicate the point of transition from laminar to 
turbulent flow in the boundary layer. Form drag co- 
efficients obtained from the local pressure as recom- 
mended by Goldstein (Fig. 161, p. 431, reference 3), 
confirm the flow over the test cylinders as being in the 
critical range. Asa matter of interest, the pressure dis- 
tribution curves given by Schmidt and Wenner‘ were 
integrated also, and the form drag coefficients are also 
shown in Fig. 6. As can be noted, they do not check 
with the accepted curves. The reason for this is that 
their test arrangement did not provide the effect of a 
cylinder in an infinite stream. A curve typical of their 
data at higher values of Reynolds Number is included 
in Fig. 4. This is the type of curve obtained pre- 
viously by the author! and, when considered in the 
light of the experience of this investigation, indicates 
that even at the maximum Reynolds Number 426,000 
obtained by Schmidt and Wenner, laminar separation 
without any turbulent boundary layer probably oc- 


curred. 


Relative Effect of Free-Stream Turbulence 


Comparison of the data taken in the stream behind 
the damping screen to that taken at comparable values 
of free-stream velocities in the stream behind the cord 
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net shows that the general effect of the net is to produce 
results characteristic of velocities higher than those 
indicated by the free-stream velocity. This is in agree- 
ment with results obtained by previous investigators, 
particularly in regard to the marked reduction in form 
drag produced by the increase in free-stream turbu- 
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lence. The average heat-transfer coefficient reveals a For the case of a heated cylinder normal to an air 
similar but not as significant a trend. Fig. 7 shows that stream, where the velocity distribution at the edge of 
the increased turbulence caused the average heat- the boundary layer near the forward stagnation point is 


transfer coefficient to increase between 10 and 20 per 
cent over the low turbulence level values. 

The average heat-transfer coefficient does not indi- 
cate, however, the effect of the increased turbulence on 
the change in the relative amounts of heat transfer from 
the front and back halves of the cylinder. In the low 
turbulence level air stream, the ratio of front- to back- 
half heat transfer is approximately 0.85. In the stream 
behind the net, around 1.1. The 
change is due both to an increase of the heat transfer on 
the front half and to a decrease on the back half of the 


this increases to 


cylinder. 


Heat-Transfer Coefficient at the Stagnation Point 


Fig. 3 reveals that, for comparable values of free- 
stream velocities, the local heat-transfer coefficients at 
the stagnation point as obtained in the stream behind 
the cord net are approximately 25 per cent higher than 
those obtained in the stream behind the screen. 









































a linear function of the distance from that point, so that 
lL’ = Bx, the value of the Nusselt Number at the stag- 
nation determined analytically by 
Squire*® (p. 631) to be: (Nyuw)sr = a3(.Vp-)(8,D?/») 

3, is the slope of the velocity distribution at the edge of 
the boundary layer, D the cylinder diameter, v the 


point has been 


kinematic viscosity, and a;(.Vp,) a known function of 
the Prandtl Number, which is tabulated in reference 3. 
The velocity distributions at the edge of the boundary 


‘layer as determined from the experimental pressure 


data were linear with x. The slopes of these distribu- 
tions were determined graphically and used to deter- 
mine the theoretical value of the Nusselt Number at the 
stagnation point in accordance with the above equa- 
tion. These values as well as the experimentally ob- 
served values of (Vy,,)s7 are listed in Table | and the per 
cent deviation noted. 

The deviations from the theoretical predictions are 
larger for the case of the stream behind the rope 
netting, but a conclusive explanation for the effect 1s 
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TURBULENCE 


TABLE 1 
Comparison of Theoretical and Experimental Values of Vy, at 
the Stagnation Point 


LEVEL OF 





[NOTE: as(.Vp,) = 0.499] 
(Nyu)sr7, Per 
Bi, Theo- N(yvu)s7T, Cent 
Stream (ft. retical Experi- Devia- 
Condition Nre sec. )/ft (Squire) mental tion 
Damping 
screen 99 , 300 565 314 317 9.95 
upstream 167,500 987 407 438 7.6 
Stream be- 95,100 530 307 384 25.1 
hind net 171,500 1,030 414 552 33.3 
Normal 
stream 213 ,000 1,260 460 525 14.0 


lacking at the present time. In this regard, the 
anomalous behavior of the shear stress at the wall as 
indicated by the surface tube is notable. At the stagna- 
tion point and up to approximately S°, in the case of 
the stream behind the net, the reading of the surface 
tube was negative with reference to the local static pres- 
sure. This has been indicated in Fig. 5 by the omission 
of skin-friction magnitudes between 6 equal to 0 and 
about S°. 
phenomenon or neglected to mention it. 


Fage and Falkner® did not observe this 
A possible ex- 
planation of this occurrence consists in the hypothesis 
that the original rise in dynamic pressure in a turbulent 
air stream due to the turbulent velocity fluctuations was 
registered on the static pressure taps but not on the 
surface tube. 
magnitude of the turbulent velocity fluctuations was 


As mentioned previously, the order of 


probably about 4 per cent of the free-stream mean 
velocity, which corresponds to a maximum velocity of 4 
ft. per sec. in these tests. Negative pressures were re- 
corded, however, corresponding to velocities of 10 to 15 
ft. per sec., indicating that some additional or separate 
effect was present. 


The Laminar Boundary-Layer Region 


The increase in local heat transfer, with about the 
same pressure distribution and free-stream velocity, 
produced by the more turbulent free stream in the tests 
made with the net, which has been noted for the stagna- 
tion point region, extended also over the entire reach 
of the laminar boundary-layer flow. Because a com- 
parable effect is not revealed in the skin-friction results 
(Fig. 8), it may be possible that the turbulent flow 
in the region adjacent to the hydrodynamic boundary 
layer was sufficiently nonisotropic to augment the 
molecular conduction of heat in this external region, 
which occurs with fluids of Prandtl Number less than 
unity. Such 
the high heat-transfer coefficients observed. 


account for 
Also, it 
may be possible that the observed phenomenon was due 


an effect would, in part, 


to the action of the turbulent stream near the edge of 
the laminar boundary layer. However, if this were 


true, a slight increase in the skin friction should also 
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result. Referring to Fig. 8, which shows the distribu- 


tion of the skin friction coefficients based on the local 
velocity at the edge of the boundary layer, it can be seen 
that the data for the tests behind the cord net are, in 
general, below or approximately equal to those for the 
stream behind the damping screen. This is felt to be in- 
consistent and could be properly explained only with 
further information about the boundary-layer charac- 
teristics in the turbulent stream. 


The Turbulent Boundary-Layer Region 


A comparison of the local heat transfer and skin 
friction in Figs. 3 and 8 also reveals effects inexplicable 
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in terms of available rationalizations. Fig. 3 shows, for 
the region from about 100° to 130°, substantially 
similar heat-transfer coefficients for the tests at values of 
Reynolds Numbers equal to 172,000 behind the net and 
213,000 in the normal tunnel stream. The skin-friction 
magnitudes, however, are different. This comparison 
is shown in Fig. 9. Thus, for the turbulent flow region, 
a consistent relation between the local Stanton Number 
and the skin-friction coefficient, such as obtains for flat- 
plate flow, does not appear to exist in simple form for a 


cylinder.* 


The Wake Region 


In general, all the heat-transfer coefficientst in the 
wake region have a minimum value at the point of 
separation, whether it be laminar or turbulent, and in- 
crease rapidly to a maximum at the rear stagnation 
point. From the shapes of the curves, it seems possible 
that, due to reverse flow, a boundary layer is formed at 
the rear stagnation point and builds up in the direction 
of the point of separation. This possibility was exam- 
ined by means of surface tube measurements taken with 
the cylinder reversed in the air stream, so that the sur- 
face tube faced the rear stagnation point. Because of 
the small response of the tube in this region, conclusive 
evidence was not obtainable. Some indications tended 
to confirm a backward flow in this region, while others 
indicated fluctuations of the velocity from the forward 
to the backward direction. 
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APPENDIX 


Measurement of Local Heat-Transfer Coefficients 


The heat-transfer cylinder was constructed from 4-in, 
outside diameter and '/s-in. wall thickness Bakelite 
tubing. It was 21 in. overall in length and was mounted 
for testing between side plates 20 in. apart as shown in 
Fig. 1. The model was divided into three sections, of 
which the center 6'/, in. constituted the test section, 
Five separate, 1 by 0.002 in. in cross section, nichrome 
heating ribbons were fitted around this section, the 
center one being the test strip and the two on either side 
acting as guard heating strips. Measurement of current 
input to the ribbons and temperature variation along 
the ribbons by means of extremely fine thermocouples 
located immediately under the center ribbon enabled 
the calculation of the local heat-transfer coefficients, 
hg, due to the relationship 


hy = (g ‘A)/(t — tf) ~ (?R)/(t — b) 


The actual equation used, accounting for conduction 
and radiation losses, is derived and discussed in detail in 
reference 1. The final accuracy of h, is considered to be 
well within 10 per cent except possibly at the maximum 
and minimum points. 


Skin-Friction Measurements with the Surface Tube 


The surface tube located in the wall of the 4-in. lucite 
cylinder is shown in Fig. 2. The cylinder was mounted 
in the 3-ft. wind tunnel in the same manner as the heat- 
transfer cylinder shown in Fig. 1. The surface tube 
itself projected about 0.0035 in. above the cylinder sur- 
face and had an opening of approximately 0.08 by 0.003 
in. <A static pressure tap was located 1'/, in. on either 
side of the surface tube. Experimental measurements 
consisted of determining the difference in pressure be- 
tween it and the average of the static pressure taps as 
the cylinder was rotated about its axis. Several types 
of micro-manometers were tried for this purpose. How- 
ever, it was eventually found that due to the air-stream 
fluctuations, an ordinary inclined manometer with a 
scale having a least count of 0.01 in. of water provides 
sufficient accuracy and was easier to use. Readings 
were estimated to 0.002 in. of water, but an accuracy of 
greater than 0.01 in. cannot be justifiably claimed. 

A considerable amount of time and effort was devoted 
to the calibration of the tube. Space does not permit 
describing this in detail, but, in general, results indi- 
cated that the accuracy claimed in previous work was 
open to question. In addition, the “‘effective distance” 
type of calibration®:* was considered misleading from a 
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On the Stability of Two-Dimensional Laminar 


Jet Flow of Gas 


S. IL. PAIt 
Institute for Fluid Dynamics and Applied Mathematics, University of Maryland 


SUMMARY ‘ 


The stability of two-dimensional laminar jet-type flow of both 
an incompressible fluid and a compressible fluid has been investi- 
gated by the method of small perturbations. Both the sym- 
metrical and the antisymmetrical disturbances are investigated 
separately. 

In the case of incompressible fluid, the solutions of neutral 
stable disturbances in an inviscid fluid are first obtained. Then 
the case of viscous fluid is analyzed. Both the symmetrical and 
the antisymmetrical disturbances are unstable for large Reynolds 
Numbers, but the region of instability of symmetrical disturb- 
ances lies inside that of the antisymmetrical one. The general 
shape of the neutral stable wave number against Reynolds Num- 

curve of the jet-type flow is different from that of the velocity 
profile of ordinary boundary-layer flow or Poiseuille flow. First, 
the effect of viscosity on the neutral stable disturbance in inviscid 
fluid is always stabilizing for jet-type flow. Second, when the 
wave number tends to zero and the Reynolds Number is large, 
the jet-type flow is always unstable. Hence, for the upper branch 
of the neutral stable wave number against Reynolds Number 
curve, wave number increases with Reynolds Number, and there 
is no lower branch of this neutral stable curve for the jet-type 
flow. 

In the case of compressible fluid, the asymptotic behavior of 
the solutions at infinity is discussed. The stability of a two- 
dimensional jet of inviscid fluid is analyzed, and the cases of vis- 
cous fluid are indicated. There is some similarity in the solutions 
of the case of compressible fluid and of the case of incompressible 
fluid. 

So far as laminar stability is concerned, an important difference 
exists between the case of supersonic and subsonic disturbances 
relative to the velocity of the surrounding stream of the jet. The 
amplitude of the subsonic disturbance dies out rapidly with the 
distance from the axis of the jet. The neutral supersonic disturb- 
ances are actually progressive sound waves whose wave length 
and phase velocity are completely arbitrary and that have no 
significance for stability of the jet flow. 

When the Mach Number of the maximum velocity of the jet 
relative to the velocity of the surrounding stream is high enough, 
the jet-type flow tends to be stable with respect to small disturb- 
ances. 


List of Symbols 


The following symbols are used in both Part I and Part II. 


a = velocity in the jet at y + ~ 
“ A nn B., 
ae D,, 


c = 


t functions in asymptotic expression of ¢ at large y 


Cr + tc; complex velocity of propagation of dis- 
turbances 


*F,, F, = determinantal equations for disturbance 

*han = functions in expansion of ¢,° 

1 = / —1 
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functions in expansion of ¢2'” 


*Ron . = 

i = characteristic length 

* Don = functions in expansion of @,''? 

*don+1 = functions in expansion of ¢-‘! 

R = UL/v = Reynolds Number 

t = nondimensional time 

u’ = nondimensional x-wise perturbation velocity 

v’ = nondimensionai y-wise perturbation velocity 

U = characteristic velocity (velocity at the axis of the jet) 

w = nondimensional x-wise main velocity 

x = nondimensional distance along axis of jet 

a = dimensional x-wise distance 

y = nondimensional distance perpendicular to axis of 
jet 

a = wave number 

v = coefficient of kinematical viscosity 

p = nondimensional density 

o = amplitude of disturbance function of y’ in Part I; of 
vin Part II 

*d, +) = functions in asymptotic expression of ¢ at large y 


* The exact expressions for these symbols for Part I and Part 
II are not the same but they are similar so that same symbols 
are used. 

Subscripts 


R 


i = 


real part of the expression 
imaginary part of the expression 


The following symbols are used in Part II only: 


Cp = specific heat at constant pressure 

Pf = amplitude of x-wise perturbation velocity 

M = Mach Number of the velocity at the axis of the jet 
r = amplitude of perturbation density 

p = nondimensional pressure 

z = nondimensional temperature 

Zs = perturbation variable in general 

4 = Qy 

¥ = ratio of specific heats 

6 = amplitude of perturbation temperature 

1, #2 «=60= Ist and 2nd coefficient of viscosity, respectively 
7 = amplitude of perturbation pressure 

a = (1 — {[M%a —c)?]/Ta})'/* 


All the nondimensional quantities are in terms of those values at 
the axis of the jet—i.e., y = 0. 


INTRODUCTION 


OQ OF THE MOST INTERESTING PROBLEMS in fluid dy- 
namics is the problem of hydrodynamic stability. 
It is concerned with the origin of turbulence. Extensive 
theoretical investigations of this problem were made by 
Tietjens,' Heisenberg,® and Tollmien* * in the 1920's. 
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In recent years much investigation,* both theoretical 
and experimental,® has been carried out for various 
types of laminar flow. However, most of these investi- 
gations deal with those classes of velocity profile that 
involves at least one solid boundary—e.g., boundary- 
iayer flow, flow through a channel, etc., but little has 
been done for the class of velocity profile of jet type 
i.e., the range of the velocity perpendicular to the 
main flow extending to both plus and minus infinity. 

The stability of velocity distribution of jet-type flow 
of incompressible fluid was first treated by Hollingdale® 
and Savic.’ Their results are limited to inviscid fluids 
and their assumption that symmetrical disturbances do 
not produce instability is not correct. Lessen* and 
Chiarulli® studied the stability of the half-jet-type flow. 
Chiarulli has developed a simplified method of calcula- 
tion which is very useful in the study of stability of jet- 
type flow. But he did not carry out the calculations. 
Little has been done about two-dimensional laminar 
jet-type flow, including viscosity effect. 

In the first part of the present paper, the stability of 
two-dimensional laminar jet-type flow of an incompres- 
sible fluid has been investigated by Heisenberg-Lin" 
method with the help of Chiarulli’s simplified method of 
calculation.’ Both the symmetrical and the antisym- 
metrical disturbances are investigated separately. The 
solutions of neutral stable disturbance in an inviscid 
fluid are first obtained. Then the cases of viscous fluid 
are analyzed. The eigen value equations are obtained in 
closed form and the calculation of the eigen values is re- 
duced to straightforward evaluation of certain definite 
integrals. The neutral stable curve in the wave number 
a and Reynolds Number R plane can be thus calcu- 
The general shapes of the neutral curve are dis- 
cussed. It is shown that the general shape of the neutral 
curve of the jet-type flow is quite different from that of 
the velocity profile of ordinary boundary-layer flow or 


lated. 


Poiseuille flow. 

The velocity profile of the main flow used in the calcu- 
lation is that of the two-dimensional jet issuing from a 
small slit into fluid at rest. We assume that the fluid is 
incompressible and that the Reynolds Number is large; 
hence the velocity in the jet must be small and the dis- 
tance from the nozzle must be large. As a result, the 
velocity profile of jet from small slit is justified to be 
However, the case for velocity profile of finite 
The essential char- 


used. 
opening of the nozzle is discussed. 
acteristics of the neutral curve will be the same for the 
finite opening case as that for the small slit but the 
numerical values of the neutral aR curve will be dif- 
ferent in the two cases. 

Because of the increasing importance of high speed 
flow from jet engine, it seems natural that the above in- 
vestigation should be extended to cover the case of 


* For additional references see Lin’s paper (reference 10) and 
also Dryden’s paper, Recent Advances in the Mechanics of Bound- 
ary Layer Flow, Advances in Applied Mechanics, Vol. I, 1948. 
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compressible fluid taking into account the effects of 
compressibility and heat transfer. 

The interest in the problem of compressible fluid js 
further enhanced by the fact that the laminar flow of 
very high Mach Number tends to be completely stable 
with respect to small disturbances. Lees'! has shown 
that for Mach Number greater than three, the laminar 
boundary-layer flow for thermal equilibrium is com- 
pletely stable at all Reynolds Numbers under free-flight 
conditions if the free-stream velocity is uniform. 
Landau'*? and Hiroshi Hatanaka!*® both showed that if 
the relative Mach Number of a surface of discontinuity 
in an inviscid compressible fluid is large enough the 
vortex sheet will be stable with respect to small dis- 
turbances. Hence, it will be interesting to determine 
whether the jet-type flow will be stable at high Mach 
Number for all Reynolds Numbers. If so, what would 
be the approximate value of this critical Mach Number? 

As in the incompressible fluid case, only small two- 
dimensional wavy disturbances will be considered for 
the compressible fluid case. Unfortunately, there is no 
rigorous proof in the compressible fluid case that these 
disturbances are more unstable than the three-dimen- 
sional ones. But the results of the incompressible fluid" 
together with some physical considerations seem to 
justify such a treatment that naturally simplifies the 
mathematical analysis to a large extent. 

One essential difference in the case of compressible 
fluid from that of incompressible fluid is the presence of 
an appreciable interchange of mechanical and _ heat 
energies. Another is the fact that the flow velocity is of 
the same order of magnitude as the velocity of sound. 

The method of analysis used for compressible fluid is 
similar to that used in the incompressible fluid 
case. Special reference should be made to the work of 
Lees and Lin who studied the stability equations of 
two-dimensional parallel flow of compressible fluid aud 
their solutions in great details. Most of their results 
will be used in this paper. Here our emphasis is placed 
on the case of laminar jet flow. 

A numerical integration method is devised to com- 
pute the velocity and the temperature distributions in 


very 


the two-dimensional jet-type flow without any specific 
assumptions about the value of Prandtl Number and the 
law of variation of viscosity with temperature. 

The method described in this paper may be also used 
to study the stability of two-dimensional wake. 


Part (I) Jet FLow or INCOMPRESSIBLE FLUID 


Orr-Sommerfeld Equation and Its Solutions 


The fundamental equation of stability of two-dimen- 
sional parallel flow and its solutions have been investi- 
gated extensively by Lin.” We are going to use this 
method to study the stability of the velocity profile of 
jet-type flow. First we will summarize some of Lin's 
results which are useful for present analysis as follows: 
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Consider the stability of the velocity distribution 
w = w(y), by superimposing upon the steady flow a 
small, hydrodynamically possible disturbance obtained 
from a disturbance stream function. 


v’ = oye" -” (1.1) 


where wave number a is a real and positive number and 
velocity of propagation of disturbance c is in general 
complex. 

If we substitute the resulting disturbed flow into the 
equations of motion and linearize them, we obtain the 
stability differential equation that is known as Orr- 
Sommerfeld equation as follows: 


(w —c) (¢” — a’) — wd = —(i/aR) X 
(¢” — 2a*o” + atd) (1.2) 


In this equation all velocities are referred to a charac- 
teristic velocity L’, all lengths to a characteristic length 
L; Uand L may be functions of x. The Reynolds Num- 
ber R is equal to L/L/v where v is the coefficient of 
kinematic viscosity. Prime represents the derivative 
with respect to y. 

The stability differential Eq. ( 
system of four solutions that are analytic functions of 
(wherever w(y) is analytic) and are entire functions of 
Two of the solutions, called by Lin 


.2) has a fundamental 


a, c, and aR. 
“viscous solutions,’’ should be rejected in the study of 
flow of jet type because of the bound of the solution. 
The other two solutions, known as “inviscid solutions,’’ 
can be obtained by developing the solution in powers of 


(aR)~!. 


o(v) = 6 (y) + (aR)~'o@ (y) + 


(aR)~*o (vy) +... (1.3) 


Following Heisenberg’s procedure, we develop ¢“ (y) 
in power of a®. In terms of Lin’s notation, the resulting 


solutions for ¢ (y) are 


¢o: = (w — c)[ho(y) + a®ho(y) + athy(y) +... ]C1.4) 


do = (w — c)[Rkily) + a*ks(y) + atks(y) +... J(1.5) 


where /ip = 0, he, +2’s and k»,+ 1's are functions of y, w, and 
¢ which are given in reference 10. 

For velocity distribution of jet-type flow, we must in- 
clude the ¢) solutions in the eigen value problem in 
order to take the effect of viscosity into account.'® De- 
veloping the solutions in powers of a’, we obtain 


o: = (w — c)[ poly) + a®pol(y) + atps(y) + ...] (1.6) 
g:' = (w — c)[qily) + a?qs(y) + atgs(y) +...] (1.7) 


where p»,’s and qz,+1’s are functions of y, w, and c. In 
reference 9, Chiarulli showed how to calculate these p», 
and q2,+, but some of the formulas given in his report 
seem to be incorrect. 

Eqs. (1.6) and (1.7) are not unique in that any linear 
combination of ¢; and ¢2 can be added to them. But 
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as shown by Chiarulli,’ it is sufficient to use these ex- 
pressions in the eigen value problem. 

The solution of the stability differential equation for 
the velocity distribution of jet-type flow is 


v1) = (19) + Coe (1.8) 


where ¢; and c, are arbitrary constants. 


Asymptotic Behavior of the Solution at Infinity 


In the boundary value problem, we have to know the 
value of ¢, for large values of y. By the method of 
Chiarulli,? we obtain the value of the ¢,” for y ap- 


proaches to — © and w approaches to a as follows: 
dD? = 4%" + BV (1.9) 


where 7 = 0 and 1 and k = 1 and 2. 


c-¢ 
@,.°° = — “ [Ay + aA; + aA» + @'®Az +... ] 
(1.10) 
c-—a By 
4.9% = — = + +B, +aB,4+... (1.11) 
~ Q@ 
has 
$.° = — | [Fy + ak, + a®k; + | (0.12) 
c-—a 
o, = 2 [Dy + aD, + a*D.+...!} (1.13) 
where Ay = 1, A,, B,, /,, and D, are definite integrals 


from minus infinity to some reference point yo of the 
function of w, a, and c. Since the velocity distribution 
w(y) which we are dealing with approaches to a con- 
stant value at infinity in an exponential manner," all 
the infinite integrals converge. 

Similar expressions for @’s can be obtained for y> + ~ 


and w—da. 


Boundary Value Problems 


We are going to study the boundary value problem 
for the case of two-dimensional jet flow. The velocity 
distribution w(y) of the main flow is symmetrical with 
respect to y. Any function ¢(y) satisfying Orr-Sommer- 
feld equation may be separated into two parts, re- 
spectively odd or even in y, each of which separately 
must be a solution of the differential equation. The 
even solution corresponds to antisymmetric oscillation 
(an alternate vortex sheet), the odd solution, symmetric 
We are going to consider these two types 
The symmetry condition 


oscillation. 
of disturbances separately. 
allows us to consider only negative value of y. 

(1) Symmetrical Disturbance.—lf $(y) is an odd func- 
tion of y—i.e., symmetrical oscillation—the boundary 


conditions are 


¢(—~) = ¢(0) = 0 (1.14) 


The determinantal equation corresponding to these con- 


ditions is 
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, ae 65) sis ; 
F\(a,c,aR) = + op eee = 0 (1.15) 

¢1(0) o2(0) 


(it) Antisymmetrical Disturbances.—-If $(y) is an even 
function of y, antisymmetrical oscillation, the boundary 
conditions are 


o(—~) = ¢/(0) = 0 (1.16) 


The determinantal equation corresponding to these 
conditions is 


2( 
¢2'(0) 


@ ) 


gi(— ©) 


= ) Pf 
$1'(0) O (1.17) 


F(a, c, aR) = 


Velocity Distribution of the Main Flow 


(t) Small Slit-~—The problem of the jet going out 
from a long narrow slit in a wall and mixing with the 
surrounding fluid at rest is one of the rare cases where 
the differential equation of boundary layer may be 
integrated exactly. The calculations were performed 
by Schlichting" and Bickley." 


Let 
lL’ = characteristic velocity = 0.4543(Mo?/po’vX)'* 
L = characteristic length = (\M/o/pov?X?)~'' 
where 
M, = rate of momentum across a section of the jet 


* 


= 2U* ool / wdy 


po is the density of the fluid. 
The nondimensional velocity distribution in the jet 


from a small slit will be 


w = sech? by (1.18) 
where 
9\' 
hb = 
3\16 
(17) Finite Opening.—In case the jet discharges 


from a nozzle of finite opening into the surrounding fluid 
at rest or in uniform motion, the velocity distribution 
can be numerically integrated by the method given by 
the author in reference 17. It is noticed that when the 
distance is small from the exit of the nozzle, the velocity 
distributions at various sections are not similar. Hence 
in the study of stability of flow near the exit, it is neces- 
sary to find the whole series of distribution first and to 
make the stability calculations for each velocity dis- 
However, the procedures of calculations for 
The problem 


tribution. 
each velocity distribution are identical. 
should not introduce any complications. 
paper the stability calculations are based on the velocity 
distribution of small slit in order to bring out some 
essential features of stability of velocity distribution of 
The calculation for velocity profiles of finite 


In the present 


jet type. 
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opening will be carried out in the future in order to de. 
termine the influence of change of velocity profile from 


uniform one to that at a large distance. 


Neutral Stability in an Inviscid Fluid 


The study of the stability problem in an inviscid fluid 
may be regarded as the limiting case where the Reynolds 
Number becomes indefinitely large. The stability dif- 
ferential Eq. (1.2) reduces to 


(w — c) (6 — a®d) — w"d = 0 (1.19) 


This equation has a singularity at w = c. It is well 
known that in an inviscid fluid, all velocity profiles with 
a point of inflection are unstable and tend to exhibit os 
cillation.‘ Both increased and neutral oscillation (c en 
tirely real) exist and for the neutral ones the wave ve 
locity is equal to the stream velocity at the point of in- 


flection of the velocity profile. From Eq. (1.18) we have 


w"(y,) = w,” = 0, wy.) = Ww = 2/3 (1.20 


Putting ¢ = w, and Eq. (1.18) into Eq. (1.19) we have 


(1.21 


o” — a’d + 6b*o sech? by = 0 
For symmetrical oscillation, the solution of Eq. (1.21 
is 


(1.22 


@ = k sech by tanh by 
where & is an arbitrary constant and the eigenvalue oa 
will be 
= b: (1.23 


a Ms 0.27516 


For antisymmetrical oscillation, the solution of Eq 
(1.21) is 

¢ = k sech® by (1.24 
with the eigenvalue 
= 2b 


ie Tas = (0.55082 


Eq. (1.24) was first found by Savie in 1941.7 


Stability in a Viscous Fluid 


Now we are going to study the stability of jet-type 


flow in a viscous fluid—i.e., for finite Reynolds Number 
Here we should deal with the boundary value problems 
15) and (1.17). We shall find the 
We take a as positive real 


formulated in Eqs. (1 
eigenvalues of a, c, and R. 
number; c may be complex. R must be real and positive 
in order to have some physical significance. As shown 
before, we treat the symmetrical and the antisym 
metrical oscillations separately. Furthermore, a 1s 
assumed to be small and a@R is large so that all the terms 
of the order of 0(1/a°?R*) may be neglected. 

(i) Symmetrical Disturbance.—After neglecting the 
higher order terms, the determinantal Eq. (1.15) be 


comes 
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(0 + 


F(a,c,aR) = , 


1 : w”(O) a? 
P, . = .% e + 
aR (1 — c)* l _— 


Je. | (1.26) 
c 


If we put c{to be real value, the eigenvalue a obtained 
from Eq. (1.26) in order to get positive real value of R 
gives the corresponding values for neutral aR curve. 
In the study of hydrodynamical stability, it is usual to 
plot the neutral aR curve to show the region of stability. 
It is interesting to examine the general shape of the 
neutral aR curves for jet-type flow and determine 
whether or not they are similar to those of ordinary 
boundary-layer flow or Poiseuille flow. For ordinary 
boundary-layer flow or Poiseuille flow, there are two 
branches of the neutral curve. For the lower branch, 
one has a—Oand R—~ ~. Let us substitute the rela- 
tions a > 0 and aR is finite but large in Eq. (1.26) we 
have 
ony” 
inp. (1.27) 


(1 — c)? 

Eq. (1.27) represents the characteristic of the flow as 
a>0,and R—~ o. Since w”(0) is a real and negative 
value, for any real value of a and c, R will be purely 
Therefore there will be no lower branch of 


imaginary. 
neutral aR curve in the present case.* For real value of 
R, we must have complex c. Put 
C = Cr + IC; (1.28) 
into Eq. (1.27), and we have for real value of R 
(1 — cr)? — c? = 0 (1.29) 
aR = —[w"(0)]/[2¢,(1 — cr)] (1.30) 


If cp < 1, ¢ > O, and if cp > 1, c; < 0. Hence for first 
approximation, the flow is unstable for symmetrical dis- 
turbance at very small a when the velocity of propaga- 
tion of the disturbance is smaller than the maximum 
velocity in the jet but stable when the velocity of propa- 
gation of disturbance is larger than the maximum veloc- 
ity in the jet. Ordinarily ce < 1. We may conclude 
that there is no lower branch of neutral aR curve for 
symmetrical disturbance of jet-type flow and that, as a 
is sufficiently small and R is large, the jet-type flow is 
always unstable with respect to symmetrical disturb- 
ance. 

Next we investigate the effect of viscosity on the 
neutral stable disturbance in inviscid fluid. If there is 


* Prof. C. C. Lin pointed out to the author that according to 
Synge’s criteria,”’ the lower branch of the neutral stable curve 
might still exist for such a low value of aR that the asymptotic 
no longer holds true. This point should be 
However, if it were true, then 


solution, Eq. (1.3), 
clarified by further investigation. 
even for moderate Reynolds Number, @ of the lower branch 
neutral curve is so small that it would practically coincide with 
the R-axis. Hence for first approximation, we may still consider 


that such a lower branch does not exist. 
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an inviscid neutral disturbance with c = c, and a = a, 


we have 
F,(as, CG, R= o~) = 0 


When a = a, ¢ = ¢, + Ac, differing but slightly from c, 
for finite but large R. We have from Eqs. (1.26) and 
(1.31) after neglecting the higher order terms 


(Cx, @ ) 


P. (1) 
a, R/O’, id a s 


Ac = Ace, + 1Ac; (1.32) 


In Eq. (1.82), for first approximation, we have 


Cc = 


l ae _ Biz’) + | (1.398 
3) 


R (Bo’ = on Byipr')? + B,"? a2 
hayeet kw ——— eed (1.34) 
R (Bo —- Bir x + By; 2a,” c 


From numerical calculations of mean velocity of jet 
from small slit, it is found that Ac; is a negative quan- 
tity. Hence the viscosity effect is stabilizing. The 
change in cz is rather small. It shows the difference of 
the upper branch of the neutral a, curve of the jet-type 
flow from that of the ordinary boundary-layer flow. An 
inviscid neutral symmetrical disturbance with non- 
vanishing wave number and phase velocity will be 
stabilizing by the effect of viscosity in jet-type flow. 

For the actual shape of the neutral aR curve we 
should calculate the eigenvalue a from Eq. (1.26) for 
various real values of c. First we should evaluate the 
integrals B,, Dy, etc. Since there is a singularity at 
w = c, usually it is necessary to perform the integration 
in complex plane. Chiarulli® suggested a method in 
which the numerical integration along the complex 
contour may be avoided. We extend his method to 
perform these integrations. 

(it) Antisymmetrical Disturbance. 
the higher order terms, the determinantal Eq. (1.17) be- 


After neglecting 


comes 
F(a, ¢ ak) = @,— + 


| ~(1 — ocw"(0) + ae - 

?; yp—} . Pp,“ = (0 (1.35) 
aR (1 — c)’ 

First we examine the lower branch of the neutral aR 
Put 


curve and see whether or not it exists. a> 0, 


aR is finite in Eq. (1.35), and we have 
aR = i }[w"(0)]/(1 — c)*} — & (1.36) 
For real values of c from 0 to 1, it is found that R is al- 
ways complex if a is real. Hence there is no lower 
branch of this neutral aR curve.* 
If we put Eq. (1.28) into Eq. (1.36), we have 


Eo(Cr + ic = —c Kor, + 15; 
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aR = » | a? I ; 
pli(w — c)p] = — + 7 A(Z) (9: 
2w”(O) (1 — cp) ? vill? oR ul (2.9 


= + Fan, | = ¢K 


[(1 — cr)? — ©;7]? + 4c? (1 — cr)? 


For ce < 1, it was found that K is always positive. 
Hence if aR is real and positive, c; must be positive. As 
a result the flow is unstable with respect to antisym- 
metrical disturbance for a ~ 0 and R is large. 

Next we investigate the effect of viscosity on the 


neutral disturbance in inviscid fluid. If there is an in- 


viscid neutral disturbance with c = c, and a = a, we 
have 

FilG,, a, R= oo) =O (1.38) 
When a = a, and c = c, + Ac differing but slightly 


from c, for finite but large R, we have from Eqs. (1.35) 
and (1.38) after nelgecting the higher order terms 


@,—“(C,, as) 


Ac = —- = Ace + tAc; (1.39) 
aR[OP— /OC]a, ; 
In Eq. (1.39) for first approximation, we have 
I Kui © OF : 

Ac; = — (1.40) 

aR a(0A1/Oc),. 

1 For ir 
a R asliin (1.41) 


aR a(O0A,/Oc),, 


For mean velocity distribution of jet from small slit it is 
found that Ac; is a negative quantity. Hence the effect 
of viscosity is stabilizing. An inviscid antisymmetrical 


disturbance with nonvanishing wave number and 
phase velocity will be also stabilizing by the effect of 


viscosity in jet-type flow. 


Part (II) Jet FLow or COMPRESSIBLE FLUID 


Lees-Lin Equations and Their Solutions 


The general equations governing one normal mode of 
the small velocity and temperature disturbances of two- 
dimensional laminar flows of gas and their solutions are 
studied in details by Lees and Lin.” We are going to 
use their method to study the stability of the velocity 
profile of jet-type flow. First we summarize some of 
their results which are useful for the present analysis 
and then derive some equations from the general equa- 
tions which are important for the problem of jet-type 
flow. In the present analysis, the effect of gravity is 
neglected, hence all the terms containing the Froude 
number in Lees-Lin equations are omitted. 

The final dimensionless and linearized equations for 
the amplitudes of small disturbances in a parallel or 
nearly parallel main flow are as follows: 

ir l 


i(w—c)f+w'o] = —- + 


(w—-ort+ pd +i/)+o0¢=0 (22 

pli(w — «6 + T’o] = — (y — 1)X 
pl (o’ + if) + (/aR)F i (Z,;) (24 
w/p = (r/p) + (0/T) (2.5 


where the nondimensional disturbance variables are de- 
fined as follows: 


u’ = x-wise perturbation velocity component 
nr” ee 
v’ = y-wise perturbation velocity component = 
aie’ (29 
perturbation density = r(y)e*" 8 
perturbation pressure = (y)e’“" (2.9 
perturbation temperature = @(y)e’**~ (2.10 


F,(Z;), Fry(Z;) and F,,;;(Z;) are linear combinations of 
i.e., /, d, 7, r, 8, and their 
first derivatives with respect to y. 


the perturbation variables Z; 


The system of five equations, Eqs. (2.1) to (2.5), is 
actually equivalent tosix homogeneous linear differential 
equations of the first order and there are six linearly in- 
dependent solutions. The system of equations can be 
solved by the method of successive approximation by 
means of asymptotic series. Two asymptotic series are 
obtained by 
(aR) ~'—1.e., 


Z, = 2; 


expanding the solution in powers of 


+ (aR)2Z,;0 + (aR) -?Z;° = 2 (2.11 


The system of Eq. (2.1) to (2.5) may be written as 
a differential equation for each of the independent 
For example, the differential equation of the 
i.e., @ will be 


variables. 
amplitude of velocity component v’ 


d[(w—c)d’ — “4 a(w—c)do | 
dy T — Al*(w — c)? ti aR 


where G(Z;) is a linear combination of Z; and their first 


G(Z 


(2.2 


derivatives. 
Substituting Eq. (2.11) into Eq. (2.12) we have 


d[(w— c)d®’ — w’d a*(w — cid 
= (2.13 
dy. T — M*(~w — cc)? T 


dl (w— c)d” — w'd” a*(w — cd” 
dyL T — Mw — c)? T : 


GZ") + a?®Gil(Z,"—-) + atGo(Z,"— (xn = 1 


where G; is independent of a and R, and is a linear com- 
bination of Z;"~! 
As soon as ¢ is obtained, the other perturbation vari- 


and its first derivatives. 


ables can be calculated from the following equations: 
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Mw — c)w’o — Te’ 
=~! T — M*w — c)? bi 
i [ Mew — OF; Fi; | (943 
aR p 7 Hw — ~ sie 
r= 1[ p(o’ + if) + p’o|/(w — c) (2.16) 
r= iy? }oli(w — c)f + w’e| — (1 aR)F,} (2.17) 
6= 7|(x p) — (r/p)] (2.18) 


Similar to the incompressible case, Eq. (2.12) can be 


solved by developing ¢ 


in power of a?—1.e., 


6) = (w — c)[ho(y) + ahe(y) + atha(y) +...) (2.19) 
Qo» y= (ty — ¢ ki(y) +- a*k3(y) — @ 'ks(y) + rose (2.20 
where iy = 1, fo, and k»,+; are functions of w, c, 17°, 
and 7 which are given in reference 15. 

For the @ ' solutions we have 
6°) = (w — c)l poly) + a*po(y) + atpa(y) +...] (2.21 
d.') = (w — c)laily) + a7qs(y) + atgs(y) +...) (2.22) 


where 
which are given in Appendix A. 


p», and gn+1 are functions of w, c, M*, and T 


Four other asymptotic solutions are obtained by put- 


ting 


lo 
tw 


ii = I (aR) Sg dy| \< 


exp 
where g is independent of aR. For a similar reason as 
in the fluid 
should be rejected in the study of stability of jet type 


incompressible case, these solutions 


flow. 
Finally we have 


6 

8 

t 
tw 
to 


Asymptotic Behavior of the Solution at Infinity 


We are going to investigate the behavior of $”’ for 
large values of vy. By similar procedure as that for in- 


compressible case, we find the following expressions of 


¢."" and ¢, ' for large negative y. 
Atyv—> ©, we assume that 
wy=—o)=a, Tiy=—0)=T 


Then for large y 


where 


Ma — 0)2/T,]} "vy (2.26) 


s=Qy = ;,1- 
may be written in the same 
(1.10) to (1.13) except that 
D,, are different. 


and the expressions of ®,_ 
form as those given in Eqs. 
the expressions of A,, B,, /,, and 
These expressions for compressible fluid are given in 
Appendix B. 


LAMINAR JET 


-— 
ay 


FLOW OF GAS 


Similar expressions for @”’ can be obtained for 


y>+o andw>a,T>T,,. 


Subsonic, Sonic, and Supersonic Disturbances 


It is convenient in the study of the stability of flow 
of compressible fluid to classify the disturbance as ‘‘sub- 
sonic,’ “‘sonic,’’ and “‘supersonic’’ relative to certain 
local velocity u* when the x-component of the phase 
velocity of the disturbance relative to the local velocity 
is less than, equal to, or greater than, the mean speed of 
sound at that point. The conditions for this classifica- 


tion are 


< supersonic 
i M?(w — Cr)? sonic (2:27) 
> ] subsonic 


If w = a, T = T,, the subsonic disturbance relative to 


a must satisfy the following inequality 


¢ 
— = 


a — (T,/M?)"? < ce <a t+ (T,/M?)"? (2.28) 


For the case of no energy transfer—.e., the relation be- 


tween the temperature and velocity is given below 
Ta + [Cy — 1)/2]M*a? = 1 + [(y — 1)/2]M? (2.29) 


Eq. (2.28) becomes 
a- = = (1 — a’) 
(2.30) 


= (1 — a’) | 


The limiting curves of Eq. 


_ 


a may vary from 0 to 1. 
(2.30) are plotted in Fig. 1, which shows the region of the 
classification of the 
It is interesting to see that no disturbance 


disturbances relative to local 
velocity. 
may be considered as subsonic relative to the entire flow 
field if the maximum Mach Number of the flow field is 
greater than 2.5 for the case of no energy transfer. 

0? will be greater than, equal to, or less than zero ac- 
cording to whether the disturbance is subsonic, sonic, or 
supersonic with respect to the velocity a at infinity. 
Since the value of 2°? determines the nature of the solu- 
tions of the disturbance differential equation we must 
study its value in detail. For subsonic disturbance 
relative to the velocity at infinity, the amplitude of the 
disturbance will die out rapidly with distance from the 
axis of the jet and discrete eigenvalues may be found. 
The neutral subsonic disturbance is an ‘‘eigen-oscilla- 
tion’ in the jet. These disturbances are important in 
the study of stability of jet flow of compressible fluid 
and determine the limits of stability of the jet. For 
supersonic disturbance relative to velocity at infinity, 
Q is imaginary, the eigenvalues are continuous and not 
discrete. this type the 
length and phase velocity are completely arbitrary be- 


For disturbances of wave 
cause the boundary condition at infinity is missing here. 


These disturbances behave as progressive sound waves 
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- 16 The determinantal equation is then 
\ SUPERGONIC | o(— »(— 
Ph F,(a, c, R, M?) = oi(— ©) do(— ©) = 0 (2.34) 
¢:(0) $2(0) 
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Fic. 1. Classification of disturbance for jet without heat 
transfer. 


in the jet and they have no significance for the stability 
of the jet. 

For jet with the surrounding stream at rest, i.e., a = 
0, the line a = O in Fig. 1 separates the subsonic dis- 
turbance and the supersonic disturbance for the case 
without energy transfer. Any disturbance with 


: fey /ape /Q]\ V2 

Ce > {(1/M?) + [(y — 1)/2]} 
should be regarded as supersonic disturbance. 
In this paper we are concerned with only the sub- 
sonic disturbance relative to the velocity at infinity. 
For the case with surrounding stream at rest, the 


physical condition for subsonic disturbance for the case 
of no energy transfer is 


Cr < {(1/M?) + [(y — 1)/2]} (2.32) 
In general, only those disturbances satisfying the con- 
dition of Eq. (2.28) are considered in the present paper. 


a oa 


Boundary Value Problems 


Similar to the incompressible fluid case, we treat the 


symmetrical and the antisymmetrical disturbances 
separately. 

(7) Symmetrical Disturbance.—Here the boundary 
conditions are 


¢(—@) = g(0) = 0 (2.33) 


(17) Antisymmetrical Disturbance.—Here the bound- 


ary conditions are 
o(—-) = $'(0) = 0 (2.35) 
The determinantal equation is then 


: 2» _ | (—-)gde(—-) _ — 
F(a, c, R, M?) = $:'(0) (0) = 0 (2.36) 


Velocity and Temperature Distribution of the Main Flow 


The velocity and the temperature distributions of the 
main flow can be calculated either by the method given 
by the author in reference 17 with the assumption of 
Prandtl’s Number to be unity or by the following step- 
wise numerical integration method. 

In terms of stream function ~ and x-wise distance, the 
equations of motion without pressure gradient of the 
two-dimensional jet flow and the equations of con- 
tinuity give the following nondimensional equation 


Ow ie o ( ~ we) (2.37) 
Oxo Wye ” dv iis 


The nondimensional energy equation will be 


OcpTy 0 — oT Ow \* oF 
Oxy ak vais a) om Kat sa 
and the initial conditions are 
x = 0, To = Tolvo), wo = wolWo) (2.39) 
The equation of state 
po = 7) (2.40) 
and the law of variation of viscosity is 
Mo = po(To) (2.41) 


where the nondimensional quantities—i.e., 
ture 7, velocity wo, etc., are in terms of certain constant 


tempera- 


quantities at x) = 0. 

Eqs. (2.37) and (2.38) can be integrated numerically 
by stepwise finite difference method starting from the 
initial conditions, Eq. (2.39), with the help of Eqs. 
(2.40) and (2.41). Any law of variation of viscosity 
may be used in this integration without introducing any 


complication. 


Stability in an Inviscid Fluid 


(1) The Differential Equation of Disturbances.—In the 


_ limit of infinite Reynolds Number, the differential equa- 


tion of @ becomes 








nd- 


33) 
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the 
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2 P Tl’ —2M*(w—cjw’ , J 
T’ — 2M*(w — c)w’ , an , w—-w . 
” - o+| T— Mw — 0)? =0 (2.42) 


T — M°*(w — c)* 


In order to remove the singular point w = c for neutral 
stable disturbance, we must have 

ka Leg (2.43) 
where 7,, = T(w = c,), ete., and c, is the velocity of 
propagation of the neutral stable disturbance in the in- 
viscid fluid. For given velocity and temperature dis- 
tributions of the main flow Eq. (2.43) determines the 
velocity of propagation of neutral stable disturbances, 
both symmetrical and antisymmetrical, in the inviscid 
fluid. In case of incompressible fluid 7,’ = 0, 7, 
constant, hence the condition reduces to 


WwW... = 0 (2.44) 


It is interesting to note that if Eq. (2.43) is satisfied 
the imaginary parts of A,,’s and B,’s are all equal to 
zero, provided that 2 is a real number—..e., for subsonic 
disturbances relative to the surrounding stream. 

For a given velocity distribution as Mach Number ./ 
increases, 2? decreases. If A/ is large enough such that 
for c = ¢, 2? < 0, there will be no subsonic neutral 
stable disturbance relative to the surrounding stream 
possible in the jet flow. The flow is then stable with re- 
spect to small disturbance for such large Mach Num- 
ber. 

(ii) Symmetrical Disturbance.—The determinantal 


equation for this case reduces to 
®, “= — (Bo Qa) + B, _— Boa + oe & QO (2.45) 


By = 7,2/c,* is real for subsonic disturbance in the 
jet relative to the velocity of the surrounding stream. 
If we put c = c, determined by Eq. (2.43) into Eq. (2.45) 
all the imaginary parts of B, are identical to zero; hence 


we have 
Bole.) — asBin(Cs) + a.*Bor(es) —...= 0 (2.46) 


Eq. (2.46) determines the wave number a, for the neu- 
tral stable disturbance which can be solved by succes- 
sive approximation. 
We see that as Mach Number J/ increases both Q(c,) 
and a, decrease. If .\/ is equal to the critical value 
Mee = VT,,/c.2 (2.47) 
a, = 0. As we shall show later a > a,, the disturbance 
is stable. Hence \/ > Mer, the jet-type flow will be 
stable with respect to all small symmetrical disturbances 
in the inviscid fluid. 
For the case of no energy transfer, the critical number 


is about 2.5. 


= £ 


a” TM 2) ' 
—— — M2(w —c)? 
7 t £ f 


If we vary a slightly from a,, in order to satisfy Eq. 
(2.45), ¢ must be different slightly from c,. Put 
a=a,+ Aa, c=c,+Ac=c,+ Acre +idc; (2.48) 


where Aa is a real number and Ac may be complex. 
Substituting Eq. (2.48) into Eq. (2.45) and neglectin 

g | | g g 
the higher order terms we have 


a Ac; 
vi a [OBy(c ) 0c]: Bir(c;) 


(= Bur)’ + ( oe)" (2 19 
~” ‘ a, <.td) 
| \ Oc ° Oc : Oc ee 


where B,,(c,) is a positive real number 


Aa 


[OB,(c.)/Oc),; = —x(T.,” /wWe,."*) (2.50) 


Both 7...” and w,’* are positive. Hence, | [0B,(c,)]/dc}, 
is a negative number. Hence, from Eq. (2.49) we see 
that Aa and Ac; have opposite signs. 

As we vary a from a, if a < a, c; being positive, 
the flow will be unstable with respect to such a dis- 
turbance; while a > a,, c being negative, the flow will 
be stable with respect to such a disturbance. Hence 
when the Mach Number is large enough such that 
a, = 0, the jet flow will be stable with respect to all 
small symmetrical disturbances. 

(111) Antisymmetrical Disturbances.—The_ determi- 
nantal equation for this case reduces to 


g,- = Ay — a.Ay + a°A, —...=0 (2.51) 


where Ay = 1. For subsonic disturbance relative to the 
velocity of the surrounding stream, A; is a positive 


real number and 


A, = Aar + 1Axn (n = 2) 
If we put c = c, given by Eq. (2.48), all the A,, 
vanish, and we have 
®,—© = 1 — a,Ai(c,) + a;*Aor(ce,) —... = O (2.52 


Eq. (2.52) determines the wave number a, for the neu- 
tral stable antisymmetrical disturbance. For a low 
Mach Number case, A; is in general larger than Ao,; 


for first approximation 


a, = 
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Here we find that for low subsonic main flow / < 1, a, 
for antisymmetrical disturbance in general increases 
with the Mach Number if the effect of energy 
transfer is small. Hence, for low subsonic flow with 
small energy transfer, the increase of flow velocity 
has a destabilizing effect on the flow. As M increases 
near the critical value M/¢r, A; approaches to zero, and 
Eq. (2.54) will no longer be a correct approximation for 
a;. The most suitable approximation will then be 

(A; — VA," — 4Aopg)/2A or (2.54) 


Qa = 
and as A; ~ 0 


a, 1/v —Aor (2.55) 


As Mach Number J/ is larger than M,, there will be no 
discrete eigenvalue a,, and all neutral disturbances will 
be supersonic disturbance. The flow is believed to be 
stable though it is not possible to show this as clearly as 
in the case of symmetrical disturbance. 

If we vary a slightly from a, in order to satisfy Eq. 
(2.52), c must be different slightly from c,. Put Eq. 
(2.48) into Eq. (2.52), and, neglecting the higher order 


terms, we have 
| ( OA ) OA | . 4 | (~ ‘) | 2 
7 Oc /; oc h. om Oc J Jc, 


] OA» 7 
a,\ OC Les A ae (Cs) 3 


(OA»/ Oc); is in general positive, A,, is a negative number. 
Hence, Aa and Ac have oppositive signs. Hence, the 
disturbance a > a, is damped, and a < a, is self- 


(2.56) 


Aa = Ac 


excited. 


Stability in a Viscous Fluid 


(1) General Considerations.—-In the case of incom- 
pressible fluid, we found that the viscosity has a stabi- 
lizing effect on the disturbance in inviscid fluid. Since 
the solutions of stability equation of compressible fluid 
are similar to those of incompressible fluid, it is be- 
lieved that the study of stability of inviscid fluid is suf- 
ficient for jet-type flow, and hence the computations for 
viscous fluid are not carried out. Of course it would be 
interesting if such computations were carried out. Such 
computations are much more complicated, for the sta- 
bility calculations must be performed for each value of 
Mach Number and of each condition of heat transfer. 
In the following we summarize the determinantal equa- 
tions for the viscous fluid. The results of numerical cal- 
culations will be reported later. 

The determinantal 


equation for symmetrical disturbance is 


(11) Determinantal Equations. 


(2.57) 


®, e + (] aR), C2)’ ee 0 


The determinantal equation for antisymmetrical 
disturbance is 
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i ; 
p +4 (b,—}F1 + W721 — c)?]F,(03 4 
aR 
F,(0) | + @—'") = Q (2.58) 


where F',(0) and F,(0) are the values of F, and F, at 
y = 0, respectively. 


. w"T Pee 7 
F = fy a ae M-? _— = 
(1 — c)° l—c 


Fy = (y/ov) [—mi(y — APT" + 2ui(0"/p) X 
{17/1 — o)?] — MP} = (Cy — DAPT (dus /dT)T") 


(2.59) 


CONCLUSIONS 


From the method of small perturbation, the following 
conclusions may be drawn for the stability of two- 
dimensional laminar jet-type flow of gas: 

(A) For incompressible fluid, we have that 

(1) Both the symmetrical and the antisymmetrical 
disturbances are unstable for large Reynolds Numbers 
but the region of instability of symmetrical disturbances 
lies inside that of an antisymmetrical one. 

(2)'The effect of viscosity on the neutral stable dis- 
turbance in inviscid fluid is always stabilizing for jet- 
type flow. 

(3) When the wave 
Reynolds Number is large, the jet flow is always un- 


number is very small and 
stable. 

(B) For compressible fluid we have that 

(1) The solutions of compressible fluid are similar to 
those of incompressible fluid. 

(2) So far as laminar stability is concerned an im- 
portant difference exists between the case of supersonic 
and subsonic disturbances relative to the velocity of the 
surrounding stream of the jet. For subsonic dis- 
turbance, discrete eigenvalues may be found and the 
limit of stability can be determined. For supersonic 
disturbance the eigenvalues are continuous and not 
discrete; no limit of stability can be found. Usually it 
is assumed that supersonic disturbances have no sig- 
nificance for the stability problem. But further research 
must be made to clarify this point. 

(3) For inviscid compressible fluid 

(a) Increase of Mach Number has stabilizing effect 
on the symmetrical disturbances. When Mach 
Number reaches the critical value, the jet-type flow 
will be completely stable with respect to small sym- 
metrical disturbances. 

(b) When Mach Number is sufficiently small, in- 
crease of Mach Number has a destabilizing effect on 
the antisymmetrical disturbance. But as Mach 
Numberissufficiently high, the wave number of neu- 
tral stable antisymmetrical disturbance increases 
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with Mach Number toward a limit. If the Mach F, = Fy + a’Fp (A.1) 
Number is larger than the critical value, no discrete 
value of wave number of neutral disturbance could ; . .i 

Fy = aFy + a'® Fx (A.2) 


be found. 
APPENDIX A Fry = Fa + a? Fy (A.3) 
tions p», and q:,.; for Compressible Fluid Case . oe ‘ ‘ : ss 
re P ¢ : P where the expressions of F;; are given in reference 15, 
The functions F;, F;,, and F;,;, of Eqs. (2.1), (2.2), pages 13-14 (except term —(7/p)(u)'w’ + ww") in Fy, is 
and (2.4) may be written as follows" missing there, which should be added). 


= ‘ 
bo = — / _ wlio + Mw — c)?]Fuo + Fao} dy (A.4) 
y w—c)?* 


> 


ok, T . T 
po = = / Mad Tal [1 + M?(w — c)? | Fire + F 20 | + Fy29 | dy + if ! — ur |x 
ye ~~ = ¢}* \ iw — ¢7* 
vm “ T *” (w — c)%pp 
dy | Fy) dy & ig / | ae u lay / = : P dy (A.5) 
y Ss (w — c)? J) 7 ; 


n2Q 
a 7 ° j 9 » . . . . { 
Pr» 4 = ( a iT ’ [1 aa M?(w —_ €)?7)[ Fuert 4 a“ Fy2(2n+2) | + Fs, 2n+4 F302, +2) 4 dy + 
e w— ¢j* 
1 =e M? | dy | Foy Qnt+2 + Fx (2n) ldy + ae M? dy a pe +2 dV A.6) 
a (w — c)? ° WS y (w — c)? J 7 : 
where F,,2.) = F,; with @ = (w — c)ho, and f, r, 7, @ calculated from this ¢ from Eqs. (2.15) to (2.18) 


Vo ii = 
= =_ / | _ Ta (1 + M?(w a o)7* Fin + Pent dy (A.7) 


(w— c)° 


, r ms 7m rT : 
oo = / | 7 ue | | 1+ M*(w — c)*| [Fig + Fis } + Fa33 + F330 dy + / —o | Ps 
Jy (w=— ¢€)* . w= < 


“~ ; : "” (w — cn 
dy Fy dy + -— M?|dy ———— iy (ae 
a ' » L(w — co)? Js 1 


n= 0 
Sa ‘A oi , oa i : ; ; — 
n+, = — 2 M? {1 + Al?(w — c)’ Firents + Fie 2n+3)| + Fy) n+5) f 2(2n+3)4 AV F 
. = CF J 
7 i T 7 Le 7 T : "” (w — c)? 
l = M- dy [Fy 2n+ 3 + Foe Int] |dy + a. M? dy mos Ge ay 
(w — c)? eS i J; (w — c)* . 7 
(A.9) 
where Fij,+1) = F;, with @ = (w — c)ke,+, and corresponding values of /, r, 7, 6 obtained from this @ from eqs. 


(2.15) to (2.18). 
APPENDIX B 


Definite Integrals in Asymptotic Expressions of ¢' of Compressible Fluid Case 


Ay = ] (B.1) 
QT, “i (w — c)? —c)? 

Die he J pat. I, ah, (B.2) 
(a — c)*J-« | 7 4 
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where 


° — 


a = Qo 


TQ 


7 (a — c)? 


f vo 
—o 


Bo 


B, = 


; 12 — ype. ” (w — c)? 
hy oes . Fy dy = - Po oa _ 
fa — c)? b (a — c)*? J-« 
= 7 9 j 9 ” . . ' 
Dy = ~ M? ; {1 + M?(w — 0)" | Fin + Fy} dy 
» L(w—c)° 
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Interrelationship Between Boundary Layer 
and Base Pressure 


H. H. KURZWEG* 


U.S. Naval Ordnance Laboratory 


SUMMARY 


Experiments conducted ia supersonic wind tunnels during the 
past 3 years (N.A.C.A. and N.O.L.) show that the pressure at the 
base of bodies is essentially a function of parameters that govern 
the boundary layer. Hence, the base pressure is closely con- 
nected to surface friction and heat-transfer phenomena. 

Experimental results obtained in the N.O.L. supersonic wind 
tunnels in a Mach range 1.5—5.0 on cylindrical bodies with and 
without boat-tails with conical heads under various systematic 
mechanical and thermodynamical variations of the boundary 
layer are presented and compared with theoretical results 


INTRODUCTION 


fo INVESTIGATIONS CARRIED OUT during the past 
years in supersonic wind tunnels and in free flight 
concerning the under pressure at the blunt or flat base 
of missiles and projectiles led to the conclusion that base 
pressure and boundary layer are closely interrelated. 
The present importance and status of the base-pressure 
problem can best be explained by answering the fol- 
lowing two questions: (a) What do we expect from 
base-pressure investigations, and (b) how far can avail- 
able base-pressure data be used by missile and projectile 
designers? 

The base drag caused by the base pressure is rela- 
tively small with respect to the total drag when the 
aerodynamic shape of the body is plump and its nose 
is blunt. It is, however, of the same order of magni- 
tude as wave and friction drag when the bodies are 
sharp-pointed and slender. Long-range projectiles 
and missiles must have these last features. 
the favorable shape, the total drag of large guided 
missiles can amount to several thousand pounds, of 
which the contribution of the base drag is then of the 
order of 1,000 or 2,000 Ibs. ! 

This is a heavy factor in the performance, trajectory, 
It would be helpful, 


In spite of 


and power-supply calculations. 
therefore, for the design engineer to have some formulas, 
or at least some reliable test data, relating to a number 
of basic bodies with the necessary directions for the ap- 
plication on a particular missile. 

In order to get such information, the Naval Ordnance 
Laboratory (N.O.L.) carried out systematic base-pres- 
sure measurements on cone-cylinder bodies of 5-cm. 
diameter in the 40- by 40-cm. Aeroballistics Tunnel over 
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a Mach Number range between 1.5 and 5.0 inclusive, 
as well as in its: Pressurized Ballistics Range with 30- 
caliber projectiles up to Mach Number 3.5 over a pres- 
sure range between 3 and 120 in. Hg. 

The various phases of these investigations are de- 
scribed in the following paragraphs. 


THE BASE PRESSURE AS FUNCTION OF MACH NUMBER 


In supersonic flow the dependence of the base pres- 
sure on Mach Number (1/) seemed to be of principal 
interest. As a matter of fact, all early base-pressure 
studies were made with the Mach Number as the only 
variable. 

If the available base-pressure data from Ames Aero- 
nautical Laboratory, National Advisory Committee for 
Aeronautics (A.A.L.-N.A.C.A.),'! N.O.L.,? Applied Phys- 
ics Laboratory, The Johns Hopkins University (A.P.L.- 
J.H.U.),* and Ballistics Research Laboratories (B.R.L.)* 
are plotted versus Mach Number, the chart on Fig. | is 


obtained. It can be seen that for each Mach Number a 
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Fic. 2.. Wake shapes of projectiles with laminar or turbulent 
boundary layer 
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wide region of measured values exists. In such a re- 
gion the base drag obviously can vary more than 50 per 
cent. 

Detailed investigations, however, soon showed that 
the test points could be separated into those points 
associated with laminar boundary layer on the bodies 
and those points associated with turbulent boundary 
layer. The turbulent data black 
signs) are more or less restricted to the shaded area 


(circles and_ short 
which represents the average scatter of the base-pressure 
ratio (ratio of base pressure to free-stream static pres- 
The 
laminar values (long vertical lines) show much larger 
variation. A.P.L.’s tests with rockets 
show especially low values in a Mach Number region 


sure) obtained in the B.R.L. atmospheric range. 
fin-stabilized 


between 1 and 2.5. Interference of body and _ fins 
might cause these low pressures. 


character of the boundary layer on the fins in this region 


Furthermore, the 


is not known. 
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This chart indicates clearly that the Mach Number 
is certainly not the only variable in the base-pressure 
function. Since the measurements and corresponding 
schlieren pictures (Fig. 2) indicate that the boundary 
layer is involved in the problem, the Reynolds Number 
must play an important part. 


PRESSURE AS FUNCTION OF THE FREE-STREAM 
REYNOLDS NUMBER 


BASE 


The Mach Number was held constant and the base 
pressure resulting from variation of the Reynolds Num- 
ber (Re) was measured. This variation was made in 
the N.O.L. wind-tunnel measurements by varying the 
lengths of the constant-diameter cylinders, as well as 
by reducing the diameter and length in proper scale. 
In the N.O.L. Pressurized Ballistics Range and in the 
A.A.L.' and Princeton University® wind tunnels, the 
Reynolds Number was varied there by changing the 
density of the air. One example of the wind-tunnel 
investigation at Mach 3.24 is presented in Fig. 3, where 
N.O.L. wind-tunnel data are plotted versus Reynolds 
Number. 

It can be seen in this figure that in the low Reynolds 
Number range, where the boundary layer is fully or 
partly laminar, a large variation of the base pressure 
occurs. It was determined experimentally that, by 
simple modifications of the contour of the model nose 
(roughness has similar effect), the base pressure can be 
changed more than 50 per cent. In this case, at a Rey- 
nolds Number of approximately 3.5 million, a constant 
value seems to be reached. Here, the schlieren pic- 
tures show wake shapes as associated with turbulent 
boundary layer, as does the picture (Fig. 4 left below) 
of the model with disturbances at the nose already at 
1.3 million. Large variations of the model surface in 
this constant region have relatively small effects on the 
base pressure. This apparent constant value is in 
agreement with B.R.L. firing-range data at the same 
Mach Number but at the much higher Reynolds Num- 
ber, 8 million. 

It is still debatable as to whether the turbulent value 
remains constant at extremely high Reynolds Number, 
as indicated by A.A.L.-N.A.C.A. measurements up to 
16 million and by A.P.L. rocket tests up to 100 million 
(7 = 1.5), or if a slight decrease takes place, as indi- 
cated in Princeton University’s wind tunnel, Fig. 6. 

Fig. 5 shows three more examples of Reynolds Num- 
ber dependence, at Mach 1.5 and 2.0, reported by Chap- 
man,! and at Mach 2.86 measured at N.O.L.? It can 
be seen in all three examples that the values become 
fairly constant as soon as turbulent boundary layer is 
developed whether or not this was initiated by artificial 
disturbance or by natural transition on the smooth 
models. 

There has been insufficient time to obtain as much 
base-pressure data as desirable in the N.O.L. Pressur- 
ized Ballistics Range since its recent completion. On 
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BOUNDARY LAYER A 


Fig. 6, therefore, are plotted single data points obtained 
in a Mach Number range between 3.00 and 3.28 at 
pressures between 3 and 120 in. Hg with 30-caliber, 25° 
cone-cylinder bullets (length diameter ratio L/D = 4). 
Here, the laminar and transition range seems to be ex- 
For 


tended to larger Reynolds Number (5-7 million). 
comparison, there are plotted the asymptote of the 
N.O.L. base-pressure curve for Mach 3.24 for fully de- 
veloped turbulent boundary layer, which is reached 
already at a Reynolds Number of 2.5 million, and the 
curve obtained in the Princeton University pressurized 
wind tunnel at slightly lower Mach Number (between 
295 and 2.98). Both wind-tunnel curves indicate a 
transition occurring at a much lower Reynolds Number 
than the free-flight data suggest. This might either be 
due to the larger cone angle (60° in the tunnel tests, 
25° in the firing range tests) or to the degree of dis- 
turbance in the wind-tunnel air. 

How much the shape of the body nose, the turbulence, 
and other irregularities of the air stream, such as shocks 
and pressure gradients, influence the transition quanti- 
tatively has not yet been sufficiently investigated. For 
short models (L/D = 2 and smaller) the variation of 
just the cvlinder length, as well as the alteration of the 
nose shape, is accompanied by complicated changes of 
Low base pressures are measured 
with extremely small 1 D ratios. This fact 
seen on Fig. 3, where at the low Reynolds Numbers the 
scale model (L/D = 


higher base pressures than the short bodies. 


the boundary layer. 
can be 


9.36) (denoted by x) has much 


Conclusions 


From the described data it must be concluded that 
wind-tunnel base-pressure measurements on models 
with fully developed turbulent boundary layer and cy- 
lindrical bodies that are not too short give, in general, 
results that are in agreement with free-flight turbulent 
data within +5 per cent. 

To predict base pressures of bodies with laminar or 
partly turbulent boundary layer seems to be much more 
difficult, since the surface can hardly be imitated on the 
small models. In spite of the same Mach and Rey- 
nolds Number being maintained, the data for a model 
can deviate considerably from full-scale data when acci- 
dental disturbances cause some degree of turbulence in 
the boundary layer. 

Uncertainty in base-pressure determination becomes 
still greater in compressible flow, at high Mach Num- 
bers, where temperature plays an important part in the 


boundary-layer characteristics. 


INFLUENCE OF MODEL AND BOUNDARY-LAYER 
TEMPERATURE ON BASE PRESSURE 


So far it can be seen that the Mach Number and the 
free-stream Reynolds Number do not constitute all the 
base-pressure parameters. The change of the velocity 
profile causes base-pressure variations in spite of con- 


ND BASE PRESSURE 745 


stant Mach and Reynolds Numbers. It should also be 
expected that temperature changes of the boundary 
layer and, in general, the temperature distribution in 
the boundary layer, characterized by Prandtl Number, 
alter the base pressure. The function p,/p = f(v, 7, 
T,) (bb = pressure, p = 
pressure) (with v as the velocity profile, 7 as the tem- 


base free-stream static 
perature profile, and 7°, as the profile of the total tem- 
perature of the boundary layer), therefore, was investi- 
gated with respect to the temperature profile. 

Two cone-cylinder models [one with laminar (L/D = 
3) the other one with turbulent (L/D = 10) boundary 
layer| were heated to approximately 400°F. and were 
cooled with dry ice to approximately —90°F. The re- 
sults for Mach 3.24 and for Re = 1.1 million and 3.5 
million, respectively, can be seen on Fig. 7a, where in 
both cases the base pressure is low at low model tem- 
perature (heat flowing from boundary layer to model) 
and increases with increase of model temperature. The 
effect is larger in the laminar than in the turbulent 
case, where the base pressure changes approximately 15 


per cent in this temperature range. Fig. 7b shows the 
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same effect at Mach Number 4.24. The 45-cm. model 
has in this case not a fully turbulent boundary layer. 

In Fig. 4 are four schlieren pictures, which show, at 
constant Reynolds Number and Mach Number, the 
various wake shapes due to differences in the velocity 
and temperature profiles of the boundary layer. 

Since in free flight at high Mach Numbers the hot 
boundary layer always loses heat by conduction to the 
missile from where part of it is radiated, we should ex- 
pect that wind-tunnel measurements give higher base 
pressures than free-flight measurements at correspond- 
ing Mach Number. The temperature of the wind- 
tunnel models is approximately the same as that of the 
boundary layer. 

Wind-tunnel measurements at Mach Numbers 4.24 
and 5.0 seem to point in this direction. 


THE EFFECT OF BOAT-TAILS ON THE BASE PRESSURE 


So far, the possibilities of predetermination of base 
pressures for new missile designs have been discussed. 
In this connection also, the influence of the shape of 
the small L/D ratio bodies on the base pressure was 
briefly mentioned. Furthermore, it has been found 
that cylindrical models show a large variation of the 
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base pressure when the tail is tapered and the boat-tail 
angle is varied. 

Fig. 8 shows results from N.O.L. measurements at 
Mach 3.24 with two cone-cylinder models (one with 
laminar and the other with turbulent, or partly turbu- 
lent, boundary layer) and conical boat-tails. The 
boat-tail angle was varied between 0° and 15°, and the 
boat-tail length was 5 cm. 
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BOUNDARY 


In both cases the base pressure increases with in- 
creasing boat-tail angle—in the laminar case to 7°, in 
the turbulent case to 15°. Fig. 9 shows that the lami- 
nar boundary layer is fully separated at 15°, whereas 
the turbulent boundary layer seems to be just de- 
tached at the same boat-tail angle. One would con- 
clude, therefore, that bodies with turbulent boundary 
layer may be fitted out with larger boat-tail angles. 


This advantage, unfortunately, is partly compensated 
by the increasing under pressure of the annular ring on 
the boat-tail, which gets larger and larger with increas- 
ing boat-tail angle. The net effect is the sum of base 
pressure times base area and of rim pressure times an- 
This pressure distribution on the tail 
Fig. 10 shows the re- 


nular-ring area. 
constitutes the total tail drag. 

sulting total tail-drag coefficient. 
in the turbulent case, as well as in the laminar case, an 
is useless in this 


It can be seen that 


increase of boat-tail angle above 7° 
example. 


A HeEvuRISTIC CONCEPT OF BASE PRESSURE 


The described phenomena suggest a theoretical con- 
cept that should help the better understanding of the 
base-pressure should make future 
calculations of the base drag possible. 


mechanism and 


It can be concluded from the experimental facts 
that the base-pressure problem is closely connected 
with boundary-layer characteristics, with surface fric- 
tion, and with heat transfer. Therefore, theoretical 
calculations must take in account not only Mach Num- 
ber and free-stream Reynolds Number but also the 
parameters of the boundary layer. 


If one makes, as a first approximation, the crude 
assumption that the velocity of the air which is re- 
duced in the boundary layer has to be restored after the 
boundary layer has separated from the body, the pres- 
sure at the end of the necessary expansion process can 
be calculated. The calculation, which requires the 
knowledge of the velocity and temperature profile and 
the stagnation-temperature distribution in the bound- 
ary layer, is based on the following function: 


7 ys 
Vy Py P 
Po Jy vp p 


— 
4 Vy Py 
| ydy 
vp 


where v, is the velocity and p, is the density in the 
boundary layer at a distance y from the model sur- 
face; the boundary-layer thickness is 6, and the body 
diameter is 27. (The letters without subscript denote 
free-stream conditions.) The pressure p; would be 
that of the layer y after the expansion if the air of the 


layer y could expand independent of the other layers y 
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from the stagnation pressure p,) and regain its original 
velocity v. Schlieren pictures indicate that this veloc- 
ity is almost regained at approximately 2 or 3 di- 
ameters behind the base. 

The pressures p; have to be connected with the vari- 
ous masses of the layers y. The integrated products of 
pressure and mass have to be divided by the total 
amount of air contained in the boundary layer. 

Replacing the density profile p,/p by the temperature 


£ 
Po 0 


v,T | 


Fully realizing that this concept can give only ap- 
proximately correct values, especially since the bound- 
ary-layer profiles are not sufficiently known (in par- 
ticular at higher Mach Numbers), we calculate this 
quantity. Fig. 11, for example, shows the evaluation 
for the base pressure of a cylindrical body with a bound- 
ary-layer thickness 6/r = 0.2 and a Prandtl Number 
Pr = 1. The velocity profile is approximated by an 
exponential function v,/v = (y/6)'/", with m taken for 


several values between | and 7. 

The average experimental data, represented by the 
curve, are located between the calculated 
curves. The effect of Prandtl Number not equal to | 
is estimated at Mach 5, resulting in higher base pres- 


dashed 


sures. 


RESULTS 


The results of base-pressure measurements in wind 
tunnels and firing ranges show that the base pressure is 
not only a function of the Mach Number and free- 
stream Reynolds Number but is also dependent on the 
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profile 7',/7° and expressing Pp, p by the velocity term 
(given by adiabatic expansion), 


2 eae l Pw 


where a,) = sonic speed at the stagnation temperature 


of layer vy and py = total pressure of the layer y, we ob- 


tain 


y dy 


6 Vy ‘a 


vi 


y dy 


specific boundary-layer conditions of the body as given 
by surface friction and heat transfer. Velocity and 
temperature profiles must be known before correct pre 
dictions of the base drag can be expected. Experi 
mental data of boundary-layer profiles, especially at 
the higher Mach Numbers, are still missing. At pres 
ent, only base-pressure data obtained by measurements 
with fully developed turbulent boundary layer can be 
used for full-scale base-drag estimations with reason- 


able accuracy. 
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Shear Center for Thin- 


Walled Open Sections 


Bevond the Elastic Limit 


G. H. HANDELMANt 


Carnegie Institute of Technology 


SUMMARY 


The location of the shear center of a thin-walled open section 
stressed beyond the elastic limit is determined by methods similar 
to those used in the elastic case. Through an analysis of the 
stress and strain components, the problem is reduced to determin 
ig one constant in the case of symmetric sections. This can be 
found from the magnitude of the applied force. The location of 


Calculations for channel 


shear center can then be found 
ious are carried out in detail, and the influence of the load and 
Formulas are 


tress-strain law is shown in a numerical example 
The shear center, in the 


also given for nonsymmetrical sections 
nonelastic case, is not a property of the section alone but depends 


the stress-strain law and the applied force. 


on 


INTRODUCTION 


HE PROBLEM OF DETERMINING the shear center for 
T. thin-walled open section has been thoroughly 
treated for those cases in which the stresses remain be- 
low the elastic limit.'. The present discussion is an 
attempt to extend this analysis beyond the elastic 
The question will be considered as one of ‘‘tech- 
rather than a problem in the exact 


limit. 
nical plasticity”’ 
theory of metals stressed beyond the elastic limit. 
Various simplifying assumptions, of the order of those 
appearing in the theory of elasticity used in studying 
strength of materials, reduce the complexity of the 
mathematical analysis considerably. It is hoped that 
these results will give at least a good first approximation 
to the phenomena under study. In particular, a thin- 
walled open section with one axis of symmetry, clamped 
at one end and loaded with force F at the free end acting 
perpendicularly to the symmetry axis will be considered 
first (see Fig. 1). This beam can be conveniently de- 
scribed in terms of a right-handed rectangular system 
of axes (x, ¥,. 2) with the origin 0 located in the cross 
section at the free end. 
along the axis of symmetry of the cross section, the x- 
axis is taken parallel to the axis of the beam, and the 
y-axis is perpendicular to both the x and z axes. The 
shear center s:(O, 0, e) is the point at which F must be 
applied such that the torques produced by the shear 
flow, on one hand, and the force F, on the other, will be 


In the elastic case, the coordinate e is 


The z-axis will be directed 


in equilibrium. 
independent of F and the elastic constants of the ma- 
terial; the present paper will consider the change in e 


Received August 9, 1950) 
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with F when the material is stressed beyond the elastic 
limit. 
fa 


A second system of coordinates, consisting 
reference line C in the cross section, the normals to C, 
and the x-axis, will be found useful. In this system, a 
point is given by the arc length s along C measured from 
a fixed point O, the normal distance n from C, and the 
coordinate x. The wall thickness /(s) is measured along 


the normals to C and is not a function of x. Precise 
statements as to the assumption that the wall thickness 


is small will be given in later sections. 


THE EQUATIONS OF EQUILIBRIUM 


Only two stresses are assumed to act in the beam —an 
axial stress o in the x-direction and a shear-stress 7 
parallel to the reference line C. This stress distribution 
will not fulfill the boundary condition of vanishing 
stress on the surface of the beam unless the inner and 
outer walls are also parallel to C. If the thickness is 
small, however, and the walls are almost parallel, the 
normal component will be small compared to the tan- 


gential component. Furthermore, 7 and o may be as- 
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sumed to be uniformly distributed across a normal to C 


if t(s) is small. Under these conditions, 


r= 7(s, x), o = a(s, x) 


The shear flow, T(s, x), is defined by the relation 
T(s, x) = t(s)r(s, x) 


The equations of equilibrium can be obtained by con- 
sidering the stresses acting on an element of a beam with 
sides ds, dx, t(s) as shown in Fig. 2. For equilibrium in 
the x-direction, 


T dx + tads — tla + (Oa/0x) dx] ds — 
[T + (O07 /ds) ds] dx = 0 


or 
(O7/d0s) + [t(00/dx)] = 0 (1) 
In the s-direction, 


T ds — |T + (OT/0dx) dx| ds = 0 


or 
oT /dx = 0 (2) 


Eqs. (1) and (2) constitute the equations of equilibrium 
which must be fulfilled by the stresses 7 and o. In 
particular, Eq. (2) implies that the shear flow is a func- 
tion of s alone. 


INTRODUCTION OF THE STRESS-STRAIN LAW 


It has been found experimentally that when a pris- 
matic beam is loaded in the direction of one of the 
principal axes of inertia of the cross section, the cross 
section will remain plane during bending.? Under this 
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hypothesis, the axial strain ¢ will be a linear function of 
y—that is, 


€ = yelx) (3) 


where the function g(x) is still to be determined. If the 
dimensions of the cross section are small compared with 
the length of the beam, it is well known in the elastic 
case that the shear strain is small when compared with 
the axial strain.*? The same assumption will be pre- 
served in the case treated here, in the sense that the 
axial stress o is supposed to be only a function of the 


axial strain e—that is, 


o = o(s, x) = E(e) = Ele(s, x)] (4) 


The function £ can be evaluated from a tensile test. In 
the elastic case, 


E(e) = Eve 


where /y is Young’s modulus. Since most materials be- 
have almost the same in compression and tension, /(e) 
will be considered as an odd function of «. 

Eq. (2) can be integrated directly with the result 


it Oa(s, x) ’ 
_ t(s) ds + f(x) 
0 Ox 


where f(x) must be determined from the boundary con- 
ditions. Since the edges of the section are free of stress, 
7T(0, x) = 0. Consequently f(x) must vanish also, and 


: 5 Oa(s, x) ? OE 
T(s,x%) = — t(s) ds = — Ks) 1s 
(sx) f (s Ox f ts Ox 


(5) 
According to Eq. (2), 7(s, x) is a function of s alone; 
from Eq. (1) then 


T(s,. x) = 


0°0/0x* = 0 


Consequently, o(s, x) must be a linear function of x; or 
in other terms, 


a(s, x) = Ele(s, x)] = Elye(x)] = Ai(s)x + Ad(s) 


(6) 
The two functions, A;(s) and A:2(s), can be found 


from the condition of equilibrium for the whole beam. 
Let M = M(x) be the bending moment at any section x 


and / the total length of the reference line C. Now M = 
— Fx, but 
l 
M=- f y(s)t(s)o(s, x) ds 

0 
Thus, 

l 

Fx = j y(s)t(s)a(s, x) ds (7) 
0 


and differentiation with respect to x yields 
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" Oa(s, x) l 
F= y(s)t(s) ds = y(s)t(s)Ai(s) ds 
0 Ox o 


Now for x = 0 
A.(s) = E[y(s) ¢(0)} 
and for x = 1 


A,(s) + Ax(s) = E[y(s)¢(1)] 
or 


Ax(s) = Efy(s)e(1)] — Efy(s)¢(0)] 


Eq. (7) can thus be rewritten in the form 


Fx =x / y(s)t(s) } ELy(s)e(1)] — Ely(s)¢(0)]} ds + 


| 
/ y(s)\i(s\E 
0 


In other terms, the functions A,(s) and A2(s) can be de- 
termined if two constants ¢(0) and ¢(1) are found such 


that 
l 
/ y(s)t(s)E[y(s)e(1)] ds = F 
0 


l 
f y(s)t(s)E[y(s)e(0)] ds = 0 


The function y(s) must be antisymmetric and ¢(s) sym- 
metric about s = //2 from the assumed symmetry of 
the cross section. Since E[y(s)¢(x)] is an odd function 





y(s)¢(0)] ds 





o (ib. in.-2) 
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of ¢, it will also be antisymmetric about s = //2. Fur- 
thermore, f(s) = 0; therefore, 


l 
f y(s)t(s)Ely(s) ¢(0)] ds 
0 


is strictly positive or strictly negative according to 
whether ¢(0) > 0 or ¢(0) < 0. It will vanish only if 
¢(0) = 0; thus, ¢(0) = O and A.(s) = 0, Ai(s) = 
E[y(s)g(1)]. The stress distribution can be completely 
determined then by finding one constant ¢(1) such that 


l 
f y(s)t(s)ElLy(s)e(1)] ds = F (8) 


In practice, it is somewhat simpler to assume various 
values for g(1) and compute F rather than apply the 
direct method. Since 


OF /dx = Ai(s) = Ely(s)¢(1)] 


the shear flow 7 can be found from Eq. (5)—namely, 


T=- [ t(s)E[y(s)e(1)] ds 
0 


Os 


LOCATION OF THE SHEAR CENTER 


The shear center e can now be determined by com- 
puting the torque generated by the shear flow T. This 
torque must be equal in magnitude and oppositely 
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directed to the torque Fe produced by the force F ap- 
plied at the shear center. Since the shear force is as- 
sumed to be uniformly distributed along the normal to 
the reference line, the torque due to the shear forces can 
be found by replacing these shear forces by the shear 
flow concentrated on the reference line. The torque Q 
is given by 


*] 
Q= / r(s)T(s) ds, 
« 0 
where 7(s) has been defined in the first section. There- 
fore, 
l 
Fe = — f r(s)T(s) ds 
0 
or 
a a 
e=--— y(s)T(s) ds = 
Ff 0 


*! "4 
/ r(s) if t(é)E[y(E) e(1)] dé } ds 
7 = 0 
/ y(s)t(s)Ely(s)¢(1)] ds 


Consequently, when the constant ¢(1) has been de- 
termined for a given value of F according to Eq. (8), 
the shear center corresponding to this value of the ap- 
plied force can be found from Eq. (9). 

In the elastic case, E[y(é)¢(1)] = Eoy(é)¢(1), and Eq. 
(9) becomes 


(9) 


/ yv*(s)t(s) ds 
0 
| a 7s 
] | | f vieye) a | ds (10) 


where /, is the polar moment of inertia of the cross 
section about the z-axis. Eq. (10) contains the known 
result that, in the elastic region, the shear center 1s inde- 
pendent of the material—i.e., independent of )—and 
the applied force; it is only a function of the geometry 
of the cross section. ' 


APPLICATION TO CHANNEL SECTIONS 


The computation of the shear center as a function of 
the applied force, Eq. (9), is relatively simple in the im- 
portant case of a channel section. Consider, for ex- 
ample, the section shown in Fig. 3 with the dimensions 
as indicated there. The origin of coordinates can be 
taken most conveniently at the center of the web. The 
constant ¢(1) can be found from Eq. (8) in the following 
way. Since y(s) = b for 0 S s S a and the section is 


symmetric about the z-axis, Eq. (8) becomes 


AERONAUTICAL 
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F = 2abt\E 





be(1)] + 


a+b 
mf (a+b—s)Ella +b—s) elds (1 
a 


Now the integral appearing in Eq. (11) can be simpli- 
fied, for 


a+b 
[ (a + b— s)El(a +h s)e(1)} ds = 
a 





=~at+b-—vs. 


dimensions of strain; it is convenient to set 7 = fy(1 


Then 
I *by(1) 
: / nE(y) dn 
¥' ] ia 0 


where & 


*b 
| tE[Eg(1)] dé = 


Furthermore, let n* = bg(1); Eq. (11) then becomes 


a fT . 
F = 2abt; E(n*) + —,, / n(n) dn (12) 
i 0 


The integral appearing in Eq. (12) has the property 
that it is independent of the dimensions of the channel 
section. It can be evaluated as a function of the upper 
limit, n*, once for all channel sections made of the same 
material provided that the stress-strain curve in tension 
is known. This integral has been evaluated for the 
stress-strain curve of Fig. 4, and the results are plotted 
in Fig. 5. For a section of given dimensions, a and 3, it is 
more convenient to compute F as a function of 7* from 
Eq. (12) rather than find the value of n* which corre- 
sponds to a given value of F. The quantity 7* is a 
useful parameter in this connection as well as in the 
subsequent evaluation of the position of the shear cen- 
ter, e, asa function of F. 

The equation defining the shear center, Eq. (9), is 
also somewhat simplified in the case of channel sections 
since r(s) = bforO S s Saandr(s) = Ofora<s<at 
b. Since 


a+b, 


the integrand is symmetric about 5s 


9 "a *s abt, fe(n* ) 
of b | / ti E(n* us| ds = : 13 
F Jo 0 F 


Eq. (13) can be most readily evaluated by considering e 
as a function of n* and using the corresponding values of 
F found from Eq. (12). The only difficulty in the com- 
putation occurs when F = 0; then 7* vanishes also and 
the expression for e is indeterminate. This problem is 
only superficial, however, since the material is elastic 
for small values of F and /(n) = 


Kon. In this case, 


F = 2abt;Eun* + (2/3)b%t2Eon* 


and 
(a7bt;) 


e= Pabt, + (2 3)b7t.] 14) 


By definition §g(1) has the 
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SHEAR CENTER FOR 


As previously remarked this equation is independent of 
Fand Eg; it can be used to find the position of the shear 
center for small values of the applied force. 

These results have been applied to the case of a sec- 
tion with dimensions a = } = 1.8in., 4) = te = 0.15 in. 
The final curve of e vs. F for this example is shown in 
Fig. 6. It is interesting to note that the change in the 
position of the shear center is less than 10 per cent of the 
minimum distance of the shear center from the center 
For small values of F, the function e is 
this result arises from the fact 


of the web. 
practically constant; 
that within this range of F the material behaves elas- 
tically and the shear center can be found essentially 
from Eq. (14). As the applied force is increased be- 
yond this range, the position of the shear center shifts 
toward 0, reaches a minimum distance, and then starts 
to increase. This increase is caused by strain hardening 
in the material in the sense that the stress-strain curve 
(Fig. 4) is almost linear for strains larger than 0.01. If 
the material could be strained indefinitely without 
failure and if the basic assumptions were valid for such 
strains, the shear center would return to the position 
given by the elastic theory. Roughly speaking, for 
large strain values the nonlinear portion of the stress- 
strain curve has little influence, and the material be- 
haves as an elastic body with a different Young’s 
modulus. For an elastic body, however, the position of 
the shear center is independent of Young’s modulus; 
consequently, the limiting value of the shear center will 
be independent of the tangent modulus in the strain- 
hardened range. This result can be found analytically 
by replacing E(n) by Lyn + k for n > m in Eqs. (12) and 
(13). If n* is then allowed to increase indefinitely, Eq. 
(14) will arise as a limiting case. This limiting case 
should not be taken too seriously since it exceeds the 
working range of the material and the original assump- 
tion of small deflections. On the other hand, such con- 
siderations explain the final increase in the distance of 
the shear center from 0. 


EXTENSION TO NONSYMMETRICAL SECTIONS 


The previous analysis can be extended to the problem 
of determining the shear center of a nonsymmetric 
thin-walled Unfortunately, the resulting 
numerical calculations for a special case are not nearly 


section. 


so simple as those of the preceding discussion or the 
equivalent elastic problem. This difficulty is caused by 
the essential nonlinearity of the problem. If the applied 
force has components F, and F,, the problem for the 
elastic state is solved roughly in the following way.° 
The shear flow 7, corresponding to the force F, alone 
can be found by a technique similar to that used in 
Section 4 specialized to an elastic stress-strain law. 
The requirement that the torque produced by F, be in 
equilibrium with the torque produced by the shear flow 
T, will determine a line y = const. along which F, must 
act. Similarly, the shear flow 7: can be found corre- 
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sponding to F,; and a similar line, z = const., along 
which F, Since the shear 
flow 7’ which is produced by the force with components 
F, and F, is given by T = T, + 7), the intersection of 
these two lines will determine the shear center. In the 
nonelastic case, however, the shear flows can no longer 


must act, can be determined. 


be linearly superposed; the effects of both components 
F, and F, must be studied simultaneously. As will be 
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seen presently, the resulting numerical work is much 
heavier than that encountered in the previous section. 

Consider an open thin-walled prismatic beam as 
shown in Fig. 7. Coordinate systems similar to those 
used in the preceding section will be used. Let G, the 
centroid of the undeformed cross section, be taken as 
the origin of a rectangular Cartesian coordinate system 
with the y and z axis lying in the plane of the cross sec- 
tion and the x axis parallel to the generators of the 
cylinder. Again r(s) will denote the distance from a 
fixed point 0 to a tangent to the reference line C at the 
point s where, as usual, s denotes the arc length meas- 
ured along C. Asis customary in the theory of strength 
of materials, the parallel to the generators passing 
through G will be assumed to be the neutral fiber. In 
addition, the cross sections are supposed to remain 
plane. Consequently, the axial strain « must be of the 
form 


€ = yo(x) + 2H(x) 


where the functions g(x) and (x) are still to be deter- 
mined. According to the stress-strain law, the axial 
stress o can be written as 


o = E(e) = Elye(x) + sy(x)] 


Now the equilibrium equations, Eqs. (1) and (2), are 
independent of whether or not the section is symmetric. 
Consequently, 

0*a/0x? = 0 
or 


o = xA;(s) + Ad(s) (15) 


The functions A;(s) and A2(s) are still to be found. 
Since o = E[ye(x) + 2z¥(x)], 


Ax(s) = Elye(0) + 2¥(0)] 
A,(s) = E[yg(1) + 2¥(1)] — Elye(0) + 2¥(0)] 


On the other hand, the beam is loaded only in the plane 
of the cross section at x = 0; consequently, the axial 
stress o must vanish for x = 0. Therefore, Ao(s)= 
E[y¢(0) + z¥(0)] = 0 for all values of y and z. Since 
the function E vanishes only when the argument is 
zero, this implies that ¢(0) = ¥(0) = O and 


Ax(s) = Elyg(1) + sy(1) ] 


The bending moment with axis in the z-direction is 
given by 


l 
—-Fyx = - [ y(s)t(s)o(s, x) ds = 
0 
l 
= ri y(s)t(s)xAy(s) ds 
0 


where / is the total length of the reference line. 
larly for the moment in the y-direction, 


Simi- 
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l l 
Fa = [ 2(s)t(s)o(s, x)ds = Z z(s)t(s)xAy(s) ds 
Consequently, 
1 
Ff, = rf y(s)t(s)A4(s) ds 
0 


l 
n= f 2(s)t(s)Ay(s) ds 
0 


l 
[ y(s)t(s)E[ye(1) + 2¥(1)] ds 


l | 
Fi 2(s)t(s)E[ye(1) + 2y(1)] ds 


Given values of F, and F,, the constants g(1) and (1) 
can be found by solving Eqs. (16). The problem of de- 
termining ¢(1) and y(1) is thus reduced to solving two 
simultaneous equations that, in general, are not linear 
in g(1) and ¥(1). This is one of the basic numerical dif- 
ferences between the nonsymmetric and the symmetric 


or 


>) 
>| 
Il 


(16) 


a 
II 


problem in which only one constant must be found 
[see Eq. (9)|. Should it be desirable to find the shear 
center for a continuous range of /, matched with all 
values of a continuous range of F,, the computation can 
be reversed by picking a suitable set of values of ¢(1) 
and ¥(1) and determining the corresponding F, and F, 
by means of Eq. (16). 

Assuming that the constants ¢(1) and ¥(1) have been 
found, the position of the shear center can be deter- 
mined by a procedure analogous to that used in the 
symmetric case. According to Eqs. (1) and (15), the 
shear flow 7(s, x) is given by 


a "e Oa(s, x) “s 
T(s,x) = - t(s) ds = — t(s)A,(s) ds 
0 Ox 0 


= — fi t(s)E[yg(1) + 2y(1)] ds 


Let e, and e., respectively, be the y and z coordinates of 
the shear center with respect to a fixed point 0 in the 
plane of the cross section. Then the torque Q produced 
about 0 by the forces F, and F, is 


Q = —Fe. + Fe, 


On the other hand, 


*l 
o- f r(s)T(s) ds 


therefore, 


Fie, — Fe, = 


l 8 
[ r(s) if tE)Ely(E)e(1) + s(e)v(1)] de ds (17) 


Eq. (17) then determines the line of action of the force. 


(Continued on page 766) 
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The Effect of Flight and Configuration 
Variables on Thermal Stresses in 
Diamond-Shaped Supersonic Wings 


FRANKLIN P. DURHAM? 
Unwersity of Colorado 


SUMMARY 


For aircraft traveling at high supersonic speeds, approximate 
stagnation temperatures exist in the boundary layer surrounding 
the aircraft. If the aircraft has accelerated rapidly to a high 
supersonic speed, large transient temperature gradients will be 
introduced, which may, in turn, cause high thermal stresses. 

Methods for the determination of the transient temperatures 
and the stresses resulting from these temperatures are presented 
for solid diamond-shaped airfoils. Numerical results are given, 
showing the effects of several flight and configuration variables 
based on a wing of given specifications. 

It is shown that the maximum thermal stresses occurring in the 
chosen wing are compressive stresses at the leading edge. Itis 
further shown for the severe flight conditions assumed that these 
stresses are of large magnitude. 


SYMBOLS 


Following is a list of the principal notations used: 


@ = thermal diffusivity, k/we; 

c = chord length, ft. 

G4 = specific heat of wing material, B.t.u. per lb.°F, 

Cp = specific heat of air, B.t.u. per lb.°F. 

E = modulus of elasticity, Ibs. per sq.in. 

£ = acceleration due to gravity, 32.2 ft. per sec.? 

h = convective heat-transfer coefficient for air, B.t.u. per 
ft.2 sec. °F. 

J = mechanical equivalent of heat, 778 ft.lb. per B.t.u. 

k = thermal conductivity of wing material, B.t.u. per ft. 
sec, °F, 

k, = thermal conductivity of air, B.t.u. per ft. sec.°F. 

j = beamwise thickness of wing at any chordwise station, 
ft. 

M = Mach Number 

p = air pressure, lbs. per sq.ft. 

r = recovery factor, dimensionless 

t = maximum beamwise thickness, ft. 

T = average beamwise temperature at any chordwise 
station, °R. 

T, = initial temperature of wing, °R. 

Taw = adiabatic wall temperature, °R, 

Tm = mean linear chordwise temperature variation that 
would produce deflection but no stress, °R. 

Ty, = mean linear beamwise temperature variation that 
would produce deflection but no stress, °R. 

free-stream temperature, °R. 


Io = 


Received May 10, 1951. 

* This paper represents a portion of the work carried out for the 
U.S. Air Force under Contract No. W33-038-ac-17240 and is 
an advance release of a portion of U.S.A.F. Technical Report No. 
6351. The project was directed by Prof. H. W. Sibert. Members 
of the University of Colorado Engineering Experiment Station 
staff, who assisted in the work, are listed in references 3 to 16. 

t Assistant Professor of Aeronautical Engineering. 


T; = stagnation temperature, °R. 

T, = temperature at any beamwise point, °R. 

u = speed, ft. per sec. 

w = specific weight of wing material, Ib. per ft.* 

x = distance along chord, measured from leading edge, ft. 
y = distance along span, ft. 

2 = distance in beamwise direction, ft. 

Nu = Nusselt Number 

Pr = Prandtl Number 

Re = Reynolds Number 

a = linear coefficient of thermal expansion, ft. per ft.°F. 
6 = time, sec. 


4“ =6 =-—i Viscosity 
Poisson’s ratio, ft. per ft. 
thermal stress due to combined chordwise and beamwise 


x 
Ml 


%Yy = 
temperature variation, lbs. per sq.in. 

oy» = thermal stress due to beamwise temperature variation 

Cyc = thermal stress due to chordwise variation of average 
beamwise temperature 

T = (Taw =r T)/( Taw = Ta) 


INTRODUCTION 


OR AIRCRAFT TRAVELING at high supersonic speeds, 
F approximate stagnation temperatures exist in the 
boundary layer surrounding the aircraft. For example, 
the boundary-layer temperature at a Mach Number of 
4 at sea level is approximately 2,000°F. The effect of 
this high temperature boundary layer is to heat the 
material that it surrounds, primarily by forced convec- 
tion. 

If the aircraft has accelerated rapidly to a high super- 
sonic speed, large transient temperature gradients will 
be introduced, which may, in turn, cause high thermal 
stresses. If the acceleration is slow, the aircraft will 
heat more or less uniformly, causing thermal expansion 
but only small thermal stresses. 

The problem of thermal stresses due to aerodynamic 
heating is particularly important in supersonic wings, 
which are thin and thus subject to being heated quickly 
and, moreover, are highly stressed due to the aerody- 
namic loading without the addition of thermal stresses. 

Before thermal stresses can be found, the temperature 
distribution in the wing under consideration must be 
obtained. The determination of the temperature dis- 
tribution, in general, involves the solution of three- 
dimensiona! heat flow in the transient state. 

This report contains the results of an investigation of 
temperature distributions and thermal stresses in a 
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wing in supersonic flight. The program undertaken 
was composed of two parts, analytical and experimental. 
Analytical studies of the various problems involved 
were made. These studies were followed as far as 
possible by experimental tests to verify the theoretical 
work. 


In order to reduce the number of variables, a solid, 
thin, diamond-shaped wing was taken as the basis of 
study. This basic wing was considered to be initially 
at uniform temperature and then instantaneously sub- 
jected to a steady supersonic flight Mach Number. 
The instantaneous attainment of the flight Mach Num- 
ber would impose the most severe thermal stresses pos- 
sible for the given flight Mach Number and is thus the 
limiting condition approached by rapid acceleration of 
an aircraft or missile. 


The importance of the various types of heat flow 
chordwise, beamwise, spanwise, and the combinations 
of chordwise and beamwise and chordwise and span- 
wise—were investigated in order to obtain the simplest 
solution of reasonable accuracy. Investigations were 
also made of the effects of flight and wing-configuration 


conditions on thermal stresses. 


Experimental tests were conducted to verify the re- 
sults of the analytical investigations insofar as possible. 
The experimental program consisted of simulating the 
conditions of heat input to model wings and determining 
temperature distributions and thermal stresses. 


The description and results of these individual 
analytical and experimental investigations are listed as 


references 3 to 16 inclusive. 
Basic CONDITIONS AND ASSUMPTIONS 
Wing 


The basic wing used in these studies had the following 
characteristics. 


(1) Configuration 

(a) Solid, diamond-shape with '/j¢-in. leading 
and trailing edges 

(b) Chord, 4 ft. 

(c) Span, 16 ft. 

(d) Maximum thickness, 5 per cent of chord 
(0.2 ft.) 

(e) Taper, none 


(2) AMaterial—mild steel 

(a) Specific weight, 492 Ibs. per ft.* 

(b) Specific heat, 0.122 B.t.u. per lb.°F. 

(c) Thermal conductivity, 0.0072 B.t.u. per ft. 
sec. °F. 

(d) Coefficient of thermal expansion, 6.5 X 10~§ 
in. per in.°F. 

(e) Young’s modulus, modulus of elasticity in 
tension or compression, 28,000,000 Ibs. per 
sq.in. 
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Flight Conditions 


Although the method adopted for the determination of 
temperature distribution and thermal stresses is a step- 
by-step process and applies equally well to flight with 
finite acceleration, the most severe flight conditions 
That is, the 
wing under consideration was assumed to be initially at 


were chosen for all problems analyzed. 


a low uniform temperature and then instantaneously 
subjected to a steady supersonic Mach Number 
Simplifying Assumptions 


The following assumptions apply to all solutions ob- 
tained. 
to particular problems and will be listed in the sections 


Additional assumptions were made with regard 


where the problems are described. 

(1) Heat transfer by radiation to or from the wing 
surface is negliglble. 

(2) The specific heat, specific weight, thermal con- 
ductivity, coefficient of thermal expansion, and the 
modulus of elasticity of the wing material are constant. 

(3) The wing does not buckle and is not restrained. 


Heat-Transfer Coefficient 


Considerable uncertainty exists in the evaluation of 
convective heat-transfer coefficients in supersonic flows. 
The following relations have been widely used for pre- 
dicting supersonic heat-transfer coefficients: For lami- 
nar boundary layers 


Nu = hx/k = 0.33Re,°Pr'” (1 
and for turbulent boundary layers 
Nu = hx/k = 0.030Re,’’Pr'’* (2 


in which the Reynolds and Nusselt Numbers are 
evaluated in terms of the distance x from the leading 
edge. 

The choice of the air temperature in evaluating Eqs. 
(1) and (2) is subject to considerable question for super- 
sonic flows and has been the subject of numerous in- 
vestigations. The air temperature at the wing surface 
has been used as well as the temperature in the free 
stream just outside the boundary layer. 

References 1 and 2 recommend the use of the surface 
temperature. In the transient state, the surface tem- 
perature varies with time, consequently the heat- 
transfer coefficient based on the surface temperature 
would vary with time. In order to simplify the prob- 
lem, it was felt desirable to use a heat-transfer co- 
efficient that was constant with time. For this reason, 
the so-called “adiabatic wall’’ temperature, which all 
surface temperatures approach with increasing time, 
has been used throughout this report. This temperature 
will be referred to as the adiabatic wall temperature. 

In the early time of the transient state near the lead- 
ing and trailing edges, the surface temperature ap- 
proaches the adiabatic wall temperature while the sur- 
face temperature near the mid-chord is much lower than 


the adiabatic wall temperature. Consequently, the re- 
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FLIGHT AND CONFIGURATION 


sult of using the adiabatic wall temperature instead of 
the surface temperature is that values of heat-transfer 
coefficient near the leading and trailing edges will be 
approximately the same as those based on the surface 
temperature, while, in the region near the mid-chord, the 
values will be appreciably lower than those based on 
the surface temperature. 

On the other hand, the use of the free-stream tem- 
perature just outside the boundary layer will result in 
values of heat-transfer coefficient near the leading and 
trailing edges which are much greater than those based 
on the surface temperature, while, in the region of the 
mid-chord, the values will also be somewhat greater 
than those based on the surface temperatures. 

These variations are due to the fact that the adiabatic 
wall temperature is approximately the stagnation tem- 
perature, while the free-stream temperature outside the 
wing boundary layer is not greatly different from the 
ambient temperature in front of the wing (even though 
oblique shock waves are attached at the leading edge). 

The value of heat-transfer coefficient is important in 
the determination of transient-state temperature dis- 
tributions and thermal stresses, and, in light of this dis- 
cussion, the numerical results presented here may be 
subject to change when more is known about the de- 
termination of supersonic heat-transfer coefficients. 

The transition from laminar to turbulent boundary 
layers for supersonic flows is believed to occur at a 
Reynolds Number of about 500,000.' For all problems 
considered in this investigation, the transition would 
then occur within a few inches of the leading edge of the 
wing, where the wing is extremely thin. In order to 
simplify the calculations, the flow over the entire wing 
was considered turbulent, even though the methods of 
solutions used apply equally well when both laminar and 
turbulent boundary layers exist on the wing. 


Recovery Factor 


The temperature in the boundary layer near the wing 
surface has a maximum value of approximately the 
stagnation temperature but varies with distance from 
the surface depending on the rate of heat transfer. In 
the correlation of heat-transfer data, the adiabatic wall 
temperature is defined as the temperature that would 
exist at the surface if the surface were perfectly insulated 
rather than the temperature that actually exists at the 
surface. The recovery factor relates the adiabatic wall 
temperature to the free-stream temperature ahead of 
the airfoil and the stagnation temperature in the follow- 
ing manner: 


r= (Tax To)/(T; — To) (3) 

Recovery factors have been determined experi- 
mentally and appear to be more dependent on the type 
of boundary layer (laminar or turbulent) than on the 
free-stream speed. The values of recovery factors 
which have been reported range from 0.8 for laminar to 


0.95 for turbulent boundary layers. For all numerical 
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solutions presented in this report a turbulent boundary 
layer was assumed, and a recovery factor of 0.92 was 
used even though there is experimental evidence that 
the recovery factor increases with increasing Mach 
Number. One result of the use of a constant recovery 
factor is that the calculated boundary-layer tempera- 
ture is the same for both upper and lower surfaces of 
the wing regardless of the angle of attack. 


INVESTIGATION OF HEAT FLow 


The analytical determination of thermal stresses de- 
pends upon the knowledge of the internal temperatures 
of the wing under consideration. The determination of 
these temperatures is, in general, a three-dimensional, 
transient-state, heat-transfer problem whose solution is 
extremely complex. 

In order to obtain an engineering solution combining 
reasonable simplicity and accuracy, an investigation 
was undertaken to determine the importance of the 
three possible types of heat flow—chordwise, beamwise, 
and spanwise. The following specific problems were 
studied: 

(a) Chordwise and beamwise heat flow assuming no 
temperature gradients in the beamwise direction. This 
method is presented in references 12 and 13 and verified 
experimentally in reference 16. 

(b) Chordwise, beamwise, and spanwise heat flow 
assuming no temperature gradients in the beamwise 
direction. This method is presented in reference 4. 

(c) Beamwise heat flow only assuming symmetrical 
boundary conditions. This method is presented in 
reference 7 and verified experimentally in reference 14. 
flow only assuming unsym- 
This method is pre- 


(d) Beamwise heat 
metrical boundary conditions. 
sented in reference 5 and verified experimentally in 
reference 14. 

(e) Chordwise and beamwise heat flow assuming 
symmetrical beamwise boundary conditions. This 
method is presented in reference 8. 

(f) Chordwise and beamwise heat flow assuming 
unsymmetrical beamwise boundary conditions. This 
method is presented in reference 10. 

The results of these investigations showed that the 
chordwise variation in temperature was always greater 
than the variation in any other direction, indicating 
that the thermal stresses of greatest magnitude would 
be produced by chordwise temperature gradients. 

The investigation also showed that the chordwise 
variation in temperature could most easily be obtained 
with a high degree of accuracy from the method listed as 
(a) above. The simplifying assumptions made for this 
method reduce it to a one-dimensional numerical step- 
by-step process that is easily solved. All the above 
methods were developed in terms of a dimensionless 
temperature parameter 7 such that 
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average beamwise temperatures. half of a diamond-shaped wing. 

Chordwise Temperature Distribution Ni = 2-1/2; (where 2; =beamwise height of right 
= on : face of element 7 (9) 
The method adopted for the determination of the 

chordwise temperature distribution consists of dividing A, = cjw(Ax)?(1 + A;)/(R: AA) (10) 

the wing cross section into a number of elements m of ” 

8 D; = (hui + Ini) (Ax)?/ (Riz) (11) 


equal chordwise length, Ax, as shown in Fig. 1, and 
choosing a small time interval, A6, where the tempera- 
tures are obtained at the end of successive time inter- 
vals. The dimensionless temperatures of any element 
i are then given by the following equation: 


° ° ° 
Ti = QTj-1 + bit; —1 + dit; + CiTi+1 » 
+<2t1<sm (5) 


where 
aq = (A; — Lf; 
b = (i + Df, 
d; = (Ai —-3 —X — Dif l<i<m (6) 
e = ef; 


f, = 1/(A, —- 1+ + D) 


a = by = () 
d, = (A, —2- BoD hi ’ jo 
: = |] 

a" 1 (7) 

fy = 1/(Ai + BoD) 

a m = a, = YW 

dm = (A . ie pe fy 5 aS me . . 
; 1=>m (S) 

Cm = 0 


Son = l (A m + Xm + Bm Pin ) 


and 


Bo = ] a 20 2Ax (12) 
Bn a | + 2m 2Ax (13) 


The superscript ° with a temperature indicates the 
temperature at the end of the time interval immediately 
preceding the time interval under consideration. Thus 
all r°’s in Eq. (5) are known. Temperatures are deter- 
mined successively from left to right so that 7;_,, which 
refers to the element preceding element 7, will also be 
known for the determination of 7; in Eq. (5). 

Fig. 1 shows the wing cross-section layout for this 
method, as well as typical chordwise temperature dis- 
tributions. A sample calculation illustrating the use of 
this method is given in reference 12. 

It was found that this method yielded results of com- 
parable accuracy to those of the more complicated 
methods (b), (e), and (f) previously described, and, for 
the region near the leading and trailing edge, the results 
of this method were much more accurate than those of 
methods (c) and (d). A discussion of the limitations of 
beamwise solutions for the determination of chordwise 
temperatures is given in reference 5. 

Another advantage of this method of solution is that 
it is applicable in the case where heat-transfer co- 
efficients and adiabatic wall temperatures vary with 


time as well as chordwise station. Thus, it can be used 
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FLIGHT AND 


for any prescribed flight plan including take-off, ac- 
celeration, steady flight, deceleration, and landing. 

An experimental program was undertaken to verify 
the accuracy of this method for determining chordwise 
temperature distributions. A mild-steel, diamond- 
shaped wing model with thermocouples installed in- 
ternally along the chord and connected through a gal- 
vanometer to a recording oscillograph was suddenly 
immersed in a water bath that was at the boiling point. 
Transient temperatures were measured, and, based on 
calculated heat-transfer coefficients for the water bath, 
good agreement with the theoretical solution was ob- 
tained. The details of this program are given in ref- 
erence 16. 


INVESTIGATION OF THERMAL STRESSES 


The determination of thermal stresses requires that 
the internal and surface temperatures throughout the 
wing be known at any time for which the stresses are to 
be found. Due to the wing configuration and the man- 
ner in which the temperature gradients occur, the pri- 
mary thermal stresses act in the spanwise or y direction. 
The stresses over a cross section at any spanwise station 
may be thought of as composed of two parts: (a) a 
stress produced by the chordwise variation of the 
average beamwise temperature which will be referred to 
as the average stress; (b) a stress produced by a beam- 
wise temperature variation at the point in question. 
Since both of these stresses act in the spanwise direction, 
they may be added algebraically, yielding the re- 
sultant stress at any point. 


Stresses Due To Chordwise Temperature Variations 


Numerical methods of solution for the stresses due to 
the chordwise temperature variation are presented in 
references 3 and 6. The chordwise temperatures con- 
sidered are the values of the average beamwise tem- 
peratures at any location along the chord. The argu- 
ment for the solution of reference 6, which is based on 
the assumption of plane stress in the plane of the wing 
cross section shown in Fig. 1 resulting from chordwise 
temperature variation and is limited to symmetrical 
diamond-shaped wings with blunt leading and trailing 
edges, may be summarized as follows :* 

If the chordwise temperature variation is nonlinear 
and is not symmetrical about the mid-chord, in addition 
to thermal stresses, a deflection of the wing will occur 
along the span in the direction of the chord. This same 
deflection could result from a mean linear chordwise 
temperature variation that would cause no stress, with 


* For wings having cross sections other than a symmetrical 
diamond shape, the method of reference 3 should be used, since it 
is not limited to any particular cross-section configuration. Both 
the methods yield identical results for diamond-shaped wings. 
The method of reference 6 is presented here since the concept of 
the mean linear temperature variation is helpful in visualizing 
thermal stresses. 
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the amount of the deflection roughly proportional to the 
slope. The difference between the actual chordwise 
temperature and this mean linear temperature at a 
point causes thermal stress. 

In the case of temperature variations that are sym- 
metrical about the midchord, no deflection results and 
the slope of the mean linear temperature variation that 
would cause deflection but no stress is zero. Thus, the 
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temperature difference causing a stress at a point is 
simply the difference between the actual temperature 
and the average chordwise temperature. 

The resulting solution for stresses due to chordwise 
temperature variations is 


Oye = —ak(T — T.,) (positive when tensile) (14) 


where 
T = the average beamwise temperature at the 
chordwise station in question 
T,, = a+ bx = the mean linear chordwise tempera- 
ture variation that would cause deflection 
but no stress, evaluated at the chordwise 


station in question (15) 


The constants a and b above are evaluated as follows: 


[2mcl, + 4(m — n)(I2 — 13)]/[e2(m + n)] (16) 


c= 


b = [48mcJ, + 96(m — n)(Jo — Jz)]/[c4(m + 3n)] 


(17) 
where 
c = chord length 
m = ('/. maximum thickness of wing cross section) 


+ (chord length) 
n = ('/, leading-edge thickness of wing cross sec- 
tion) /(chord length) 


By we [ Tie (18) 
0 

I, = f . i xax (19) 

i / Tee (20) 
Jf 0 

i Tele wm ht he (21) 
0 

J, = J T x? dx (22) 


J3 = f T x? dx (23) 
J0 


The temperatures 7) may be determined from the 
solution for combined chordwise and beamwise heat 
flow outlined in the preceding section, and the integrals 
Eqs. (18)—(23) may be evaluated by Simpson's rule or 
some other method of mechanical integration. 

An experimental program was undertaken to verify 
the accuracy of this method for determining stresses in a 
solid diamond-shaped wing. A mild steel model wing 
with surface strain gages and thermocouples installed 
along the chord at the mid-span was used for the ex- 
periment. A nonsymmetrical chordwise temperature 
distribution was obtained by passing boiling water 
through a channel surrounding the leading edge, with 
cold water passing through channels at the mid-chord on 
the top and bottom of the wing and water at an inter- 
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mediate temperature passing through a channel suyr- 
rounding the trailing edge. Based on the measured 
temperature distribution and strain-gage readings cor- 
rected for temperature variation, excellent agreement 
between the stresses measured experimentally and 
those calculated from the theoretical solution was ob- 
tained. The details of this program, including the 
strain-gage calibration, are described in reference 15, 


Thermal Stresses at a Sharp Leading or Trailing Edge 


In the case of a wing having a sharp leading or trailing 
edge, the maximum stress at the leading or trailing 
edge may be easily calculated from Eq. (14), where 7,, 
is equal to the initial metal temperature. This is due to 
the fact that the leading or trailing edge would in- 
stantly attain the adiabatic wall temperature while 
the rest of the wing remained at the initial wing tem- 
perature. This would be the most severe thermal stress 
condition possible, and the effect of leading-edge thick- 
ness on stress will be discussed in a later section of the 


report. 


Stresses Due to Beamwise Tetnperature Variations 


A numerical method of solution for the stresses at any 
point due to the beamwise temperature variations is 
presented in reference 6. This solution, which is based 
on the assumption of plane strain in elements normal to 
the chord shown in Fig. 1, is 


oy» = —lak/(1 — v)|(T — T,,) 
(positive when tensile) (24) 
where 
T = the temperature at the point in question 
Tn = 8+ yz = the mean linear beamwise tempera- 


ture variation that would cause deflection 
but no stress, evaluated at the point in ques- 


tion (25) 


The constants 6 and y have the following values: 


1 1/2 
rt dz 
idm 
12 {'/ 
T 2 dz (27) 
[3 P 1/2 


/ = beamwise thickness at the point in question 


(26) 


where 


The temperatures 7 may be determined from the 
solution for unsymmetrical beamwise heat flow® and 
the integrals, Eqs. (30) and (31), may be evaluated 
by Simpson’s rule or some other method of mechanical 
integration. 

In the case when the beamwise temperatures are 
symmetrical about half the beamwise thickness, the 
temperature 7,, becomes the average beamwise tem- 
perature and may be evaluated analytically.’ 
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FLIGHT 


The resulting thermal stress due to both chordwise 
atid beamwise temperature variations is the algebraic 
sum of Eqs. (14) and (24)—that 1s, 

(28) 


CO, = Oy + Tyb 


Fig. 2 shows a typical stress distribution for the front 
half of a wing. The “average stress’ shown in Fig. 2 is 
the stress due to the chordwise temperature variation 
described in the preceding section. The stress distribu- 
tion for the rear half of the wing would be similar to 


that shown. 


Discussion 


The results of a number of calculations have shown 
that the maximum stress in the wing is a compressive 
stress occurring at the leading edge. At or near the 
leading edge, the beamwise temperature variation is so 
small as to be insignificant in the determination of 
thermal stress, indicating that Eq. (14) is sufficient for 

the maximum This is 
> 


shown graphically in Fig. 2. 


the determination of stress. 


In the thickest portion of the wing —that 1s, near the 
mid-chord——the beamwise temperature variation may 
cause appreciable stress but the magnitude of this stress 
is always much less than the stress at the leading edge. 
Further, the stresses due to the chordwise variation of 
the average beamwise temperatures result in a tensile 
stress in the region near the mid-chord, while the beam- 
wise temperature variation itself produces compressive 
stresses at or near the wing surfaces and tensile stresses 
in the region of the center of the beamwise thickness. 
The algebraic sum of these two stresses as given by Eq. 
(28) reduces the magnitude of the surface stresses while 
increasing the magnitude of the stresses near the center 
of the beamwise thickness. This effect is also shown in 
Fig. 2. Due to the nature of the temperature gradients, 
the compressive stresses are always greater than the 
tensile stresses, and it is felt that the determination of 
stresses based on the chordwise variation of the average 
beamwise temperatures, Eq. (14), is adequate for the 
overall of thermal stresses from the 
These stresses are shown as the 


determination 
standpoint of design. 
average stresses of Fig. 2. 

Due to the fact that the heat-transfer coefficient at 
the leading edge is greater than that at the trailing 
edge, for wings having finite leading- and trailing-edge 
thicknesses, the maximum transient stress will occur at 
the leading edge. At some later time, however, the 
trailing-edge stress may be greater than the leading- 
edge stress at the same time. The effect of sharp leading 
or trailing edges has been previously discussed. 


INVESTIGATION OF THE EFFECT OF FLIGHT CONDITIONS 
ON THERMAL STRESSES 


The methods of calculating temperature distributions 
and thermal stresses previously discussed were used in 
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an investigation to determine the effect of variation in 
important flight conditions on thermal stresses. 

The stress magnitudes presented in the following 
sections are extremely high, and, in general, buckling of 
the leading edge would occur before the stresses were 
attained. It should be remembered, however, that the 
stress conditions investigated resulted from an assumed 
infinite acceleration to the given flight condition and 
thus represent the most severe thermal stress condition 
possible. 

The stress magnitudes are significant in that they 
represent the values that the stresses would tend to 
It is be- 
in stress presented are 


approach if the wing material did not buckle. 
lieved that the 
representative in showing the general trends caused by 


variations 


variation in the various flight and configuration con- 


ditions. 


The Effect of Angle of Attack 


To determine the effect of angle of attack on thermal 
stresses, the wing was assumed to be at a Mach Number 
of 4 at 60,000 ft. standard altitude having an initial 
uniform temperature of 393°R. The chordwise varia- 
tion in heat-transfer coefficients for the upper and lower 
wing surfaces as calculated by Eq. (2) for angles of 
attack of 0°, 2°, 4°, and 6° are shown in Fig. 3. The 
chordwise temperature distributions for the various 
angles of attack were determined over a sufficient time 
period for the maximum stress to occur. These calcula- 
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tions showed the maximum stresses to be unaffected by 
a variation in angle of attack. Calculations also showed 
that the average beamwise stress at any chordwise 
station was independent of angle of attack. 

In Fig. 3 the heat-transfer coefficient is seen to vary 
by a factor of approximately 2 for an angle of attack 
variation of 0°-6°. The reason this variation in heat- 
transfer coefficient has no appreciable effect on the 
maximum stress at the leading edge may be explained 
briefly as follows: 

For wings with finite leading-edge thickness, with the 
passage of time the leading-edge stress will first increase, 
pass through a maximum, and later decrease as thermal 


equilibrium is approached. The greater the heat- 
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transfer coefficient, the sooner will the maximum stress 
occur. However, in the case of a sharp leading edge, 
the heat-transfer coefficient does not enter into the 
equation for maximum stress in any way. Thus for qa 
leading edge that is extremely thin, it is not surprising 
that the heat-transfer coefficient has no appreciable 
effect on the value of the maximum stress, although it 
does affect the time at which the maximum stress occurs, 
Further effects of time on maximum stress are discussed 
in a later section. 


The Effect of Mach Number and Velocity 


To determine the effect of Mach Number on thermal 
stresses, the wing was assumed to be at zero angle of 
attack having an initial uniform temperature of 393°R. 
at 60,000 ft. altitude. The chordwise variations in heat- 
transfer coefficient for this condition are shown in Fig. 4 
for Mach Numbers of 2.5, 4, and 6, corresponding to 
velocities of 1,658, 2,650, and 3,970 m.p.-h. 

The chordwise temperature distributions for the 
various Mach Numbers were determined over a suf- 
ficient time period for the maximum stresses to occur. 
The resulting variation in maximum stresses occurring 
at the leading edge plotted versus Mach Number and 
velocity are shown in Fig. 5. 

As shown in Fig. 5, the maximum stress increases 
rapidly with Mach Number and velocity, being 4.78 


5 


times as great at M = 6asat M = 2.5. 


The Effect of Altitude 


To determine the effect of altitude on thermal 
stresses, the wing was assumed to be at zero angle of 
attack having an initial temperature of 530°R. and 
suddenly subjected to a Mach Number of 4. In order to 
simplify the calculations, an average value of heat- 
transfer coefficient as calculated from Eq. (2) was used 
at each altitude. The chordwise temperature distribu- 
tions for altitudes of 0, 20,000, 40,000, 60,000, and 100,- 
000 ft. were calculated over a sufficient time for the 
maximum stress to occur. The resulting variation in 
maximum stress occurring at the leading edge plotted 
versus altitude is shown in Fig. 6. 

As shown in Fig. 6, the maximum stress decreases 
somewhat with an increase in altitude, being 1.47 times 
as great at sea level as at 100,000 ft., again indicating 
that the magnitude of the heat-transfer coefficient has 
little effect on the maximum thermal stress. The mag- 
nitude of the heat-transfer coefficient does have a great 
effect on the time at which the maximum thermal stress 
occurs, as will be discussed in the following section. 


The Effect of Time 


Since all problems investigated in this study are 
based on transient-state heat transfer, time is an im- 
portant variable. In general, the effect of time on the 
stresses at any point in the wing is that the stress in- 
creases with time to some maximum value and then de- 
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FLIGHT AND 


creases approaching zero at infinite time. A typical 
variation of stress with time is shown in Fig. 7. The 
conditions chosen for this figure are for a Mach Number 
of 4 at 60,000 ft. at 4° angle of attack. The stresses 
plotted are the leading-edge stresses. 

Variation in the flight conditions or wing configura- 
tion also affects the time at which the maximum stress 
occurs. As an example, Fig. 8 shows the time required 
to attain the maximum stress for various altitudes based 
on the conditions of the preceding section. As shown in 
Fig. 8, the maximum stress at sea level occurred 46.7 
sec. after the wing was subjected to the flight Mach 
Number, while the time required at 100,000 ft. was 
1,480.1 sec. This effect is due primarily to the large re- 
duction in heat-transfer coefficient with increasing alti- 


tude. 


INVESTIGATION OF THE EFFECT OF CONFIGURATION 
CONDITIONS ON THERMAL STRESSES 


The methods of calculating temperature distributions 
and thermal stresses previously discussed were used in 
an investigation to determine the effect of variation in 
important wing configuration conditions on thermal 
stresses. The purpose of this investigation was to com- 
pare the effects of certain conditions and not to es- 
tablish accurate stress magnitudes. For this reason, an 
average value of heat-transfer coefficient was used for 
all calculations except those for the effects of solidity on 


stress. 


The Effect of Leading-Edge Thickness 


To determine the effect of leading-edge thickness, the 
airfoil was assumed to be a Mach Number of 4 at 20,000 
ft. with 530°R. initial temperature. The temperatures 
were calculated over a time range sufficient to determine 
the maximum stress for an average value of heat-trans- 
fer coefficient based on Eq. (2). 

The leading-edge thickness was considered as a 
variable, and maximum stresses were determined as a 
function of leading-edge thickness for values of leading- 
edge thickness from 0—0.1 ft., corresponding to leading- 
edge-to-maximum thickness ratios of 0-50 per cent. 
The results of this investigation are shown in Fig. 9. 

As shown by Fig. 9, the maximum stress is greatly 
affected by leading-edge thickness. With a leading 
edge 20 per cent of the maximum thickness, the maxi- 
mum stress is only 38 per cent of that for a sharp leading 
edge. With a leading-edge thickness 5 per cent of the 
maximum thickness, the maximum stress is 68 per cent 
of that for a sharp leading edge. 

It should be remembered that this curve is based on a 
specific wing having a 4-ft. chord and 0.2-ft. maximum 
thickness. The results should be considered to repre- 
sent a trend rather than as a guide to absolute values. 
Subsequent investigations showed that the leading-edge 
thickness has a much greater effect on stress than do the 
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Fic. 9. Maximum stress versus leading-edge thickness 


chord length, thickness ratio, or span, as will be dis- 
cussed in the following sections. 


The Effect of Chord Length 


To determine the effect of chord length, the airfoil 
was assumed to be at a Mach Number of 4 at 20,000 ft. 
altitude with 530°R. initial temperature. The tem- 
peratures were calculated over a sufficient time range to 
determine the maximum stress with an average value of 
heat-transfer coefficient based on Eq. (2). 

The chord length was considered as a variable with a 
constant leading-edge thickness-to-maximum thickness 
ratio and maximum thickness-to-chord ratio. That is, 
cross sections remained geometrically similar. Maxi- 
mum stresses were determined as a function of chord 
length for values of chord length from 3 to 7 ft. The re- 
sults of this investigation showed the maximum stress to 
be unaffected by a variation in chord length based on the 


above conditions. The time at which the maximum 
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Time required to attain maximum stress as a function 
of thickness-to-chord ratio. 


Fic. 11. 


stress occurs, however, increases with increasing chord 
length as shown in Fig. 10. 


The Effect of Thickness-Chord Ratio 


To determine the effect of the thickness-to-chord ratio, 
the airfoil was assumed to be at a Mach Number of 4 at 
20,000 ft. altitude with 530°R. initial temperature. 
The temperatures were calculated over a sufficient time 
range to determine the maximum stress with an average 
value of heat-transfer coefficient based on Eq. (2). 

The maximum ratio 
sidered as a variable with a constant leading-edge 
thickness-to-maximum thickness ratio and chord length. 


thickness-to-chord was con- 


Maximum stresses were determined as a function of 


thickness-to-chord ratio for values of 1, 4, 5, 6, and 8 
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per cent. The results of this investigation showed the 
maximum stress to be unaffected by a variation jp 
thickness-to-chord ratio in the range of values invest- 
gated. As in the case of variable chord length, the time 
at which the maximum stress occurs increases with in- 
creasing thickness-to-chord ratio as shown in Fig. 11, 


The Effect of Span and Taper 


An investigation of combined chordwise, beamwise, 
and spanwise heat flow given in reference 4 showed that 
spanwise heat flow had a negligible effect on the chord- 
wise variation of the average beamwise temperature, 
even in the case of a wing having a two-to-one taper, 
Later, calculations for a wing with three-to-one taper 
also showed that spanwise heat flow had a negligible 
effect on the chordwise variation of the average beam- 
wise temperature. 

As a result, it was found that the maximum thermal 
stress was unaffected by span or taper for wings in the 
For the 
tapered wings, the time required to attain maximum 
stress was less in the region near the tip than that for the 
It should be remembered, how- 


taper range of one-to-one to three-to-one. 


region near the root. 
ever, that the method for solution of thermal stresses 
was limited to the region of the wing one chord length or 
more from the wing tip, since St. Venant’s principle was 
apphed in the development of the stress calculation 
methods. 


The Effect of Sweep 


All the solutions presented for stress are based on the 
simplifying assumption that the wing is unrestrained. 
Therefore, by choosing the spanwise or y direction 
parallel to the mean angle of sweep and the chordwise or 
x direction normal to the spanwise direction, the 
methods of solution may be applied directly to swept- 
back wings. 

It was felt unnecessary, however, to perform calcula- 
tions to determine the effect of sweep in view of the re- 
sults of the preceding sections on the variation of chord 
The effect 
of sweep is to change the chord length, thickness-to- 


length, thickness-to-chord ratio, and span. 


chord ratio, and span slightly, and small changes in 
these parameters have been shown to have negligible 
effects on maximum thermal stresses. 


The Effect of Solidity 


To determine the effect of solidity, the basic wing was 
assumed to be at a Mach Number of + at 60,000 ft. alti- 
tude with zero angle of attack. The solidity was varied 
by considering wings having sections removed from the 
center of the cross section which were geometrically 
similar to the wing cross section as shown in Fig. 12. 
A solidity ratio was defined as the maximum thickness 
of the cross section minus the maximum thickness of the 
section removed, all divided by the maximum thickness 


of the cross section. That is 
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FLIGHT AND 
Solidity ratio = (¢ — f’)/t (29) 
where 
tf = the maximum thickness of the cross section 
t/ = the maximum thickness of the section removed 


The solidity ratio was varied from unity (solid wing) 
down to 0.125. The temperatures were determined 
based on the assumption that the air inside the wing 
would instantaneously attain the inner surface tempera- 
ture at any chordwise station. A sufficient time range 
was covered to obtain the maximum stress at the leading 
edge. Thermal stresses were calculated by the method 
of reference 3, since the method of reference 6 pre- 
viously presented applies only to a solid, diamond- 
shaped wing. Fig. 12 shows the variation in maximum 
stress with increasing solidity ratio. 

As shown in Fig. 12, the maximum stress is greatly 
reduced for wings of low solidity ratio, being only 48 per 
cent as great for a solidity ratio of 0.125 as for a solid 
wing. Calculations also showed that the maximum 
tensile stress in the region of the mid-chord was appre- 
ciably reduced with decreasing solidity ratio. 

The results of this investigation indicate that a 
hollow wing of low solidity would be greatly superior to 
a solid wing from the standpoint of thermal stresses. 


CONCLUSIONS 


(1) The maximum thermal stress, which is com- 
pressive, occurs at the leading or trailing edge of the 
wing and can be determined from a one-dimensional 
solution for the chordwise variation of the average 
beamwise temperature. 

(2) For the limiting case of infinite acceleration to a 
steady flight Mach Number greater than two, the lead- 
ing-edge stresses for the chosen configuration would 
be prohibitive from the standpoint of exceeding the 
critical buckling stress. 

(3) The tensile stresses of greatest magnitude occur 
in the region of the wing near the beamwise center of 
the mid-chord. These stresses are much smaller than 
the compressive stresses at the leading or trailing edge. 

(4) Thermal stresses increase rapidly with increasing 
Mach Number (see Fig. 5). 

(5) Thermal stresses may be reduced appreciably by 
increasing the leading-edge (and trailing-edge) thickness 
(see Fig. 7). 

(6) A hollow wing of low solidity is greatly superior 
to a solid wing from the standpoint of thermal stress 
(see Fig. 12). 

(7) Thermal stresses decrease somewhat with in- 
creasing altitude (see Fig. 6). 

(S) The following variables have no appreciable effect 
on maximum thermal stresses: angle of attack, chord 
length, thickness-to-chord ratio, length of span, taper, 


and sweep. 
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Effect of Turbulence Level of 


Incident Air Stream on Local 


Heat Transfer and Skin Friction on a Cvlinder 


(Continued from page 730) 


physical point of view, so the curve of Fig. 10 was used 
instead. This type of a plot brings out the fact that the 
pressure reading of a surface tube is due to combined 
viscous and inertia forces, the relative contributions of 
each depending on the fluid flowing and on the physical 
dimensions of the tube. 

The points on Fig. 10 used to establish the calibration 
curve were determined from pressure measurements 


made in the laminar boundary-layer region and using 
the Blasius-Hiemenz boundary-layer solution. Agree- 
ment between the skin friction in the laminar boundary- 
layer region calculated from this curve and that deter- 
mined from the Blasius-Hiemenz method was within 12 
per cent for most of the data taken, and it is felt that this 
is about the best accuracy that should be claimed for the 
skin-friction curves in general. 





Shear Center for Thin-Walled Open Sections Beyond 
the Elastic Limit 


(Continued from page 754) 


CONCLUSION 


The methods outlined above can be used to find 
the position of the shear center for thin-walled sec 
tions. It should be emphasized that the shear center 
is no longer a section property but depends for its loca- 
tion on the load and the stress-strain curve. Further 
more, while it has been shown that in the elastic case® a 
load acting on such a section can be replaced by a shear 
load at the shear center and a torque with reasonable 
accuracy, it still remains to be studied whether this is 
true beyond the elastic limit. Consequently, it is still an 
open question whether the results given in this paper 
can be used for such a purpose. 
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A Simple Graphical Method for Constructing 
Families of Nyquist Diagrams* 


Robert M. Stewartt 
Jet Propulsion Laboratory, California Institute of Technology 


July 20, 1951 


INTRODUCTION 


(1) 


I REFERENCE 1, an outline is given of the widely used fre- 
quency-response method of analysis of aircraft-autopilot 
stability. By means of this method, it is possible to determine 
the stability of an automatically controlled aircraft from meas- 
ured autopilot characteristics and calculated aircraft character- 
The criterion used in interpreting frequency-response 
It is sometimes 


istics 
data is the Nyquist criterion or variations of it. 
necessary to make a comprehensive survey of aircraft stability 
over widely varying conditions of operation, with possible vari- 
ations in aircraft configuration and with various combinations 
of control-system parameters. It is the purpose of this note to 
describe a simple graphical aid to constructing such families of 
An example is given concerning a system whose con- 
It is as- 


diagrams. 
trol-system response is a function of two parameters. 
sumed that these parameters correspond to physical character- 
istics of the equipment which are accessible for change and ad- 
justment. It is then shown that, with other conditions remain- 
ing unchanged, Nyquist diagrams corresponding to any com- 
bination of values of these two parameters can be constructed 
using only a straightedge and dividers after four Nyquist dia- 
grams have been determined in the usual manner correspond 


ing to four combinations of the two parameters. 


(II) Form or Basic STABILITY EQUATIONS 


Following reference 1, let 


0 = (0/8)a(p)-5 (1 


define the calculated aircraft transfer function and 


56 = (6/0)c(p)-0 (2) 


define the measured control-system transfer function. The char 
acteristic equation of the system, whose roots indicate the sta- 
bility and transient response of the system, is given by the simul- 
taneous solution of Eqs. (1) and (2) 

(0/8)a(p)-(8/0)c(p) = 1 (3) 
Hence, stability is determined by the relationship of the contour 


G(jw) = —(0/5)a(jw)-(6/A)c(Fw) (4 
* This paper presents the results of one phase of research carried out at 
the Jet Propulsion Laboratory, California Institute of Technology, under 
Contract DA-94-495-ORD 18, sponsored by the Department of the Army 
Ordnance Corps 
Tt Research Engineer, Systems 
Section 


Analysis Group, Guidance and Control 


J 
N 


to the point (—1, 70). Many common types of control systems 
are so arranged that it is known that their transfer function is of 


the form 


(8/9)c(jw) = [go(jw) + Kigi(jw) + Koge(jw) + Ki Kogieljw 


(5) 
where K, and Ke are pertinent control-system parameters. Thus, 
the function that is mapped is of the form 
G(jw) = —(6 5)a(jw) - [gol jw) + Kigi(jw) + 

K2gA jw) + Ki Kogiel jw (6 
(III) Properties or G(K,, Ke) 
Now consider G as a function of K,; and Ky. If w and Ky, are 
constants and if we let 
x = [G(Jw) |reat / 
a ( 
J = [G( J) |imaginary ' ‘ 
then 
x =a, + b,K,) 
8) 
y = a2 + boKi\ 
where a, a2, >;, and b; are real constants, or 
y=ct+dx (9) 


where c and dare real constants. Thus, if K is a constant and if 
we vary K,, lines of constant w are straight lines. Furthermore, 


from Eq. (8), 


Ax/AK, = hi, Ay/ AK, = by (10) 
J + IMAGINARY a | a ae 
| 
. + REA 
Fic. 1 G( jw) for Kz = 20, 
i. 
- 


Fic. 2. G( jw) for Ke = 40 
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+IMAGINARY 
- TRACED FROM Figs. | & 2 
K, ie -- CONSTRUCTED 
| | 
a - - , ~ : —— _+ REAL | 
° ' 
? 
c 
/ P. 
x 
: x . 
| : 
ai 
aay a | 
Fic. 3. G( jw) for K; = 16. 
or if 
As = (Ax? + Ay?)'/? (11) 
then 
As = (bi? + bo?)'/2- AK, (12) 
or 
As = (constant): AK, (13) 


which means that along these lines of constant w the displacement 
of the locus of jw is proportional to the change in gain. Similar 
results are obtained if we let A, be constant and vary Ko. 

It is now possible by simple construction to obtain the stability 
diagram for any combination of A, and A» after we have com 
puted and plotted four loci. This construction is illustrated in 
Figs. 1, 2, and 3. 

In Fig. 1, we have mapped the locus ot jw for 


K». = 20, Kk, = 8 
K.=20, K,=12 


From these two curves we can readily obtain the plot correspond- 
20 and any desired value of K,. The plot corre 
16 is shown as an illustration. In Fig. 2, we 


ing to Ky = 
sponding to K,; = 
have mapped the locus of jw for 


K, = 40, K,= 8 
K,=40, K, = 12 


From these two curves we can readily obtain the plot correspond- 
ing to Ke = 40 and any desired value of Kj. 

Suppose that the plot for K, = 30 and K, = 16 is desired. 
First, the plot for K. = 40 and K, = 16 was constructed as shown 
in Fig. 2. Then the two plots (from Figs. 1 and 2) correspond- 
ing to K, = 16 were traced on Fig. 3. Now by the same pro- 
cedure as before, we obtain the plot for K, = 16 and any desired 
value of Ky. The construction of the curve corresponding to 
K, = l6 and Kz = 30 is shown. 

(1V) CoNcLuUSION 

In general, if it is desired to apply this procedure to the effect 
of » variable parameters, it is necessary first to compute and 
map 2" fundamental contours. The procedure is obviously ap- 
plicable to electronic and other systems that are susceptible to 
analysis by means of the Nyquist criterion and to the analysis 
of the effect on stability in aircraft of inertial and aerodynamic 
parameters. Some manipulation is required in certain cases in 
order to obtain a function that will determine stability and at 
the same time be of the form of Eq. (6) with respect to the lumped 


parameters that are to be varied. Typical cases of aircraft con 
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trol to which the method is particularly adaptable are (a) co; 
bination rate and position gyro control and (b) combination 


ground station and internal control. 


REFERENCE 
' Seacord, Charles L., Application of Frequency-Response Analysis to Air. 
craft-Autopilot Stability, Journal of the Aeronautical Sciences, Vol. 17, No 
8, pp. 481-498, August, 1950 


Strength of Reflected Shock in Mach 
Reflection* 


H. S. Tan 
Cornell University 
July 13, 1951 


wu A PROPAGATING PLANE SHOCK Is diffracted by a corne: 
Mach reflection is set up, the situation after passage of the 
wave being as sketched in Fig. 1. The reflected front R is an 


acoustic wave, according to first-order theory.'~* However, 
results of first-order theory can be used to determine more accu 
rately the strength of this wave, and it will be shown here that 
this is in fact a shock of second-order strength. The case of ex 
tremely weak incident shock was studied by Bargmann,! where 
vorticity in region O7'P was neglected. The present investiga- 
tion, however, is not limited in incident shock strength. Light 
hill’s first-order pressure solution® and his technique‘ for finding 
the shock strength have been used. The latter 
process of pushing away the apparent singularity appearing in 
the differential equation of motion due to its linearization 


Since for small 3 the flow behind straight shock DT and slip 


is essentially a 


stream 5S is irrotational [vorticity 2 < O(#*)]|, we can introduce 
a velocity potential @ which, because of its couical property in 
time-space, further reduces to the form 


’(7, 0, 1) = a,*tf(r, 0) l 


in a coordinate system fixed to flow 1, with origin at the disturb 
ance center O, and a; = yPi/pi,7 r/ayt. 
Denoting partial derivatives by subscripts and the radial and 


tangential velocity components by u, v we have 


u = % = af, 
up = —(ayr/t)frr 
v = (1 na de = (da r) fo 
vy = —(ayr/t) [—(fo/r?) + (1/r)fol 
and 
?, = af — rf-) 
Vb = (1/t)V,*f 
O,, = (a,?7?/t) fy 


On introducing these relations into the equations of motion, 


1 a? a,? 
>, + (u? + v?) + = 9 
2 y-1 y-1!1 
a Ou; Ou, = 
av 7b = Oy + 2u; UM; 3 
ot Ox; 


and eliminating a/a;, the following result is obtained: 


11 
1-e--0(r4 it) + (0 
2 r? 


1 
(1 - 14) frr + (additional terms in f,g and fag) = 0 (4) 


~ i) 


* This investigation has been part of a research project sponsored by the 
U.S. Air Force, Office of Air Research, at the Graduate School of Aeronau 
tical Engineering, Cornell University A more detailed account has been 
submitted to that Office for publication as an Air Force report. 
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where additional terms, not explicitly required here, have not 
been written out. 
Referring to coordinates fixed to the shock, the Rankine 


Hugoniot conditions can be written 
(a q" continuous 
(b) Ag’” = [2/(y + 1)] [(a2/g’”) — gq if 
where q" and i denote velocity components tangential and 
normal to the shock, respectively, and A( ) denotes the increase 
of a quantity across the shock front. Now, since @ exists and the 
flow is stationary in the (7, @) plane, let the reflected shock front 


ber = ni @): then 
t fa + (0) (f, — n) 
q ay 
rV1l+pr-2n”? 
- —( lo ~ @) 7 9 2n'fa 
q = ay 
Vic+r- yn” 
Alfa + n’(@) Uf n)| = 0 
So 
sie Af, — 7 29’? Af, : , 
Ag ay —Vi1 + 7-44", of, 


df = ——~ vy — I) U — of) — 2f, + Un — 1)}. 
4f=0 (6) 


Following Lighthill,* we imagine both 7 and f expanded in series 
of terms of ascending order: 
r R+/7(0) = R+7,(0) + r-(@) 4+ _t = 
e (4) 
f = fi(R, 0) + fo(R, 0) 4 { 


Then f; is the familiar linearized result. The process of suc 
cessive calculation of 7, 72, should be thought of as a process 
of successive pushing away of the apparent singularity that ap 
pears at R = 1 in the approximate equations. Lighthill has 
pointed out that, physically, the process amounts to applying the 
results of linearized theory at an improved location in the field 


Since f} = fir = Oon R = 1, fir ~ A(O) V1 — Ras R—~1in 


linearized theory, and 





a 1 >, —/f-—+ rf 
p p = = _ ’ = 1(0)V1—R, 
p, U7? M;? 
R—1 (8) 
So 
fire = M7A(0) V1—R (9) 
O- 
1 D 
Rn &% 
i 
~ 
se i 
. . 
2 - 
ri 
a.t SH 
o 
rd 
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Now if Eq. (7) is put into Eq. (4), and we collect terms of the 
second order, we obtain a differential equation for f.(R, @). 
Here a term in /irr appears, and this derivative is infinite at R = 
] To avoid infinite fxr, we must choose r,(@) to make the co 
efficient of fier vanish at R = 1. The result is 7;(@) 0; Le., no 
first-order difference between r and R is required in this case 
Evaluating the limiting value of this fizgr term as R — 1, we can 
use the differential equation to calculate (fsr)r=1; the result is 


(for)r=1 = [Cy + 1)/2] 144A) (10) 


We proceed to find 7.(@) by a similar process to that used for 
r\(9); viz., to make the coefficient of f gr in the differential equa- 
tion for f;(R, 0) vanish at R = 1. We tind 


r(@) = [(y + 1)?/4] AN4A (8) (11) 


The shock wave is now to be fitted into the region of validity 
of our results, which we have established to second order. The 
stock conditions [Eq. (6)] at r = (@) are: Since f = 0 ahead of 
shock, so directly behind it 

f=0 fr = 4n — 1)/(y + 1) (12) 


The shock comes in the range R < 1, and » — 1 (= m, say) is 
of second order. The first condition |[Eq. (12)] is satisfied by Tay 
lor expansion about R = 1 to the second order and tells us only 
that f. = Oat R = 1 fas has already been assumed in getting 
Eq. (11)]. The second condition [Eq. (12)] involves two second 
order terms; namely, (for)r=1 and the contribution of fiz, which 
is ,2A(0) V rs — m. Thus, the second shock condition leads to 

m = (3/16) (y + 1)?14,14°%(0) if A(@) > O 
= () if A(9) < 0 


For the case of diffraction by a concave corner, as in Fig. 1, 
(0) > O, and the shock is located at 
r= 1 + (3/16) (y + 1)710,14°(@) (14 


Now the pressure change across the shock is 


pl? 3 ; . 
pe, = Me (rf, — f) = , vy + 1)441A1%0) 
So the shock strength is given by 
Ap 4yU? 3 
P = : (0) = yy + 1)1%,'!A%0) (15) 
pb (y+ DM a? { 
where 
A(@) = lim [((p — pri)/ml?% LF (16 
r—>l 


on linear theory. 
This limit as given by reference 3 is of the form 


A(@) = (function of 0) X 3 = A(@)d, say (17) 


There is special interest in the strength of the reflected shock at 
the triple point 7, especially in view of the well-known difficulties 
of interpreting certain experimental resultsin Mach reflection. Ac- 
cording to reference 3, K(@) in Eq. (17) is zeroat T, so the strength 
of the reflected shock vanishes there. We conclude, therefore, 
that while this shock is O(#?) over the rest of its front, it can be 
at most O(3*) at T 

It is also of interest to determine the angle at which the re 
flected shock meets the main shock at 7. The question is 
whether this angle differs, in the second order, from the angle 
formed by the main shock and the acoustic circle with center at 
0 (Fig. 1). This can be determined from Eq. (14) above, be- 
cause the second-order correction to this angle is given by the 
value of dr/rd@ at T. This quantity has been worked out using 
Lighthill’s A(@) and has been found equal to zero. Thus, the 
wave is not only of vanishing strength at the triple point, but it 
is locally tangent to the acoustic circle there. 
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A Simple Rule for the Velocity Rise with 
Subsonic Mach Number on Ellipsoids of 
Revolution 


D. Kuchemann, 
Royal Aircraft Establishment, Farnborough, Hants, England 
July 17, 1951 


— SUBSONIC COMPRESSIBLE FLOW over a slender body of 
revolution can be determined to a first approximation from 
the Prandtl-Glaue1t method in the form of the ‘‘streamline anal- 
ogy”’ by Busemann! and Géthert.2 The velocity increment 
vz, in the direction of the main stream, in compressible flow is 
1/8? times the velocity increment v,, in incompressible flow on 
an analogous body, obtained by reducing the lateral dimensions 
B= V1-— M,?, and Mois 


The problem thus becomes 


of the original body in the ratio @: 1. 
the free-stream Mach Number. 
that of finding a relation between the velocity increment v,, of 
the analogous body and the velocity increment v,,; of the given 
bedy, both in incompressible flow 

From the results of exact solutions by various authors* ¢ for 
the maximum velocity increment on ellipsoids of revolution in 
incompressible flow we find that the relation 


Uz; = 0.6458°/* 


is a reasonable approximation up to a fineness ratio 6 = D/L of 


about 0.3. This is shown in Fig. 1. Hence, 


= ().645(85)*/? = p/%y,, 


ra 


so that 


Oe oad 2\e, = 3/2 / a2), 
v, = (1/8?)v,, = (6°/*/B?)u,; 


I 
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and finally, 


Vz aa ] l 
Vz; Cri V ‘B v I Tg M,? 


assuming the usual linearized relation between the pressure eo. 
efficients Cp = (p — po)/go and the corresponding velocity in 
crements. It will be found that this relation agrees well with ex. 
perimental results. 

v, or C, can be used to find Merit., the value of the free-stream 
Mach Number at which sonic speed is reached locally on the body 
A scale for Merit. obtained in this way is added to the right of the 
figure, which thus provides a simple means of estimating the 
greatest free-stream Mach Number at which ellipsoids of revoly- 
tion of known fineness ratio will be free from compressibility 


disturbances. 
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Ebene und rdéumliche Strémung bei hohen Unterschallge- 


Rotationskorpern bet axtaler Anstré- 
lranslation, Brit 


On Slender Wing Theory 


John W. Miles* 

Fulbright Lecturer, Auckland University College, Auckland, New 
Zealand 

August 29, 1951 


oe. HAS RECENTLY PUBLISHED an analysis of the 
supersonic flow over a low aspect ratio rectangular wing. His 
results (a) prove that the slender-wing theory, first applied by 
Jones? to pointed wings, does give the correct lift coefficient for 
the rectangular wing, albeit not pointed, in the limit of zero as- 
pect ratio and (b) furnish the pressure distribution aft of the re- 
gions where the results of Gunn* and Lagerstrom‘ are appli- 
cable. However, he does not calculate the lift and moment co- 
efficients for other than zero aspect ratio. 

An alternative approach to the boundary value problem is 
to separate the Laplace transform of the potential equation in 
elliptic coordinates and exhibit the solution as an expansion in 
Mathieu functions, following the classical analysis of diffrac 
tion around a flat strip.» With the exception of minor numerical 
discrepancies, the results are in agreement with those obtained 
by Stewartson, and we present here only the expressions for the 
lift and moment coefficients. 

Let A.R. be the geometrical aspect ratio, A the effective as- 
pect ratio, viz., 

1 = (iM? — 1 A.R l 
and Cy and Cy the lift and moment coefficients (with respect to 


the leading edge and with the chord as a reference moment 


arm); then 


CL. = 7A.R./2 i1 + (1/16)A2 — (3/64) [In(8/A) — (29/24)] X 
A4—1.17 exp (—1.364/4A ) cos [(2.22 8/A) + 0.275]} 2 
Cu = 7A.R./2 (—(1/8)A2 + (1/16) [In(8/A) — (9/8) ]A4 4 
exp (—1.364/A ) (1.17 cos [(2.22 8/A) + 0.275]+ 
0.45A cos [(2.228/A) + 1.296]}) (3 


* On leave from University of California, Los Angeles 
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The expansions of which Eqs. (2) and (3) give the leading terms 
are semiconvergent and represent the true results only in an 
asymptotic sense. In addition to terms of order A‘ln?A, ex- 
ponential terms are neglected,' since only one pole of the Laplace 
transform has been determined accurately. 

The approximation given by Eq. (2) agrees with the exact 
result {Cr = 4(M? — 1)-'/2 [1 — (2A)~"]} to better than 0.1 per 
centat A = land is indistinguishable from Lagerstrom’s result at 
A =1/2. The agreement of Eq. (3) with the exact results is, of 
course, less satisfactory but is within 3 per cent at A = 1. 

The foregoing results are to be compared with those recently 
obtained by Lomax and Sluder.6 Comparison of Eqs. (2) and 
(3) with the exact results of reterence 4 in 1/2 < A < 1 estab- 
lishes Stewartson’s analysis as the more accurate. 
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Proper Use of the M.1.T. Tables for Supersonic 
Flow Past Inclined Cones 


Antonio Ferri* 

Langley Aeronautical Laboratory, N.A.C.A., Langley Air Force 
Base, Va. 

July 27, 1951 


y THE May, 1951, issue of the JOURNAL OF THE AERONAUTICAL 
SCIENCES, a letter with the same title was published from 
Van Dyke, Young, and Siska of the RAND Corporation, in 
which expressions for the values of velocity components, density, 
and pressure at the surface of a conical, circular body were given 
as function of the values given by the M.I.T. tables. 

The expressions given are not valid for the components of 
velocity and for the density at the surface of the cone but are 
values outside a vorticose layer existing near the body.! 

The entropy at the surface of the body must be constant be- 
cause the flow at the body follows a streamline, and, therefore, 
Eqs. (le) and (2e) are not physical at the surface of the body. 
Near the surface of the body, a vorticose layer exists of zero thick- 
ness in the first-order approximation and of a thickness of the 
order of e when terms in e? are retained. Across the layer a vari- 
ation of entropy occurs from a value 


S=S,+ «>, Sn cos nd + e>_ Sn cos nN 


toaconstant value 


S' = 5S. +5, + 2s; 

The entropy changes the velocity components and density at 
the surface, and the expressions (2a), (2b), (2c), and (2e) are 
valid only outside the layer. However, the layer in the first ap- 
proximation has zero thickness and, therefore, does not produce 
changes in pressure. When terms of the order of ¢? are retained, 


expressions of the form 
p=prt «>>, cos nd 4 ep, cos nd 


p=pc >) —, cos n@ + e) Pp, COS Nd 
* At the present time, Professor of Aerodynamics, Polytechnic Institute 


of Brooklyn 
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would produce an entropy distribution at the surface of the form 
S = Sy, + €S, cos @ + ey S, cos n@ 


which again is not physical. 

The expression of the pressure Eq. (2d) is valid at the surface 
of the body and is consistent with the approximation accepted 
for «— 0; however, it is not sufficient when small but finite angles 
In this case, the pressure must be ob- 
Indeed, 


of attack are considered. 
tained from the complete expression of the velocity.! 


V2 = uo? + 2ujmoe cos @ + €2w,? + u,%e2 cos? @ + 


other terms in & 


Now, in Eq. (2d), the terms e?w,? sin? ¢ and u;*e? cos? @ are con- 

sidered small with respect to 2u,%¢ cos ¢, which is the only term 

retained. However, at the surface of the body, 
ew? = —e*u,?/sin? 8, 


and for finite angles of attack, ¢/sin? 6, is of the order of one. For 
example, for@, = 10°, e/sin? 6, = 1 for a = 1.75°; therefore, for 
finite angles of attack, this term must be retained. 
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Visualization of Boundary-Layer Flow 


James S. Murphy and Ralph E. Phinney 

Research Associates, Engineering Research Institute, University of 
Michigan 

August 30, 1951 


i FOR VISUALIZATION of boundary-layer flow has 
been developed at the University of Michigan Vacuum Type 
Supersonic Wind Tunnel. The technique was developed during 
experimental investigations on location of boundary-layer transi- 
tion! by the ‘‘China-clay’’? and “‘liquid-film”* 4 techniques and 
may be designated as a “China-film”’ technique, since it is a hybrid 
of these methods. A wide variety of applications of the technique 
is evident. It has been used in shock-reflection studies* to give 
indications of boundary-layer streamline directions and informa- 
tion about the line of impingement of a shock wave on a solid sur- 
face. In addition, proper application of the technique will yield 
useful information concerning location of boundary-layer transi- 
tion, regions of separated flow, tunnel blocking, interference 
effects, etc. 

The technique has its most fruitful application where three- 
dimensional boundary-layer effects near a surface are large, such 
as in the neighborhood of a three-dimensional shock reflection 
The position of shock reflection from a surtace 
first, through the sharp change 


from a surface. 
may become visible in two ways: 
in streamline direction (unless the reflection is strictly two-dimen 
sional, in which case there is no streamline curvature) or second, 
through the existence of a distinct streak under the shock (pre- 
sumably due to sharp shear-stress change under the shock foot 
or a separation bubble). For indicating streamline direction, 
this method has the advantage over tufts in that there is no flutter 
problem and a permanent pattern remains on the surface which 
may be examined and photographed after a run. 

Examples of the patterns obtained by use of the method are 
Fig. 1 indicates the streamline pat- 
In Fig 


shown in Figs. 1, 2, and 3. 
tern on the tunnel floor associated with tunnel blocking 
2, the square body has been rotated approximately 45° from its 
position during the run. The two lines on the wedge are the re- 
reflection of the wedge shock from the body (where it splits into 
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two waves). Fig. 3 shows the indication of boundary-layer tran- 
sition on a flat plate. Oblique lighting was used to obtain this pho 
tograph. The others were taken with direct light. 

The method requires the preparation of a suspension of China 
clay particles in a bonding agent and a volatile liquid that has 
the proper drying time. A more precise description of the mix- 
ture will be given below. Shortly after preparation, a layer ot the 
A free 
application by hand to the surface, such that the surface feels 
oily but the liquid does not run, yields a layer giving a good flow 
pattern that does not interfere noticeably with the external flow. 


suspension is applied to the surface to be investigated. 


After application of the suspension, flow is established over the 
surface. During the initial period of flow, the volatile liquid is 


carried along the surface by the action of shear stresses in the 





KEY: Dork ( wet) surface shows extent of laminar layer 
White ( dry) surface indicates turbulent layer 
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boundary layer. As time progresses, the volatile liquid evapo 
rates, leaving the china-clay particles attached to the surface by 
the bonding agent. These particles are arranged in streaks hay- 
ing approximately the direction of the streamlines in the bound- 
ary layer near the surface. In order to provide for a good photo- 
graph, it is desirable that the surface be black (see Figs. 1, 2, and 
3) so as to contrast with the white China-clay particles. For stee} 
surfaces, treatment with a chemical gun blueing has proved satis 
factory. 

The suspension may be prepared in the following manner. 
The bonding agent is attached to the china-clay particles by 
making a lacquer similar to that in reference 2 or 5. Mixture of 
this lacquer with a volatile liquid in the ratio of one part lacquer 
to five parts volatile liquid by volume will give a sufficiently 
dense suspension of particles. The mixture should not be pr 
pared too far in advance of use, since the China-clay particles 
tend to coagulate and form a sludge that causes the streamlines 
to be less distinct. 

Choice of volatile liquid is dependent upon the drying time re- 
quired in a particular application. Reference 2 gives a list of 
suitable liquids covering a large range of drying times. For a 
vacuum type supersonic tunnel having an atmosphere stagna- 
tion pressure, methyl salicylate and iso-safrole and mixtures of 
these liquids have been used successfully. For example, at J/ = 
1.90, a mixture of one part iso-safrole and three parts methy] sali- 
cylate by volume gives a drying time of approximately 30 sec 
in regions of laminar flow and 12 sec. in regions of turbulent flow 
In such a tunnel, the effect of transient starting conditions on 
flow of the volatile liquid is negligible, since the transients exist 
for less than 2 sec. while the steady flowing period is of the order 
of 15 sec. For continuous tunnels, the effect of transients is more 
pronounced, and, consequently, a volatile liquid having a longer 
drying time must be used.? 
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Local Compressible Pressure Coefficient 


V. G. Szebehely 
Associate Professor of Research, Virginia Polytechnic Institute 
July 23, 1951 


T° A RECENT NOTE, Laitone! suggests the use of the local Mach 
Number for computing compressible pressure coefficient and ob- 


tains his formula [Eq. (2)] by ‘‘analytical continuation.’’ Based 
on physical considerations, the relation between compressible and 


incompressible pressure coefficients, 


was derived by Truitt.” * Ina thesis‘ under the writer’s super- 
vision, the result was compared with experimental data for air- 
foils and with the von Karman-Tsien relation. It is interesting to 


note that the method does not offer always and everywhere on the 
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surface better results than the von Karman-Tsien formula. In 
fact, good agreement was found’ with experiments around the 
nose of the airfoil using the formula above. Between 20 and 70 
per cent chord points, however, on the upper surface of the 
NACA airfoil (4412) at 1/7. = 0.603 and a = —2°, the above 
formula overestimates the experimental values. Generally, the 
method appears to predict closer to the experimental values for 


pressure coefficients that are relatively large.® 
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The Effect of Compressibility of the Material on 
Plastic Buckling Stresses 


Elbridge Z. Stowell and Richard A. Pride 

National Advisory Committee for Aeronautics, Langley 
Air Force Base, Va. 

July 23, 1951 


I THE ORIGINAL N.A.C.A. computation of plastic buckling 
stresses,! they were taken to be of the form 
cr = n} wkRE/(12(1 — v?)]} (t/b)? 

in which v is the value of Poisson’s ratio which applies in the 
elastic range of the material and 7 is the plasticity coefficient ap- 
plicable to the problem in hand. The 7-coefficients were calcu 
lated by dividing the plastic buckling stress by the buckling stress 
that would be obtained if the material remained elastic, with 
Poisson’s ratio taken equal to 1/2 in both stresses for this pur 
pose. In other words, the material was taken to be incompres- 
sible as far as the computation of 7 was concerned, with the re- 
sult that a small error was introduced into the coefficient 7. 

In order to obtain an estimate of the magnitude of this error 
in the n-coefficients, the calculation of plastic buckling stresses 
has been repeated with allowance made for compressibility of the 


material. The stress-strain relations used were those of Hencky: 


o, = [Esec/(1 — p*)]| (er + ne,) 
o, = [Exe-/(1 — p?)] (e, + we) 
tr = Eeecy/(2(1 + pw)] 
where uz is a composite Poisson’s ratio of magnitude 
zw = (1/2) (1/2) — v] (Esec/E) 
The resulting expressions for plastic buckling stresses are: for 
the simply-supported flange 
Ger = | W*REvec/[12(1 — w?)]}(t/b)? 


for the simply-supported plate (with k = 4) 


4.34. Evan | Ecce [t \? 
wales: Sad sy + = we laa — =f: 
The trial and error process now requires, in the case of the flange, 
the use of Exec and yw and, in the case of the plate, the three quan- 
tities Exec, Etan, and wu, corresponding to each assumed oer until 
the respective equations are satisfied. 
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The results of such calculations for the 14S-T6 aluminum alloy 
used in the plastic buckling tests? are shown in the figure. The 
effect of compressibility is to raise the buckling stress for the 
flange by about '/, k.s.i. in 60 and to lower the buckling stress for 
the plate by about the same amount. Thus, the error in using 
the original formulas for 7 is of the order of 1 per cent for this 
material. 

There is no error in using Fis,/F as the n-coeflicient for a long 


column, 
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An Integral Relation from the Turbulent 
Energy Equation 


Stanley Corrsin 
Department of Aeronautics, The Johns Hopkins University 
August 14, 1951 


} maton VON KARMAN,! integral relations for mean motion 
have found wide application in both laminar and turbulent 
shear flows of ‘‘boundary layer type.’’ The purpose of this note 
is to demonstrate an integral relation deducible from the ‘“‘bound- 
ary-layer form” of the turbulent energy equation 

For steady state with two-dimensional mean motion,’ 


1 — dg? , 1 — dg? ou oy fl, P\l 
i $- fo -+fg— ~ —8-¢4 , 
2 Ox 2 oF ov oy | « ep} | 


v 0°? — (25> ) 
2 oy? Ox Ox 
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where x is the principal flow direction (U, ~); 
turbulent velocity (= u? + v? + w?); pis static pressure fluctu- 
ation; p is density; v is kinematic viscosity; and the last term 
is in Cartesian tensor notation, with repeated index indicating 
summation. The bar denotes average 


g’ is the resultant 


Integration of (1) in the y-direction across a shear flow, be- 
tween boundaries a(x) and h(x), gives 


1d b(x) . b(x) _ ou 
~ Ug? dy = — uv — dy — 
2 dx a(x F a(x) y 
b(x) 
Ou; Ou; 
» = "\dy (2) 
a(x) OX, Ox; 


The first term is x-rate of increase in rate of flow of energy; the 
last two are total production and dissipation, respectively. Con- 
vection disappears, since it is zero at axes of symmetry and/or at 
edge boundaries, the latter by ordinary choice of boundary. 

For sufficiently large Reynolds Numbers, dissipation occurs 
chiefly in the locally isotropic part of the spectrum, and the last 
integrand can be approximated by Taylor’s isotropic dissipation 
expression ® 


(Ou; /Ox,) (du; /dx) = 5q?/d? 


where A is the “‘microscale.”’ 

As an example, consider the application of Eq. (2) to predic- 
tion of \ in a “‘half-jet’’ of the type studied in reference 4, with 
the additional assumption of similarity (with x) in the turbulent 
as well as the mean motion and the empirical fact that \ is inde- 
pendent of y in the fully turbulent region. 
have used mean motion integral relatiois and mean flow simi- 


Liepmann and Laufer 


larity to demonstrate the linear spread of this particular turbu- 
With the 


lent shear flow, a ~ b ~ x. similarity statements 


[U = Umaz.:f(n); q? = 5 oe 2.9(m) and uv = Umaz.?- h(n), where 
n~y/(b a) ~ y/x], and with local isotropy, Eq. (2) gives 

d = [(v/A Umar.) +x]'/* (3) 
where A is a numerical constant, dependent upon various in- 


tegral functions of f, g, and h. 

This parabolic variation of \ is in rough agreement with the 
experimental results given in Fig. 35 of reference 4. 

A rough estimate by von Karman? from the differential form 
of the turbulent energy equation gives qo’A? ~ vé, where 6 is a 
length proportional to the shear region width and qo’ is a measure 


of the root mean square turbulent velocity. For the particular 


const., so that \ ~ x'/? 


case of the “‘half-jet,”’ 6 ~ x and qo’ = 
as pointed out by Liepmann and Laufer. 

In fact, von Karman’s process of dimensional reasoning and 
estimation in reference 2 is qualitatively the same as using an in- 
tegral relation plus similarity assumptions 
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The second example is an estimate of the change in average 
turbulence level in fully developed turbulent channel flow as q 
function of Reynolds Number. With a(x) and 6(x) equal to 
+H (the channel boundaries) and the replacement of xp in Eq. 
(2) by its expression in the momentum equation [Eq. (2) in ref. 
erence 5], there results 
H @ lop (#4 H 


a Udy — » 


g dU 2 
_H d? p OX J_y 


—H \dy 


ov dy 


(4) 


which says that the total energy dissipated to heat by the turby- 
lence equals the total rate at which energy is fed into the flow 
by pressure work minus the energy dissipated directly without 
passing through the turbulent state. 
of the average turbulent dissipation across the channel follows 
directly: 


A dimensionless measure 


0H [" @ 
dy = 1;-R-e, — I, (5 
Umer.? J—-H 
where 
1 
= . f(n) dn 
1 
I, = f , (df/dn)? dn 
f = U/ Umax. 
n = 9/H 
R= Umax H/yv 
270 2H(O0P/Ox) 
i” TA tal. 
pl mar PU maz. 
7 =u dU/dy, 


y= + H 


For an estimate of the R-variation in 


H 
E = (2U mas.2H)~ r q* dy 


consider \ roughly constant across the channel® and estimate 
c,(R) from the Blasius calculation on data in tubes;® this gives 
cy, = 0.05R~ “4, 
sublayer (of thickness e), J; 


With a (1/7)-power turbulent region and linear 
1.75 and I, + 2(H/e)’? =~ 0.1R%, 
and Eq. (5) becomes 


E ~ R°(\/H)*{R'* — 2.0) (6 


With the von Karman estimate mentioned earlier, (H/A)? ~ R 
Therefore, 


E ~ R~-*/(R'/* — 2.0) (7) 


which should be roughly true for extremely large R. Note that 
E = Ofor R = 260, probably the same order as the lower critical 
R for channel flow. 
Reynolds energy method for estimating this limit.’ 


In fact this example partly resembles the 
Fig. 1, with 
ordinate normalized at R = 30,800, gives a comparison of Eq 


(7) with the Laufer data.’ 
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Concerning Hadamard’s Solution for a Source 
in Supersonic Flow 


E, V. Laitone 
Associate Professor, University of California, Berkeley 


August 8, 1951 


N ANSWER TO SERBIN’S COMMENTS on the apparent sign discrep- 
I ancy in the analogy of ‘‘The Source-Flow in Supersonic and 
Incompressible Theories,’’! it should be pointed out that in the 
yon K&arman-Moore integral the pseudo-source 
strength is proportional to S”, where S = rr? is the cross-sectional 
area of the body of revolution, while in the Hadamard solution 
the supersonic source strength is proportional to S’, analogous to 


supersonic 


that for incompressible flow. 
This is shown in reference 2 wherein Eq. (17) gives the von Kar- 
man-Moore integral for a slender sharp nose body of revolution as 


?, ] ex—Br S"(é) (& — x) dt 
5 er i, 7 (1) 
l 2rr V/ (x — t)? — Br? 
and the corresponding Hadamard solution as 
d, B2r a Br |S'(x — Br) — S'(&)] dé x S'(x — Br) 
Un J0 [(x — &)? — B*r?]*?? 2ar ~/x2 — Br? 
(2) 


Note that Eq. (2) can be obtained either by integrating Eq. (1) 
by parts or directly from the Hadamard method as followed in 
Eq. (8) of reference 2. 

The corresponding solution for linearized subsonic flow is given 
by Eq. (51) of reference 2 as 
S'(é) dé 


[(x — £)? + B%r?]*/? 


¢, 


Br fr 
U 4x Jo 


It is seen that only Hadamard’s solution is analogous to the lin- 


(3) 


earized subsonic or incompressible flow (8 = 1) solution, as would 
be expected. 
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On the Unsteady Motion of a Thin Rectangular 
Wing in Supersonic Flow 


N. Rott 
Cornell University, Ithaca, N.Y. 
August 15, 1951 ; 


geo aaaagted SUPERSONIC FLOW PROBLEMS of the rectangular 
wing have been solved by Miles,” ? Stewartson,*® and re- 
cently (following an idea of Gardner‘) by Goodman.® * The 
following treatment is restricted to the case when the prescribed 
normal velocity yg, on the wing is an arbitrary function of the 
chordwise ordinate x and of time ¢ but is independent of the 
spanwise ordinate y. A general solution for this case will be 
presented, based on the potential in the corresponding two- 
dimensional problem. 

Consider a quarter-infinite plate occupying the plane x > 0, 
y > 0. After introduction of polar coordinates r, 3 in the yz- 


plane, the potential satisfies the equation 





FORUM 
1 1 aii 2M 1 
Ce t— & + — C00 = GP = Tess + @x: + eu (1) 
r r? a a* 
The solution is obtained following a method by Lamb” ° in 


his treatment of Sommerfeld’s diffraction problem. With a slight 
generalization of Lamb’s procedure, the following expression will 


be tried for the spanwise component ¢, of the velocity: 


¢y = cos (8 2)r—'/? F(x, t, r) (2) 
On the plate—i.e., for ) = 0 and 3d = 27—this yields g,, = 0, so 


that ¢y, is a function of x and ¢ only, in accordance with the 
boundary conditions stated above. Introducing Eq. (2) in Eq 
(1), it follows that F satisfies 


F., = (M? — 1)F,, + (3) 


i.e., the potential equation in two spatial dimensions x and r. 
The spanwise flow will vanish on and outside of the Mach cone 


ax 


thus 


F(x, t, ax) = 0 (4) 
Notice that for r—> 0, g, must have a singularity of the type r~'/2, 
so that F(x, ¢t, 0) must be finite. 

The potential on the plate follows by integration of Eq. (2): 


"y my, yo « 
J vg, dy = F(x, t, y) r (5) 
— ax 0 / 


Vy 
since g, = 0 for 3d = x. 

It will be necessary also to express the potential ¢(x, /, 2) in 
the two-dimensional domain by integrating Eq. (2) along a line 
z = const. between the limits given by the two intersections of 
this line with the Mach cone r = ax. This integration may be 
performed, following Lamb,” § by the substitution ./r cos (1/2)8 
* The result is 


g(x, t,y) = 


= yn, a method that is also described by Bateman 


- 2 : (6) 
¢(x, t, 2) = f F(x, t, 7) TT 
Z Vr-—-s 
This equation implies that @¢ = 0 for z = ax, on the limit of the 


disturbed region, a permissible normalization, as the limit is an 
equipotential surface of the disturbance potential. (If the dis- 
turbance velocity does not vanish at the limiting surface, then it 
is perpendicular to the latter by a general property of discon- 
tinuity surfaces. ) 

The functional relationship between go and ¢ as established by 
Eqs. (5) and (6) is sufficient to prove a significant property of the 
spanwise integral of go between y = 0 and y = ax. (Simple ex- 


pressions in special cases where obtained for this integral by 
Miles" 2 and Stewartson* independently.) From Eq. (5) it is 
found after integration by parts: 
ax ax dy 
go dy = ax Kz, t, 9) — 7 - 
0 0 Vy 
x " = 
 % F(x, t,y) Vydy (7) 
0 
The integral in the first term is 
go (x, t, ax) = (x, t, 0) 
For the second term, it may be shown that 
(8) 


ax J 1 ax 
f, F(x,t,v) Vydy = oe on t, s)dz 


for, substituting @ from Eq. (6) into the right-hand side integral, 
the following change in order of integration will be justified 
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l ax ax dr 

odo a: f F(x, t, r) /7 aa 
1 ax y Flx. t dz 
i. arf, (x, iO eee 


The domain covered by this double integral in the rz-plane is a 
triangle with the corners r = 0,2 = 0; r = ax,s = 0; andr = 
ax, 2 = ax. The original form of the integral on the left-hand 
side of Eq. (9) requires the first (inner) integration on strips 
parallel to the r-axis, while the right-hand form represents the 
same integral first evaluated on strips parallel to the z-direction 
Integration in the latter with respect to z gives 


l ax r dz ax 
2 f dr F(x, t, r) f, /r — f F(x, t,r) UV rdr 


which proves Eq. (8). Inserting these results in Eq. (7), one 


obtains 


ax 1 ax 
J vAx, t, y) dy = ax dx, t,0) — — S ox, t,2) dz (10) 


Evidently in this equation ‘potential’? may also be replaced by 
‘pressure’; then the statement of Eq. (10) is, in words: The lift 
of a spanwise strip in the region of the side-edge disturbance is 
equal to the lift of a corresponding strip in the two-dimensional 
domain minus one-half of the pressure integral perpendicular to 
the wing at the same chordwise station in the two-dimensional 
domain. 

The ultimate aim is to find a general connection between go and 
¢, if @ is given, since the two-dimensional problem may be con- 
For this purpose Eq. (6) should be inverted, 
For Lamb’s original treat- 


sidered as soluble. 
so that F is expressed explicitly by ¢. 
ment,’ Bateman’? has solved this problem by Abel’s inversion for- 


mula. In the present case, the solution is 


dz 


4 (11) 
VWr-s 


F(x, t, r) 


Insertion of Eq. (11) into Eq. (6) gives an identity after the order 

of integration is changed in the manner described above. 
Introducing, finally, Eq. (11) into Eq. (5) and changing the 

order of integration again, one has the final result (after integra- 


tion by parts) ; 


| 


1 ax | 
go x,t, vy) = O(x, 1,0) — &(x, t, 2) V 
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Eq. (12) permits the calculation of the potential on the wing in 
the domain of three-dimensional disturbance by the use of the 
two-dimensional solution, which must be known in its whole 
field of validity. This potential, ¢(x, ¢, z), can be determined by 
straightforward extension of the various methods customarily 
used to calculate the surface values ¢(x, ¢, 0) only. Applied to 
the cases treated previously in steady and unsteady motion, the 
results of Eq. (10) and (12) are in agreement with known solu- 


tions. [Actually, the possibility of a general relationship as in 


Eq. (12) was suggested to the author by inspection of the results 
| gs j 


of references | to 6. | 
It should be noted that the steady case permits a simpler 


treatment, using a general solution of Eq. (3), 

Fix,7r) = flax — r) 
If, moreover, a constant angle of attack is required, then ¢ must 
be conical, a condition that is fulfilled in Eq. (2) only by putting 


flax — r) = const. Vax —r 


The final result in this special case follows by elementary steps 
only. 

The simplicity of the approach outlined here is obviously based 
on the fact that, by the use of Lamb’s form, Eq. (2), the depend- 
ence of g, from # is explicitly given by a factor cos (1/2)3. It 
should be mentioned that the boundary conditions would be 
satisfied also by the more general factor cos [m + (1/2) ]#, but, if 
this is used, the differential equation for F will show that the only 
cases where F may be finite for r = O aren = Oandn = —1, both 
giving Lamb’s form. Only regular and nonvanishing functions 
F(x, t, 0) will give the required singularity of ¢, for r— 0. 
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